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Abstract 



In a recent series of papers |I]-|1] a novel method based on the cod- 
ing of integer partitions has been used to derive power series expan- 
ds) ■ sions to previously intractable problems, where the standard Taylor/- 

Maclaurin method fails. In this method the coefficient at each order 
k of the resulting power series expansion is determined by summing 
k>( ' all the specific contributions made by each partition whose parts/ele- 

^ . ments sum to k. The specific contributions are evaluated by assigning 

values to each element in a partition and then multiplying by a multi- 
nomial factor, which depends on the frequencies of the elements in the 
partitions. This work aims to present for the first time the theoretical 
framework behind the method, which is now known as the partition 
method for a power series expansion. To overcome the complexity in 
evaluating all the contributions from the partitions as the order k in- 
creases, a programming methodology is created, thereby allowing for 
far more general problems to be considered than originally envisaged. 
This methodology is based on an algorithm called the bi-variate recur- 
sive central partition (BRCP) algorithm, which is, in turn, developed 
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from a novel non-binary tree-diagram approach to scanning the inte- 
ger partitions summing to a specific value. The main advantage of the 
BRCP algorithm over other means of generating partitions lies in the 
fact that the partitions are generated naturally in the multiplicity rep- 
resentation. By developing the theoretical framework for the partition 
method for a power series expansion it becomes apparent that scan- 
ning over all partitions summing to a particular value can be regarded 
as a discrete operation denoted here by the discrete operator Lp^k[-]. 
The summand inside this operator depends on the coefficients of an 
inner and outer series resulting from expressing the original function 
as a pseudo-composite function of two power series expansions. As 
a consequence, simple modifications to the program for the partition 
operator result in programs for other operators involving specific types 
of partitions such as those with: (1) only odd or even elements, (2) a 
fixed number of elements, (3) discrete elements, (4) specific elements 
and (5) those with restrictions on the size of their elements. Another 
interesting modification results in the generation of all the conjugate 
partitions for their original partitions by transposing the rows and 
columns of their Ferrers diagrams. The operator approach is then ap- 
plied to the theory of integer partitions, in particular to generalisations 
of the generating functions for both discrete and standard partitions. 
The main generalisation involves the introduction of the parameter 
u!, whose powers indicate the total number of elements in the par- 
titions, while the coefficients of the power series expansions become 
polynomials in lo. Finally, power series expansions for more advanced 
infinite products involving quotients and products of the discrete and 
standard partition generating functions are derived, culminating in 
the multi-parameter infinite product first studied by Heine. 
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1 Introduction 

This work grew out of a desire to develop a programming methodology on 
the partition method for a power series expansion, which has featured promi- 
nently in a recent series of papers [1]-[1] aimed at deriving power series ex- 
pansions for previously intractable mathematical functions. Although still a 
relatively novel method, the method for a power series expansion was first 
introduced in the derivation of an asymptotic expansion for the particular 
Kummer or confluent hypergeometric function that emerges in the response 
theory of magnetised quantum plasmas such as the degenerate electron and 
charged Bose gases [3]. Specifically, a large | a [-expansion was obtained for 
i-Fi(a, a + 1; z), which is itself a variant of the incomplete gamma function 
7(0;, z). As a consequence, the physical properties of the magnetised charged 
Bose gas were later determined in the weak field limit in Ref. [6], the first 
time ever that the limit had been investigated for a plasma. Since then, 
the method for a power series expansion has been applied successfully to 
various mathematical functions such as l/ln'*(l + z), sec* z, z'^csc'^ z, and 
the three Legendre-Jacobi elliptic integrals, F{ip, k), E{ijj, k) and Il{ip, n, k), 
while more recently, it has been applied to a finite sum of inverse powers of 
cosines |1], resulting in the development of a spectacular empirical method 
to solve the problem. 

Initially, it was observed that in order to apply the partition method for a 
power series expansion the original function needed to be expressed as a com- 
posite function. In the present work, however, this condition will be relaxed 
to quotients of "pseudo-composite" functions. The two functions involved 
in the construction of the quotients of the pseudo-composite functions must 
themselves be expressible in terms of power series expansions, which are re- 
ferred to as the inner and outer series. Neither of these series, however, is 
required to be absolutely convergent. When the conditions for the quotients 
are met, a power series expansion can be obtained in which the coefficients 
at each order k are calculated by summing the specific contributions due to 
each partition in which the parts or elements sum to k. The contribution 
made by each partition is in turn determined by: (1) assigning a value to 
each element in the partition, (2) multiplying these values by a multinomial 
factor composed of the factorial of the total number of elements, N^l divided 
by the factorial of the number of occurrences or frequencies of each element 
i in the partition, nj, and (3) if necessary, carrying out a further multiplica- 
tion with the coefficients of the inner series set at the number of elements in 



the partition. 

As explained in the introduction to Ref. [3], the partition method for 
a power series expansion is able to produce power series where the stan- 
dard technique or Taylor/Maclaurin series approach breaks down, but for 
those cases where a power series expansion can be obtained by the stan- 
dard technique, the power series expansions are identical, although from a 
totally different perspective. As a consequence, the cross-fertilisation of both 
approaches is frequently responsible for new mathematical results and prop- 
erties. A particularly fascinating property of the partition method for a 
power series expansion is that the discrete mathematics of partitions is being 
employed to derive power series expansions for continuous functions. 

Whilst Ref. |3] was primarily concerned with the application of the parti- 
tion method for a power series expansion to basic trigonometric and related 
functions, it was stated that a theoretical framework was being developed 
for situations involving more complicated functions where the assignment 
of values to the elements of the partitions is no longer specific, but quite 
general. This theoretical framework was necessary for the the development 
of a programming methodology whose purpose was to facilitate the summa- 
tion process over the partitions. Therefore, both the theoretical framework 
and development of a programming methodology represent important topics 
in the present work. However, in the course of developing the theoretical 
framework it became apparent that the theory of partitions was also af- 
fected. Moreover, with the development of the new algorithm to facilitate 
the calculation of the coefficients in the partition method for a power series 
expansion, one could tackle various problems pertaining to integer partitions 
such as the evaluation of: (1) doubly-restricted partitions, (2) partitions 
with a fixed number of elements, (3) conjugate partitions via Ferrers dia- 
grams, and (4) discrete/distinct partitions. As we shall see, such problems, 
which often require implementing different algorithms if they can be solved 
at all, can be addressed by making minor adjustments or modifications to 
the algorithm presented in Sec. 2 of this work. Therefore, by developing the 
theoretical framework and tackling the mathematical programming issues 
associated with the method, we have not only been to improve the method 
by making it more general than first thought, but we have also been able 
to solve other outstanding problems in the theory of partitions. In addition, 
the new operator approach that evolved in the process casts partitions in a 
different light and it is hoped that this approach will result in new advances 
in the future. Consequently, the work has progressed well and truly beyond 



its original conception resulting in its present size. 

Throughout this work the concept of regular isat ion of a divergent series 
will be employed. This concept is defined here as the removal of the infinity 
in the remainder so as to make the series summable or yield a finite limit. 
The finite limits obtained in this process are referred to as regularised values, 
whilst the statements in which they appear together with the series can no 
longer be regarded as equations. Instead, they are referred to as equivalence 
statements or equivalences, for short. A necessary property of the regularised 
value is that it must be unique, particularly if it is identical to the value one 
obtains when the series is absolutely convergent within a finite radius of 
the complex plane such as the geometric series. Much of this will become 
clearer as we progress further in the present work, but for those readers 
wishing to seek a greater understanding of the concept and its ramifications 
to asymptotics, they should consult Refs. |3] and j7]-pT]. 

Since all the partitions summing to a particular value are required to 
determine the coefficients in the partition method for a power series expan- 
sion. Sec. 2 examines the current state of the art for generating partitions. 
Here, various algorithms are compared with the bi-variate recursive central 
partition algorithm or BRCP for short, which has only been sketched out in 
Refs. [1] and [3] . The BRCP algorithm is based on a non-binary tree diagram 
representation for all the partitions summing to a specific value. Whilst it is 
acknowledged that some of the alternative algorithms can be faster than the 
BRCP algorithm in scanning the partitions and occasionally, when printing 
them out on a screen, the BRCP algorithm is the most efficient algorithm 
for implementation in the partition method for a power series expansion. 
This is because it generates the partitions in the multiplicity representation, 
which, as we shall see, is both the ideal and minimum amount of information 
required by the partition method for a power series expansion. 

Sec. 3 presents the general mathematical theory underpinning the par- 
tition method for a power series expansion, which is discussed in terms of 
the quotient of two pseudo-composite functions in Theorem 1. Under certain 
conditions it is also shown that the method can be used to derive a power 
series expansion for the inverse or reciprocal of the quotient. Consequently, 
a hybrid recurrence relation involving the coefficients from both power series 
expansions is obtained, while an example involving the reciprocal of a Bessel 
function to arbitrary order, viz. 1/ Ji,{z), is presented to make the preceding 
material clearer to the reader. Later in the section. Theorem 1 is extended 
in a corollary, where the quotient of the pseudo-composite functions is taken 
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to an arbitrary power, p. This extension only affects the multinomial fac- 
tor in the method for a power series expansion by transforming it into the 
Pochhammer symbol involving p and the number of elements in each parti- 
tion. 

Because the partition method for a power series expansion is based on 
summing over all the partitions summing to a specific value /c, a discrete 
operator called the partition operator or -^p,fc[-] is introduced. When this op- 
erator acts on unity, it gives the number of partitions or partition function, 
p{k). On the other hand, if the operator acts upon only the multinomial fac- 
tor mentioned above in Step 2, then it is found to yield a value of 2^~^, while 
if the phase factor of (— 1)^*^, where Nj. represents the number of elements, 
is included with the multinomial factor, then it vanishes. This means that 
the expressions for the coefficients given in Theorem 1 can be expressed in 
terms of the new operator. One interesting property of the partition opera- 
tor is that frequently the arguments inside the operator can be interchanged 
with the result the operator yields. As a consequence, not only are many 
of the results derived previously via the partition method for a power series 
expansion in Refs. [I]- [3] expressed in terms of the operator, but also the 
inverse relations are presented. The section concludes by showing that if the 
quotient of the pseudo-composite functions yields an infinitely differentiable 
function, then its derivatives can also be expressed in terms of the partition 
operator. 

Sec. 4 is devoted towards the creation of the programming methodology 
that enables the coefficients obtained from the partition method for a power 
series expansion to be calculated from the general theory in Sec. 3. This 
is necessary because beyond the first few orders, it becomes increasingly 
onerous to determine the coefficients by hand. The computational task is 
divided into two steps. In the first step a code utilising the BRCP algorithm 
of Sec. 2 is written in C/C++ so that its symbohc output can be processed 
by Mathematica |12J in the second step. The reason for the second step 
is that often the coefficients in the resulting power series expansions are 
either rational or algebraic in nature and that neither of these forms can 
be handled effectively in C/C++ with its floating point arithmetic. In fact, 
the coefficients often become so small that they would practically vanish in 
C/C++. By importing the output into Mathematica, we can exploit its 
integer arithmetic routines to express the coefficients in integer form or we 
can use the symbolic routines to express them, for example, as polynomials 
when required. 



The material in Sees. 2 to 4 serves as a platform for studying various 
problems in theory of partitions. First, we consider the issue of generating 
specific types/classes of partitions or different operators such as: (1) those 
with a fixed number of elements in them, (2) doubly restricted partitions 
where all the elements are greater one value and less than another, (3) dis- 
crete or distinct partitions where an element occurs only once in a partition 
and (4) partitions with specific elements in them. Solving these problems 
invariably means developing new algorithms or codes, but in the case of the 
BRCP algorithm we shall see that they can be solved with relatively minor 
modifications. This is due to the power and versatility of the of the tree 
diagram approach upon which the BRCP algorithm is based. Furthermore, 
the BRCP code can be adapted to exploit two-dimensional dynamic memory 
allocation when one wishes to determine conjugate partitions by means of 
Ferrers diagrams. For the benefit of the reader many of the codes discussed 
in Sees. 4 and 5 are presented in the appendix, where it can be seen that 
they are surprisingly compact. 

In Sees. 6-8 the partition method for a power series expansion and the 
operator approach are employed in the derivation of generating functions 
from increasingly sophisticated extensions of the infinite product defined by 
-f (^) = rifcLill ~ ^'^) ^^d its inverse or 1/P{z). This famous product was 
found by Euler to yield the generating function whose coefficients are equal 
to the partition function p{k). Before the study commences, however. Sec. 6 
begins by showing in Theorem 2 that the generating function of P{z) is ab- 
solutely convergent within the unit disk centred at the origin in the complex 
plane, but is divergent for all other values of z. That is, P{z) represents the 
regularised value for divergent values of z as given by Equivalence fl73l) . This 
is often postulated in the literature, but no formal proof of this important 
result has ever appeared previously. Next Theorem 1 is used to derive the 
generating function for 1/ P{z) whose coefficients are given in terms of the 
partition operator acting with each element i assigned a value of p{i). In 
this case the coefficients q{k)^ which are referred to as the discrete partition 
numbers, are only non-zero when /c is a pentagonal number, again a result 
that was first obtained by Euler. By inverting this method we obtain the 
partition function p{k) in terms of the partition operator acting with each 
element i assigned to q{i). 

Because most of the discrete partition numbers vanish, the new result for 
the partition function simplifies dramatically when a program based on Sees. 
4 and 5 is created. In fact, scanning over those partitions in which the ele- 
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ments are pentagonal numbers represents a totally different application from 
the examples studied in Sec. 5. Therefore, a code is developed which only 
determines the partitions whose elements are pentagonal numbers. However, 
if P{z) and its inverse are expressed as exponentiated double sums, then 
Theorem 1 can be used to derive different expressions for the coefficients of 
the generating functions. In this instance the elements in the partitions are 
assigned values 7j, which are obtained from summing the divisors or factors 
of i divided by i. As a result, it is found that the only difference between 
the discrete partition numbers and the partition function via this alternative 
approach is that the former set of numbers possess an extra phase factor of 
( — 1)^*= inside the partition operator, where again N^ is the total number of 
elements in a partition. 

Sec. 7 begins with an extension of the inverse of P{z), where the coeffi- 
cients of z^ in the product are now equal to Ck rather than -1. Theorem 3 
shows that a generating function can be obtained from this product where 
the coefficients are determined by applying the discrete partition operator 
LDP,k['] acting with each element i assigned a value of Cj. Conversely, this 
theorem implies that any power series can be expressed as an infinite prod- 
uct. If the Ck equal the parameter u, then the coefficients of the generating 
function for the product Q{z, u) not only become polynomials of degree k in 
u giving rise to the discrete partition polynomials q{k, u), but also the powers 
of u yield the number of elements involved in the discrete partitions. Then 
more identities involving the discrete partition polynomials are derived be- 
fore a corollary to Theorem 3 appears, the latter dealing with the case where 
the factors in the product are taken to an arbitrary power pk- Hence, one 
can derive generating functions for very complicated products with varying 
powers and/or with different factors accompanying the powers of z. 

Sec. 7 concludes by looking at the special case in the corollary to Theo- 
rem 3 where the arbitrary powers pk and coefficents Ck are set equal to the 
constant value p and the parameter u, respectively. In this case the coeffi- 
cients of the generating function become q{k,u,p) and are now polynomials 
of degree k in both p and u. Since the special cases of p = 2 and 3 feature in 
well-known products studied by Euler and Gauss, the properties and values 
of q{k, u, 2) and q{k, u, 3) are also examined in detail. 

Sec. 8 deals with the derivation of the generating functions for even more 
complicated products than those appearing in Sec. 7. As a result of the suc- 
cess of introducing the parameter u in the inverse of P{z), the first example 
deals with the introduction of u into P{z). This results in a generating func- 
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tion for the new product P{z, u) whose coefficients p{k, u) are polynomials 
of degree k in u that reduce to the partition function p{k) when u is set 
equal to unity. These partition function polynomials are expressed in terms 
of the partition operator acting with each element i assigned to the discrete 
partition polynomials or rather q{i,—u). They are found to possess many 
interesting properties, while their coefficients represent the number of parti- 
tions in which the number of elements is given by the power of u. Out of 
this analysis interesting recurrence relations are obtained for the number of 
discrete partitions or q{i,l). 

As was the case for P{k), P{k, u) and its inverse can also be expressed as 
an exponentiated double sum, both of which can be handled by Theorem 1. 
Thus, it is found that the partition function and discrete partition polyno- 
mials can be be expressed in terms of the partition operator with the main 
difference being that in the case of the former polynomials the elements i are 
assigned to the polynomials 7j(co'), while for the latter they are assigned to 
—'~fi{u). These new polynomials represent the extension of the 7j mentioned 
above. Their coefficients are equal to divisors d oi i divided by i, while each 
power of a; is equal to the reciprocal of the coefficient. Moreover, they reduce 
to the 7j when u = l. Because this is a somewhat unusual situation involving 
divisors, a program is presented that evaluates the partition function and dis- 
crete partition polynomials in symbolic form so that they can be processed 
by Mathematica. This means that the final forms for the polynomials can be 
obtained by evaluating the divisor polynomials via the Divisors [k] routine in 
Mathematica rather than having to create a separate program to solve this 
problem. 

Sec. 8 continues with the introduction of an arbitrary power p into the 
product P{z,u) and determining the coefficients p{k,u,p) of the generating 
function for this extended product. These coefficients are given in terms of 
the partition operator acting with the elements i assigned to minus the parti- 
tion function polynomials, viz. —p{i, u), multiplied by the Pochhammer sym- 
bol of {—p)Nk- Next the generating function of the product of Q{z, —(iuj) with 
P{z, aoj) is studied. This infinite product denoted by P{z, f3u, au) combines 
the properties of discrete partitions with standard partitions. The coefficients 
of the resulting generating function, which are denoted by QPk{u, a, /3), are 
polynomials of degree k in u, while the powers of a and /3 indicate the num- 
ber of elements in the standard and discrete partitions respectively. With 
this result Sec. 8 concludes with the derivation of the generating function for 
Heine's product, which can be represented as the product of P{z, ux, u) and 
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P{z,uy,u)xy). The coefficients of the generating function for this infinite 
product, which arises in q-hypergeonietric function theory, are denoted by 
HPk{(jJ, X, y) and are obtained by summing the product of QPj^oo, x, 1) with 
QPk-j{uj, y, xy) for j ranging from to k. Several coefficients in the last two 
examples are tabulated in order to display their complicated nature. 

2 Generating Partitions 

As indicated in the introduction the partition method for a power series ex- 
pansion is composed of two major steps. The ffist step involves determining 
all the partitions summing to an integer /c, which represents the order of 
the variable in the resulting power series expansion. The second and more 
complicated step is to calculate the contribution that each partition makes 
to the coefficient of the fc-th order term in the series expansion. This step 
will be described extensively in the next two sections when the theoretical 
framework and programming methodology for the partition method for a 
power series expansion are presented. For now, however, this section is de- 
voted to the problem of generating partitions in a suitable format to enable 
the second step of the partition method for a power series expansion to pro- 
ceed. This means that we shall not only be interested in determining the 
parts or elements in a partition, but also with evaluating their number of 
occurrences or frequencies. So, whilst the generation of partitions is an in- 
teresting problem in its own right and continues to be the source for new 
algorithms as evidenced by Refs. [I3]-[T6], it is required here in order to de- 
velop a programming methodology for the partition method for a power series 
expansion. Hence, we need to examine the existing algorithms for generating 
partitions to find which, if any, is the most suitable for implementation in the 
partition method for a power series expansion. Ultimately, we shall find that 
the novel algorithm sketched out in Ref. [1] will prove to be the most suitable. 
Moreover, we shall see in Sec. 5 that partitions with specific properties can be 
determined by making modifications to the algorithm, which is often difficult 
to achieve, if not impossible, with the other partition-generating algorithms. 
That is, a completely different algorithm is usually required to solve for each 
specific property of partitions, whereas only minor changes to the bi-variate 
recursive central partition algorithm presented in this section are needed to 
generate partitions with these specific properties. In addition, as a result of 
the material in appearing Sees. 3 and 4, we shall be able to re-formulate the 



11 



partition method for a power series expansion in terms of a partition opera- 
tor denoted by Lpfc[-]. The modifications to the BRCP algorithm presented 
in Sec. 5 will mean that we are effectively programming different operators, 
which will, in turn, lead to the presentation of new and fascinating results 
when we study the various generating functions belonging to the theory of 
partitions in Sec. 6. 

As described in Ref . [I] , when applying the partition method for a power 
series expansion there is actually no need to generate all the partitions at 
each order. For example, one can write down all the partitions necessary 
for evaluating the first few orders on a sheet of paper. Once they have been 
determined, one can then proceed to the second step of determining the 
specific contribution made by each partition to the coefficient of the series 
expansion. The problem occurs when we wish to evaluate the higher order 
terms, especially if our ultimate aim is to derive an extremely accurate ap- 
proximation to the original function. When we need to go to higher orders, 
the complexity increases dramatically due to the exponential increase in the 
number of partitions. Then it is no longer feasible to write down all the par- 
titions and carry out the calculations to determine their contributions to the 
coefficient at a particular order. Consequently, a programming methodology 
is required for all orders despite the fact that this will ultimately become 
very slow for very high orders of the series expansion due to combinatorial 
explosion. Nevertheless, we shall see that in developing this programming 
methodology we shall uncover very interesting results for the first time in 
the theory of partitions. For example, it has already been stated that series 
expansion obtained via the partition method produces a power series that 
is identical to a Taylor/Maclaurin series when the latter can be evaluated. 
In these situations the development of a programming methodology for the 
partition method means that the higher order derivatives in such a series can 
be expressed in terms of a sum of the contributions from all the partitions 
summing to a particular order. This has profound implications in mathe- 
matics in that we now have a means of linking the continuous/differentiable 
property of a function with the discrete mathematics of partitions or number 
theory. 

When it was mentioned above that there was no need to employ an al- 
gorithm to generate all the partitions in the partition method for a power 
series expansion, it was meant that there was no need to write them down 
in for each value of k, since only the elements in each partition and their 
frequencies are required as the input for the second step of the partition 
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method for a power series expansion. Representing a partition in this man- 
ner is known as the multiphcity representation, whereas we shall refer to 
the representation where each element in a partition is written down as the 
standard representation. As k increases, the number of partitions summing 
to k or p{k) increases exponentially. E.g., the number of partitions summing 
to 100, viz. j9(100), is 190 596 292. As a result, it is no longer practical to 
write the partitions in the standard representation. Whilst the multiplicity 
representation is sufficient for the application of the power series for a power 
series expansion, if one wishes to determine high orders of the resulting power 
series expansion via the method, then one still needs to consider the vari- 
ous algorithms for generating partitions because it could turn out that an 
algorithm generating partitions in the standard representation may require 
only minor modifications to provide them in the multiplicity representation. 
Furthermore, the generation of integer partitions continues to attract interest 
up to the present time. Therefore, we shall review the existing algorithms, 
but ultimately our aim will be to determine the most appropriate for the 
partition method for a power series expansion. On the other hand, those 
with only an interest in generating partitions may find the other algorithms 
more suitable in which case they are urged to obtain more information by 
consulting the list of references. 

Having justified the need for generating all the partitions summing to an 
arbitrary integer, we now turn to the issue of finding an appropriate algo- 
rithm that will expedite the process, but will do so in appropriate form for the 
second step of the partition method for a power series expansion. For a time 
there seemed to be only one useful algorithm for generating partitions. This 
was McKay's algorithm [17], which was basically a succession rule whereby 
partitions were generated in linear time. It was developed further by Knuth 
[18j, who used the fact that if the last element greater than unity is a two, 
then the next partition can be determined very quickly. This modification 
means that each partition takes almost a constant amount of time to be 
generated. The Knuth/McKay algorithm, which is implemented in C/C++ 
below, generates the partitions summed to a global integer n in a particular 
form of the standard representation known as reverse lexicographic order. 
Consequently, the elements in a partition, say a[l] up to a[k], are printed out 
according to a[l] > a [2] > ■ ■ ■ > a[k], while the first element of each new 
partition is less than or equal to the first element of the preceding partition. 
For example, the partitions summing to 5 in reverse lexicographic order are: 
5| 
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4|1 
3|2 
3|1|1| 
2|2|1| 
2|1|1|1| 
1|1|1|1|1| 

The rules for generating partitions in reverse lexicographic order can be ob- 
tained from Refs. [18] and [19]. Briefly, if the partition is not composed 
entirely of ones, then it ends with a value of x + 1 followed by zero or more 
ones. The next smallest partition in lexicographic order is obtained by re- 
placing the segment of the partition {. . . , x + 1, 1, ..., 1} by {■ ■ ■ ,x, ..., x, r}, 
where the remainder r is less than or equal to x. 

/* This program determines partitions in reverse 
lexico graphic order following McKay/ Knuth 
algorithm as discussed on p. 38 of Fascicle 3 
of Vol. 4 of D.E. Knuth ' s The Art of Computer 
Programming . */ 

T^include <stdio.h> 
T^include <memory . h> 
T^include <stdlib.h> 

int main(int argc , char *argv[]) 

{ 

int *a , i , m, n , q , x ; 

if (argc != 2) printf (" execution : ./knuth <partition7^>\n" 

else 

n=atoi ( argv [ 1 ] ) ; 

a=(int *) malloc ( (n + 1)* sizeof ( int ) ) ; 
PI: a[0]=0; 

m=l; 
P2 : a [m] = n ; 

q=m— (n= = l); 
P3: for(i=l; i<^m; i++) printf ("%i I" , a[i]); 

printf( " \n" ); 

if(a[q] != 2) goto P5 ; 
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P4: a[q— ] = 1; 

a[++m] = l; 

goto P3 ; 
P5 : if(q==0) goto end; 

x= a[q]-l; 

a[q]=x; 

n= m-q+1; 

ir^ q + 1; 
P6: if (n<= x) goto P2 

a [irir|-+]=x 



n — =x ; 
goto P6 ; 



end : ; 

} 
printf ("\n" ); 

free (a ) ; 

return (0) ; 

} 

By current standards the above code is considered to be slow for generat- 
ing the partitions at each order due to the excessive unconditional branching. 
From a computational point of view it is also very non-structured and hence, 
does not accord with modern programming practice. A significantly faster 
algorithm for generating partitions in reverse lexicographic order has been 
developed by Zoghbi and Stojmenovic in Ref. [13]. Actually, these authors 
present two algorithms in their paper, but the second, which generates the 
partitions in lexicographic order, is slower than the first. Nevertheless, if one 
runs the above code against a C/C++-coded version of the first algorithm, 
then one finds that it takes 1362 CPU seconds to print out the partitions 
summing to 80 on the screen of a Sony VAIO laptop with 2 GB RAM com- 
pared with 1399 CPU seconds using the Knuth/McKay code given above. 
On the other hand, if the partitions are directed to an output file, then it 
takes only 30 CPU seconds with the Knuth/McKay code compared with 28 
CPU seconds with the Zoghbi/Stojmenovic code. 

/* This program determines p artitions in ascending order 
following the algorithm given at J. Kelleher's web— page. *^ 
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T^include <stdio.h> 
:^include <memory . h> 
T^include <stdlib.h> 

int niain(int argc,char >Kargv[]) 

{ 

int *a, i, in, n, yduininy, xduminy, count; 

if(argc!= 2) printf (" execution : ./kellelier <partit ion7^>\n" 

else { 

n=atoi ( argv [ 1 ] ) ; 

a=(int *) malloc ( (n+2)*sizeof ( int ) ) ; 

for(i=0; i<=(n + l);i++) a[i]=0; 

a[l] = n; 

count =1; 

while (count != 0){ 

xduniniy=a [ count — 1] + 1; 
ydumniy=a [ count ] — 1 ; 
count=count —1; 
while (xdummy <= ydummy) { 
a [ count ]=xduniniy; 
yduniniy=ydummy— xdummy ; 
count=count +1; 

} 
a [ c o u n t ] =xdummy-|-ydummy ; 

for (i=0;i<= count ;i++){ 

printf("%i" ,a[i]); 

if (i<count) printf("|"); 

} 
printf ("\n" ); 

} 

} 
printf ("\n" ); 

free (a ) ; 

return (0) ; 

} 

Lexicographic ordering of the partitions ordering of the partitions is not 
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the only method of generating partitions in ascending order. Kelleher and 

O'Sulhvan in Refs. [13] and ^E\ have developed algorithms for generating 

partitions in ascending order, but not in lexicographic order, one of which 

has been created as the C/C++-code appearing above. When this code is 

run for the partitions summing to 5, the following output is obtained: 

1|1|1|1|1 

1|1|1|2 

1|1|3 

1|2|2 

1|4 

2|3 

5 

Now, by inverting the reverse lexicographic order of the partitions summing 

to 5 given earlier, we see that the partition {1,1,3} in Kelleher/O'SuUivan 

code appears before {1,2,2}, while in McKay/Knuth code it appears after 

the latter partition. A similar situation occurs with the {1,4} and {2, 3} 

partitions. Yet, the partitions for both codes are arranged in ascending order. 

Moreover, Kelleher and O' Sullivan are able to take advantage of the different 

ordering between to develop an even more efficient version of the above code. 

As a result, they state that their optimised version is much superior to the 

Zoghbi/Stojmenovic code. Interestingly, if the above code is run to yield the 

15 796 476 partitions summing to 80 on the same Sony laptop as the previous 

code, then it takes 1342 CPU seconds to output the partitions. Yet, if one 

does the same with the more efficient version of their code, then it takes 

the same amount of time to generate the partitions. On the other hand, if 

the partitions are directed to an output file, then it is found that the above 

code takes 29 CPU seconds, while the more efficient version takes 28 CPU 

seconds, the same time taken as the Zoghbi/Stojmenovic code. 

According to Ref. [20], Kelleher has found on his computer system that 
partitions summing to 80 are generated at a rate of 1.30 x 10^ per second using 
the first algorithm, while with the second algorithm the rate is 2.87 x 10^ per 
second. For the Zoghbi/Stojmenovic code the rate is 1.26 x 10® per second, 
while with the Knuth/McKay code the partitions are generated at a rate of 
1.73 X 10® per second. Hence, the reason why there was no marked difference 
in performance in the C versions of the Kelleher codes is attributed to the 
manner in which the partitions were printed out. 

Unfortunately, the above codes do not utilise the recursive nature of par- 
titions, which can be observed by realising that the partitions summing to 
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k+1 include all the partitions summing to k with an extra element of unity in 

them in addition to other partitions possessing elements greater than unity. 

In fact, according to p. 45 of Ref. [TH], the number of partitions summing to 

k 
k with m elements, which is denoted by 



relation: 



TTl 



, obeys the following recurrence 



k-1 
m — 1 



+ 



k — m 
m 



(1) 



In addition. 



k + m 
m 



represents the number of partitions summing to k with 



at most m elements. If P{k, m) represents the partitions summing to k with 
at most m parts, then with the aid of the above recurrence relation we obtain 



P{k, m) = P{k, m — 1) + P{k — m, m) 



(2) 

This result is derived on p. 96 of Ref. [2T] . 

To incorporate recursion into the process of generating the partitions, we 
need to construct an algorithm that utilises the special tree diagrams, which 
first appeared in Ref. |5j, but have since been applied to other problems 
or situations in Refs. [I]- [3]. These trees are different from the binary tree 
approach in the recent work of Yamanaka et al [16] , which seeks to generate 
each partition in the standard representation in a constant time rather than 
an average time. Their work represents a further development on Fenner 
and Loizou [22], who seemed to have been the first to develop a binary tree 
representation for partitions based on a lexicographic ordering. To construct 
the special tree diagrams, one begins by drawing branch lines to all pairs of 
numbers that can be summed to the seed number k, where the first number 
in the tuple is an integer less than or equal to [k/2]. Here [x] denotes the 
greatest integer less than or equal to x. For example, in Fig. 1, which displays 
the tree diagram for the seed number equal to 6, we draw branch lines to 
(0,6), (1,5), (2,4) and (3,3). Whenever a zero appears in the first element of 
a tuple, the path stops, as evidenced by (0,6). For the other pairs, one draws 
branch lines to all pairs with integers that sum to the second number under 
the prescription that the first member of each new tuple is now less than or 
equal to half its second member. Hence, for (1,5) we get paths branching 
out to (0,5), (1,4) and (2,3), but not (3,2) or (4,1). This recursive approach 
continues until all paths are terminated with a tuple containing a zero as 
indicated in the figure. 




Figure 1: Tree diagram of the partitions summing to 6 for the BRCP algo- 
rithm. 
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It is obvious that all the first members plus the second member of the final 
tuple in each path represents a partition for k. E.g., the path in the figure 
consisting of (1,5), (1,4), (2,2) and (0,2) represents the partition {1,1,2,2}, 
while that consisting of (1,5), (1,4), (1,3) and (0,3) represents the partition 
{1, 1, 1, 3}. In addition, the number of branches along each path represents 
the number of elements in each partition, whilst those tuples with zeros in 
vertical columns represent the partitions with the same number of elements 
in them. For k > 3 the last path in the figure consists of {[k/2]+l,[k/2]) 
and (0,[A;/2]) when k is odd and {k/2,k/2) and (0,^/2) when k is even. In 
both even and odd k cases the tree diagram terminates at what we shall 
call the central partition, which represents the partition of {[A;/2] + 1, [A;/2]} 
for odd values of k and {k/2, k/2} for even values of k. Unfortunately, this 
is not all that is required to produce the final tree diagram in the figure. 
Duplicated paths involving permutations of the same partition must also be 
removed so that each partition appears only once in the final diagram. When 

1 eventually end up with the tree 

k 

m 

with zero in them m branches from the seed number. Hence, 



this removal process is carried out, one wi 
diagram displayed in Fig. [H To determine 



we simply count all the tuples 



equals the 



number of tuples with zero in them in the third column, which comes to 
3. In addition, the number of partitions with exactly m parts is the same 
number as the number of partitions whose largest element is m, which will 
be become apparent when we discuss conjugate partitions. 

In Ref. [T] it was stated initially that duplicated partitions could be re- 
moved by the introduction of a search algorithm. Such an algorithm would 
result in a major increase in the complexity of the above graphical approach 
for generating partitions, which would, in turn, make the other algorithms 
presented earlier not only more appealing, but also far more efficient. To 
avoid the introduction of a search algorithm, it was stated later in the same 
reference that the generation of the partitions could be improved by viewing 
the tree diagrams from a different perspective. Thus, one notes first that the 
entire tree emanating from (1,5) in the diagram is the same tree that one 
would obtain if the graphical method was applied to 5 instead of 6, thereby 
exhibiting the recursive nature of the method. Similarly, the tree emanat- 
ing from (1,4) is the same tree diagram one would obtain by applying the 
graphical method to 4 instead of 6. This continues all the way down to the 
last partition whose elements are only composed of ones. Moreover, only the 
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partitions with unity in them emanate from (1,5) in the diagram, whereas the 
partitions emanating from (2,4) only possess elements greater than or equal 
to 2. Similarly, the partitions emanating from (3,3) only possess elements 
that are greater than or equal to 3. In fact, in the last instance since 3 is 
half of 6, there will only be threes involved along the path from (3,3). Thus, 
we see that the last path represents the central partition {3, 3}, which would 
have been {4, 3} had we constructed a tree diagram for partitions summing 
to 7. 

From the tree diagram we see that the second number in a tuple decre- 
ments with each rightward horizontal movement or right branch, while the 
first element of each tuple increments with downward vertical movement. In 
short, the trees are two-dimensional. That is, two variables are required to 
construct them, a property first observed by D. Balaic. This is particularly 
interesting as it means that we are describing a rare instance of bi-variate 
recursion. Hence, there is no need for the introduction of a search algorithm 
to remove duplicated partitions. Since it has been noted that the tree di- 
agrams terminate at the central partition, we shall refer to the algorithm 
that generates the partitions based on the tree diagrams as the bi-variate 
recursive central partition or BRCP algorithm. 

Before presenting an elementary code utilising the BRCP algorithm, let 
us investigate how the recursive properties of partitions are included in the 
tree diagrams such as Fig. [H From the figure we see that the total number 
of partitions p{k) can be obtained by summing all the partitions of k that 
can be separated into m elements, where m ranges from 1 to k and k = 6 
in this instance. Since there is only one partition with one element and one 
with only k elements, by scanning over the columns in the tree diagram we 
obtain the trivial equation of 



fc-i 



p{k) = 2+J2 



m=2 



(3) 



which is valid for fc > 3. On the other hand, if we scan the rows of the tree 
diagram, then we see that the total number of partitions can also be obtained 
by letting p(fc, m) represent the number of partitions whose elements are 
greater than or equal to m. As a result, we arrive at 

[fc/2] 

p{k) = 1 + \^p{k — rn^m) , (4) 

m=l 
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It has already been stated that 



represents the number of parti- 



where, again, [x] is the greatest integer less than or equal to x. This result is 
given in Ref. [3D] except that the variables in p{k, m) have been interchanged 
and will be used later in Sec. 4 when we introduce the partition operator. 

k + m 
m 

tions summing to k with at most m elements, i.e. P{k,m). If we put m = 2 
and k = 4, then from the diagram that there are three tuples with zeros in 
them in the column two branches away from the seed number, viz. {1,5}, 
{2,4} and {3,3}. Hence, P(4,2) = 3. From Eq. (^ we see that P(4,2) is also 



equal to the sum of 



and 



If we treat the five in the tuple {1,5} in the 



tree diagram as a seed number, then 



is equal to one corresponding to the 



tuple {0,5}. Furthermore, if we now treat the four in the tuple {1,4} from 
{1,5} as a seed number, then we find that two branches further to the right 

4 

equals 2 corresponding to the tuples {0,3} and {0,2}. Note that we could 

not have used the four in the tuple {2,4} in the figure as the seed number 
because this tree diagram gives all the partitions summing to 4, whose ele- 
ments are greater than unity. Including the first branches emanating from 
the seed number of 4, we see that these correspond to the partitions of {1,3} 



and {2,2}. Hence, summing the 



and 



we find once more that P(4, 2) is 



equal to 3, confirming that the tree diagrams do indeed possess the recursive 

properties of partitions. 

k 

also represents the number of partitions 
m 

ement is m. This connection can be observed 

by using Ferrers diagrams, which are studied later in this work. As an exam- 

6 

which can be determined by summing all those tuples 



According to Knuth 
summing to k, whose largest e 



pie, let us consider 

with a zero in the vertical column three branches from the seed number (the 
standard approach) and is, therefore, equal to 3. Meanwhile, the largest el- 
ement of a partition always appears in the final tuple ending a path in the 
tree diagram. Hence, the number of partitions whose largest element is 3 
can be determined by summing all the paths ending with the tuple of (0,3). 
There are three of these in the tree diagram with the first occurring at the 
top of the fourth column from the seed number, the second at the bottom of 
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the third column and the third at the bottom of the second column. 

With regard to the BRCP algorithm, an elementary version in C/C++, 
which first appeared in Ref. p^, is 

void idx(int k,int j) { 

printf("%d" ,k); 

k=k-j ; 

while (k > = j ){ 

printf(",%d(",j); 

idx(k , j++); 

printf(")"); 

} 

} 

Note that the order of the variables in the above code is interchanged com- 
pared with the tuples in the tree diagram. For fc = 4 the output from this 
code is 4, 1(3, 1(2, 1(1))), 2(2). By processing the commas and parentheses, 
we obtain the partitions in the order they appear in the tree diagram, viz. 
{4}, {1,3}, {1,1,2}, {1,1,1,1} and {2,2}. Ahhough the output is very com- 
pact, the code in this form is not suitable for the implementation into the 
second step of the partition method for a power series expansion. In fact, the 
above code has to be adapted in order to solve various problems connected 
with the theory of partitions studied later in Sec. 5. Even if we want to 
list the partitions on separate lines in a similar manner to the other codes, 
modifications are required. Nevertheless, the above code does represent the 
simplest implementation of the BRCP algorithm. It is not only more struc- 
tured and hence, more elegant than the reverse lexicographic algorithm of 
McKay and Knuth, but it is also more powerful or versatile. For example, 
one single call to idx(6, 1) results in all the other calls to the routine as 
shown in the tree diagram. In fact, the total number of calls to idx yields 
the total number of partitions p{k), which is an important quantity in its 
own right. By introducing a counter for the number of calls to idx in the 
above code, we obtain the total number of partitions required to construct a 
tree diagram without the need to create new routine. We shall observe later 
in this work, especially in Sec. 5, how making only a few changes can result 
in a host of the special properties being determined from partitions. 

Let us consider the generation of partitions summing to a particular value 
as we have done for the other codes codes presented in this section. To demon- 
strate the versatility of the BRCP code, we shall not print the partitions as 
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in the same manner as these codes. Instead, we shall generate the partitions 
in the multiplicity representation, which is required for the second step of 
the partition method for a power series expansion. Consequently, a func- 
tion prototype called termgen needs to be introduced into the BRCP code 
called partgen as displayed below. When it is called, it will compute the 
frequencies of the elements by counting the same elements in each partition, 
the latter being represented by the array called part. In order to facilitate 
the call to termgen, idx has undergone minor modification so that in pro- 
cessing the partitions the ones are counted first, the twos next and so on. As 
a result, the following output is obtained for k = 5: 
1 : 1(5) 
2:1(1)1(4) 
3:2(1)1(3) 
4:3(1)1(2) 
5:5(1) 
6:1(1)2(2) 
7:1(2)1(3) 

In the code given below the variable term represents a rolling count of the 
partitions as they are being determined by idx. In the above output the first 
value printed out on each line is the value of term, which is followed by a 
colon. Hence, the final value of term represents the total number of parti- 
tions for each value of k or the variable tot in the code. Each line of output 
only gives the nonzero values of the frequencies of the elements accompanied 
by the values of the elements presented in parentheses. For example, 1(5) 
denotes the partition {5} where n5 = l and all the other n^ equal zero, while 
3(1) 1(2) represents {1,1,1,2}, in which case ni = 3 and n2 = l. As expected, 
the final partition is the central partition for k = 5, viz. {2,3}. 

As a comparison, the code given below was run on the same Sony laptop 
as the other codes, where it was found that it took 1561 CPU seconds to 
compute all the partitions summing to 80 according to the above format. 
Hence, the execution time compared with the other codes has increased, 
primarily due to the extra processing of the partitions. However, if the 
output is directed to a file, which can, in turn, be used as input to the 
partition method for a power series expansion, then it only takes 26 CPU 
seconds to execute, which makes it the best performing code in this mode. 

Whilst the other algorithms/codes discussed in this section may prove to 
be faster than even an optimised version of the BRCP algorithm in other situ- 
ations, they do not possess the versatility or flexibility of the latter. We shall 
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use this versatility when we embark on programming the partition method 
for a power series expansion in Sec. 4. It should also be noted that Refs. 
\T8] and fT6] present extra algorithms for generating partitions according to 
a specific number of parts, whilst the latter reference present another algo- 
rithm which outputs doubly-restricted partitions or where the elements lie 
in a specified range. In Sec. 5 we shall see that only minor modifications to 
the BRCP algorithm are required to solve these problems. That is, there is 
no need to create an entirely different algorithm to solve such problems. 

T^include <stdio.h> 
T^include <memory . h> 
^include <stdlib.h> 

int tot ,* part ; 

long unsigned int term = l; 

void termgen () 

{ 

int freq , i ; 

printf ("%ld: " ,term++); 

for (i=0;i<tot ; i++){ 

freq=part [ i ] ; 

if (freq) printf ("%i (%i ) " ,freq , i+1); 

} 

printf ("\n" ); 

} 



void idx(int p , int q) 

{ 

part [p-l] + + ; 

termgen ( ) ; 

part [p-1] ; 

p -= q; 

while (p >= q) { 

part [q-l] + + ; 

idx(p , q 

part [q++ -1] 
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} 

int niain(int argc,char *argv[]) 

{ 

int i ; 

if(argc ! = 2) printf (" partgen <sum of the partitions >\n" ) ; 

else { 

tot=atoi ( argv [ 1 ] ) ; 

part = (int *) malloc ( tot * sizeof ( int ) ) ; 

if ( part^=NULL) printf (" unable to allocate array\n\n" ) 

else { 

for ( i=0;i<tot ; i++) part[i]=0; 

idx ( tot , 1 ) ; 

free ( part ) ; 

} 

} 
printf ("\n" ); 

return (0) ; 

} 

Another interpretation of the tree diagram in Fig. [1] is to realise that the 
first branch from the seed number is the only single element partition, viz. 
{6}. The partitions summing to n—1 or 5 in this case appear along the second 
branch to {1,5} in the tree diagram. If the partitions summing to n — 1 
have already been stored in an array, then all we need to do is increment 
the number of ones by one in all these partitions to get the partitions for 
n. The next branch emanating from the seed number represents the tree 
diagram for n— 2, but now all the elements in the partitions are greater than 
unity. So, if the partitions summing to n — 2 have been stored previously, 
then we disregard those partitions where the number of ones is non-zero 
and increment the number of twos in the remaining partitions by one to 
get the partitions summing to n. The next branch from the seed number 
represents the tree diagram for n — 3 except that those partitions, where 
the number of ones and twos are non-zero, are now neglected. To get the 
partitions that sum to n, we increment the number of threes by one in the 
remaining partitions. This process continues until [n/2] tree diagrams have 
been processed. Although the new interpretation may lead to less processing 
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of the partitions, it comes at the expense of having to store all previous 
partitions in memory. Nevertheless, we shall return to this interpretation in 
the next section. 

3 The Partition Method for a Power Series 
Expansion 

As mentioned earlier, the second step in the partition method for a power 
series expansion is the more important step. It involves coding partitions 
so that each makes a distinct contribution to the coefficient in the resulting 
power series expansion. As discussed in Refs. P] and [3] these contributions 
need not be numerical in nature as was the case in the computation of the 
reciprocal logarithm numbers in Ref. pj. They can also be functions or 
polynomials, which will become evident later in this work. 

The contribution to a coefficient made by each partition in a tree diagram 
is not only dependent upon the total number of elements or parts in the par- 
titions, N, but also on the elements in the partitions. Originally, when the 
method was devised, the elements were set equal to li, while rii represented 
the number of occasions or frequency each /j appeared in a particular par- 
tition. Therefore, if there were j elements in a partition, then X]i=i ^i ~ ^' 
while X]i=i ^ih = k, which represents the order of the power series expansion. 
Later, it was decided to let tt-i represent the number of ones in a partition, 
n2 represent the number of twos, 77,3, the number of threes and so on. Then 
Y2i=i f^i = ^ ■ The reason why the former approach was adopted initially was 
that it eliminated the redundancy caused by the fact that often, many of the 
Hi were equal to zero in the partitions. For example, ng equals unity in only 
one partition in the tree diagram, while for this partition rii to n^ equal zero. 
On the other hand, the problem with the first approach is that N varies and 
this means that writing down general formulae or expressions involving all 
the partitions is far more awkward. Furthermore, the redundancy due to the 
fact that many of the rij are equal to zero in the partitions has no effect on 
the tree diagram and consequently, on the the BRCP algorithm. Hence, we 
shall adopt the second approach when discussing partitions, which seems to 
be the generally accepted approach used by mathematicians [TS] . 

In Refs. [l]-[3] the partition method was apphed to specific instances 
where the standard method of deriving Taylor/Maclaurin series expansions 
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breaks down. However, Taylor/Maclaurin series expansions are simpler to 
derive when the original function is differentiable. In cases where both the 
partition method can be applied and a Taylor/Maclaurin series can be de- 
rived, they yield identical results even if the resulting expansion is divergent. 
In addition, it was described extensively in both Refs. |2] and [3] how the 
method can be extended to situations where the coefficients may be depen- 
dent upon a variable rather than provide a pure number. Therefore, it should 
be possible to develop a general theorem describing the partition method for a 
power series expansion as is the case with Taylor/Maclaurin series. Of course, 
such a theorem will need to indicate under what conditions the method is 
valid. Moreover, in Refs. [2] and [3] it was described how the method could 
be inverted to yield the power series expansion of the reciprocal function. 
That is, if the method can be applied to a function /(x), then it could often 
be applied to 1/ f{x). A theorem on the method would also need to indicate 
under what conditions it can be inverted. 

Before we present the theorem describing the partition method for a power 
series expansion, we need to introduce some preliminaries. First, in the lit- 
erature a composite function g o /(x) is defined as being equal to g{f{x)). 
Here we shall define a pseudo-composite function ga o f as being equal to 
g{af{x)), where a need not necessarily be a number. Next we need to explain 
the concept of regularisation. For more detailed descriptions of this process 
the reader is referred to Refs. [2]-[3] and [7]-[Tn]. According to this concept, 
when a power series representation for a function is divergent, particularly 
as happens with an asymptotic expansion, it needs to be regularised in order 
to yield meaningful values that are representative of the original function. 
As a consequence, when it is uncertain that a power series representation is 
convergent or when it is known to be divergent, we cannot use the equals 
sign in a mathematical statement. Instead, we introduce the less stringent 
equivalence symbol and refer to the resulting expression as an equivalence 
statement. For example, in Refs. [2] and [7]-[8] it is shown that the geometric 
series, i.e. YlT=o^^ ^^ absolutely convergent for \z\ <1, conditionally conver- 
gent for ^z <1 and \z\ > 1 , undefined for 3fJ z = 1 and divergent for ^z>l. 
In the last two instances it is simply invalid to say that the series is equal 
to anything. However, through the process of regularisation it is found that 
the regularised value of the geometric series is the same value of 1/(1 — z) 
that one obtains as the limit value when the series is either conditionally or 
absolutely convergent. Furthermore, this value is bijective and hence, unique 
for z lying in the principal branch of the complex plane. Therefore, when 
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the series is divergent, we can only say that it is equivalent to its regularised 
value. Because the regularised value is equal to the limit value of the series 
when it is convergent, which is not always the case as discussed in Ch. 4 of 
Ref . [7] , we replace the equals sign by the equivalence symbol for all values of 
z. That is, we can express the geometric series as an equivalence statement 
or equivalence for short by the following statement: 



oo ^ 

E^'-I%7 . v.. (5) 

k=N 

Now we introduce the main theorem in this work. 
Theorem 1. Given that the function fi^z) can be expressed in terms of 
a power series referred to here as the inner power series, in which f{^z^ = 
Y^kLoPkU'' and y = z^, and that the function g{z) can be expressed in terms 
of another or outer power series, i.e. g{z) = h{z) J^'kLo ^^^^ , where h{z) is 
an arbitrary function or number, then for non-zero values of po there exists 
a power series expansion for the quotient of the pseudo-composite functions 
gao f and haO f given by 

J, oo 

In the above equivalence the first few coefficients D^ are given by 

Do = F{apo) , Di = aF«(apo)pi , (7) 

and 

D2 = Y F^^\(^Po)pI + aF^'\apo)p2 . (8) 

A general formula for the coefficients can be derived by analysing the parti- 
tions summing to k. This yields 

fc,[fc/2],[fc/3],...,l k n, 

Dk= E a^^^^H^l^o)!!!! • (9) 

ni,n2,ri,3,...,nfe=0 i=l 

In this equation ^^^^ ni = N, while F^^^apo) represents the iV-th derivative 
of the function F{apo), which, in turn, represents the regularised value of the 
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power series expansion Y^'iLolji^PoY ■ That is, F(apo) = Yl^olJ^^Poy ■ For 
the important case where Po = 0, the coefficients are given by 

fc,[fc/2],[fc/3],...,l k n, 

D,= Y. ^^^'^^'H^ ' (10) 

ni,n2,n3,...,nfc=0 i=l 

Moreover, for Z^qt^O, the inverted quotient of the pseudo-composite functions 
can also be expressed in terms of a power series and is given by 

^-^^^^YEuy' . (11) 

Here the coefficients E^ are found to be 

fc,[fc/2],[fc/3],...,i k m 

Ek= E i-lfNlD.^H^ . (12) 

ni,n2,n3,...,nk=0 i=l 

Finally, the coefficients Dk and Ek satisfy the following recurrence relation: 

fc 
J2DjEk-, = 0. (13) 

j=i 

Remark 1. As in the case of a Taylor/Maclaurin series the power series 
given by Equivalences IQ and fITT]) can be either (1) convergent for all values 
of the variable, (2) absolutely convergent within a finite radius of convergence 
or (3) asymptotic, which is defined here as a power series expansion with zero 
radius of absolute convergence. 

Remark 2. The second result for the D^, viz. Eq. (TTOj) . is similar in form 
to the definition on p. 134 of Ref. [21] for the partial or second type of Bell 
polynomial. In fact, the latter are a special case of the above theorem, which 
can be obtained by setting qk = a = l and pk = Xk/k\. 

Proof. Since g{z) can be expressed in terms of a power series expansion 
and the function h{z), we have 

J^^[lo + Y^.a'fizr] . (14) 
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Introducing the power series expansion for f{z) into the above result yields 

9aOf / ^ " • "^" -2 

k=0 



ha of 



' oo oo „ 

qo + qia^Pk y^ + g2 a^ y^ Pk y^j 



k=0 



oo r, OO . 

+ qsa^ij^Pky'') +ii(^^{^Pky'') + 



(15) 



fc=0 



k=0 



Isolating the zeroth order term of the power series expansion for f(z) in the 
above result yields 

e ( oo oo ^ 

-^^^ - \qo + qi(^\Po + ^Pk y^j + ^2 a^ (po + ^Pky'' 

\ k=l k=l 

cc „ oo . 

fc=l fc=l 

Expanding in descending powers of po yields 



(16) 



|^^5:„(ap„y + 5:j„aVo-'EP'=!'' + 5:(2,*' "•'" 



i=o 



i=i 



k=l 



i=2 



.fc=l 



i=3 



(17) 



vfc=l 



Let us now represent the regularised value of the series '^'^Lo qki.C'Po)'' by 
F(apo)- That is, 



'^qkiapo)'' = F{apo) . 



(18) 



k=0 



We can use this definition to simplify the sums on the rhs of the Equivalence 
flTTl). but first we note that 



E 



J \ 7 J— * 

qjO-Fo 



i-l 



a' d' 



^ £*("?«)'- no -') = ^;^E*-'' 



3=1 - ' 3=\ 1=0 j=0 

Introducing Equivalence flTSj) into the above result yields 



• (19) 



z=apo 



E 



j=i 



^^•"''^°"^ " 7! ^""^^^ 



2! 



F(*)(apo 



(20) 



2=apo 



31 



Consequently, Equivalence (IT7|) can be expressed as 



oo ■ ' ~- ^ ^ 



|^-E^^'"(ap.) EmM . (21) 

■' k=0 \j=l / 

If Equivalence (|2T|) is expanded in powers of y, then we obtain 
1^ = F{apo) + a F(i) (apo)pil/ + ( y i^^'^ (apo) P? + a F^'^ (0^0)^2) 2/' 
+ (^F(3)(apo) + a2F(2)(apo)piP2 + aF«(apo)P3)y' + 0(2/^) . (22) 

From this result we see that the coefficients of the zeroth, first and second 
order terms in y correspond to the results for Dq, Di and D2 given in the 
theorem. Moreover, we see that 

^3 = ^ F^'^ (apo) pI + a'F(2) [apo] P1P2 + aF^'^ (apo) Ps . (23) 

The first few coefficients are relatively easy to write down, but beyond 
that it becomes progressively more difficult. It is at this stage we need to 
introduce partitions into the analysis to facilitate the derivation of a general 
expression for the coefficients of the powers of y in Equivalence (l22l) . 

If we look closely at the result for D3, then we see that it is the sum 
of three separate contributions, which is due to the fact that there are only 
three partitions summing to 3, namely {1,1,1}, {1,2} and {3}. As stated in 
Ref. [l]-[3], in order to evaluate the specific contribution due to each parti- 
tion, a value must be assigned to the elements appearing in the partitions. 
In particular, Ref. [2] states that these assigned values depend upon the co- 
efficients of the inner series, i.e. the power series expansion for f{x), which 
becomes the variable in the outer power series for g{z). Hence, each element 
2 in a partition is assigned a value of p,. Furthermore, there is also a multino- 
mial factor associated with each partition. This factor is not only dependent 
upon the frequencies n, of the elements in a partition, but also on their sum, 
N = Y2-^-^ Hi. It arises from the fact that the value of k on the rhs of Equiv- 
alence fl2T|) is used to render each partition. That is, k in Equivalence fl2Tl) 
corresponds to A^ in the partition method for a power series expansion. Often 
the multinomial factor simply becomes A^!/ni! ^2! ■ ■ ■ n;;.!, but in Equivalence 
fl^T]) there is an extra factor of a^ F'^^^apo) / k\ in the terms. Again, as k in 
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Equivalence (l^Tj) plays the role of A^ in the partition method for a power 
series expansion, this means that the standard multinomial factor must be 
multiplied by a factor of a^ F^^\apQ)/N\ for each partition. Therefore, the 
contribution from a partition to the overall coefficient is given by 

cU,n2,...,n,)=a^FW(apo)n^ • (24) 

We see from this result that we do not actually require the partitions them- 
selves, but their frequencies in order to evaluate the coefficients in the power 
series for the quotient of the pseudo-composite functions. Nevertheless, it 
should be noted that each set of frequencies identifies a distinct partition. 

To make the preceding material more understandable, let us evaluate the 
fourth order term in y in Equivalence (12T1) or D4, which is determined by 
considering all the partitions summing to 4. There are 5 of these: {1,1,1,1}, 
{1,1,2}, {2, 2}, {1,3} and {4}. For the first partition, ni=4, while the other 
Hi are zero. Therefore, from Eq. ( l24l) we have 

C(4,0,0,0) = a^F(^)(apo)pt/4! • (25) 

In the case of the second partition rii = 2, ^2 = 1 and the other nj vanish, 
while for the third partition, only n2{=2) does not vanish. According to Eq. 
021]), the contributions due to both partitions are 



C(2,l,0,0) = a3F(3)(apo)p?P2/3! , (26) 

and 

0(0,2,0,0) = a^F^^\apo)pl/2\ . (27) 

For the fourth partition ni = 713 = 1 with ^2 and n4 equal to zero, while in 
the final partition only 714 (=1) is non-zero. Hence, these partitions yield 

C{1,0,1,0) = a' F^'\apo)piP3/2\ , (28) 

and 

C(0,0,0,l) = aF«(apo)P4 ■ (29) 

Hence, 1^4 is given by the sum of five or p{4) contributions. 
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To derive a general formula for the coefficients D^, "we need to sum over 
all the partitions summing to k. This entails summing over all values that the 
frequencies of the elements can take. Each frequency, rij, can only range from 
to [k/i], since this is the maximum number that the element i can appear 
in a partition summing to k. In addition, the partitions are constrained by 
the condition that 

k 

y riii = 77-1 + 2ri2 + 3ri3 + ■ ■ ■ + kn^ = k . (30) 

4 = 1 

From this result we see that n^ equals zero for all partitions except {k}, 
in which case it will equal unity, while the other rii will equal zero for this 
partition. That is, the rzj are more often zero than non-zero. Consequently, 
much redundancy occurs when summing over the allowed values of rij. More 
succinctly, the sum over all partitions summing to k can be expressed as 

k,[k/2],[k/3],...,l 

Dk= ^ C{ni,n2,...,nk} ■ (31) 

ni,n2,n3,...,nfc=0 

where N = Y2i=i'^i- If Eq- ^SM is introduced into the above equation, then 
we obtain the general formula given by Eq. (II ip . 

Now we consider the case of po = 0. This means that Equivalence (1TB]) 
reduces to 



9a o f 



k 



!kiL^,„ + Y.''-''i'£ny') ■ (32) 

"'^ •^ fc=i i=i 

Expanding the first few powers in y yields 

9a ° f ( 1 iX 1 ( 

T—J = Qo + qi apo y + [qiap2 + q2a Pijy + [qiaps 

+ 2q2a'piP2 + q3aM)y^ + 0{y') . (33) 

Hence, we obtain a power series expansion in y except now the coefficients 
Dk take a different form. In particular, we see that the first few coefficients 
are given by 

Do = qo , Di = qiapo , D2 = qiap2 + q2a^pl , (34) 
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and 

D3 = gi ap3 + 2q2 a^i P2 + Q^a^pl ■ (35) 

Furthermore, we can introduce the partition method to derive a general for- 
mula for coefficients since Equivalences fl2T|) and ( 132|) are isomorphic. In 
fact, the only difference between the two equivalences is that the terms in 
the outer series in the latter equivalence are multiplied by g^a'^ as opposed 
to a^ F^^\apQ) / k\ in Equivalence (121 p . This means that the only change to 
the partition method will occur in the multinomial factor, which will now be 
q^a^ N\ / ni\n2\ ■ ■ -rikl. Consequently, the contribution by each partition to 
the Dk is given by 

c(ni,n2,...,n,)=g^a^iV!]^^ . (36) 



Introducing this result into Eq. (ISTj) gives Eq. Qj. 

Although it has been stated that the expansion given by Equivalence 
([6]) can be asymptotic, this does not mean that it will be divergent for all 
values of y. For if it were, then inverting the equivalence implies that the 
inverted quotient of the pseudo-composite functions vanishes for all values 
of y. Therefore, there must be a region in the complex plane where the 
equivalence symbol can be replaced by an equals sign. In Ref. [2] it was found 
that an asymptotic series possesses zero radius of absolute convergence, but it 
is either conditionally convergent or divergent depending upon which sector 
in the complex plane the variable is situated. For the values of y where 
Equivalence ([6]) is convergent, we can invert the quotient of the pseudo- 
composite functions Qa ° f and ha o /. Therefore, provided Dq 7^ 0, we find 
that 

The rhs of this equation can now be regarded as the regularised value of the 
geometric series where the variable equals to — Yl'k'=i(^k/Do)y''. According 
to Refs. and [I|-[9], the geometric series or 1 + z + z'^ + ■ ■ ■ is either 
conditionally or absolutely convergent for 3fJ 2; < 1. For all other values of z in 
the principal branch of the complex plane, it is either divergent or undefined 
and must be regularised. Hence, treating the rhs of the above equation as 

35 



the limit of the geometric series means that 3? YlT=i(^k/ Do)y'' > —1. For all 
other values of y the series will be either divergent or undefined, the latter 
case occurring when 3fJ '^'^LiiDk/ Do)y'' = —1. Therefore, Eq. f l37|) can be 
expressed as 



hgOf _ 1 

gao f ~ Dq 



k=0 \ j=l 



Equivalence (l38l) is isomorphic to Equivalence (pTl) . which means in turn 
that we can apply the partition method again. In this instance the coefficients 
of y^ in the inner series or rather the values to be assigned to the elements 
j in the partitions are equal to —Dj/Dq instead of pj. In addition, the 
multinomial factor becomes the standard value of A^!/ni!ra2! ■ ■ -Ukl for the 
partitions summing to k. Hence, the contribution from each partition is 
given by 

c(ni,n2,...,n,) = (-l)^iV!Do^n^ • (39) 

1=1 

Introducing the above into Eq. ( l3T|) yields the result given by Eq. ( fT2|l . 
To derive the final result, we use Equivalences (E]) and ( ITTj) to obtain 



oo 



ga° f hgO f _ 



^ oo 

^'^y'i^T.^'^y') ■ (40) 



fc=0 ° fc=0 

Since Eo = l from Eq. fll2p . we can separate the zeroth order term in y, which 
cancels the term of unity on the Ihs of the equivalence. After multiplying the 
series together, we are left with 

oo k 

5^/5^Z},E,_, = . (41) 

fc=i i=o 

As indicated earlier in the proof, there will be some values of y, actually a 
region in the complex plane, where the equivalence symbol can be replaced 
by an equals sign. Otherwise, either Equivalence iQ or fITT]) would be zero 
for all values of y. Given that there will be an infinite number of values of 
y where an equals sign applies, the Ihs of the above result will also vanish 
for these values of y. For this to occur, it means that the inner series, which 
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is independent of y, must also vanish. Hence, we arrive at Eq. (IT^ . thereby 
completing the proof of the theorem. 

It should be noted that Theorem 1 has avoided the issue of determining 
the radius of absolute convergence for the power series on the rhs of Equiva- 
lence ([6]). This is because although one can determine a value for the radius 
of absolute convergence when deriving the resulting power series, it is often 
only an estimate, not the supremum as demonstrated by various examples 
in Ref. [3]. Furthermore, one can introduce a divergent power series as the 
inner series and despite the fact that the outer series may also possess a finite 
radius of absolute convergence, the resulting power series appearing in Equiv- 
alence (E]) niay be convergent over the entire complex plane. For example, 
the power series expansion for cosecant derived via the partition method for a 
power series expansion in Ref. |2] has a radius of absolute convergence equal 
to IT. This inner series became the variable for the outer series, which was 
given by the geometric power series, whose radius of absolute convergence is 
unity. Yet, the resulting power series in Equivalence ([6]) was merely another 
representation of the standard Taylor/Maclaurin series for sine, which is con- 
vergent for all values of the variable. So, this is an example where both the 
inner and outer series were both absolutely convergent within different radii 
in the complex plane, but the resulting series expansion obtained via the 
partition method for a power series expansion was convergent for all values 
of the variable. 

To make the preceding material clearer, let us consider an example. Power 
series expansions for esc z and sec z have already been obtained in Ref. [3] via 
the partition method for a power series expansion. Because cos z and sin z can 
be expressed as Bessel functions of half-integer order, namely J-i/2{z) and 
•^1/2(^)5 we can also say that power series expansions have been developed 
for the reciprocal or inversion of these special functions for u = —1/2 and 
i/ = l/2. Whilst various extensions of these results are presented in Ref. p], 
one which has been overlooked is the derivation of a power series expansion 
for the reciprocal of a Bessel function to arbitrary order z/. Therefore, if we 
introduce the standard power series expansion for Bessel functions given by 
No. 8.440 in Ref. ESI into the denominator, then we obtain 



"^ ^ (l-zV4(z/+l) + 2Vl6-2!(z/ + 2)(z/ + l) )' ^ ' 

The above result is not valid for 3fiz/ = — 1 since in this case the leading term 
in the series expansion for Ju{z) vanishes. Then we would need to examine 

37 



this case by itself, which is left for the reader to consider. From Theorem 1 
we see that the the inner power series is simply the Taylor/Maclaurin power 
series expansion for Bessel functions. Hence, y = z'^, and po = 0, while for 
k > 1, Pk = (— l)*^/2^*^(i^ + l)A:fc!, where (i/ + l)fc is the Pochhammer notation 
for r{k + z/ + l)/r(z/ + 1). Moreover, we can regard the denominator on the 
rhs of Eq. (j32]) as the regularised value of the geometric series, which means 
that coefficients of the outer series, viz. q^, are equal to (—1)'^ for k > 1, 
while the function h{z) = {2/zYV{y + 1). Since a = l, the pseudo-composite 
functions become composite functions with the quotient in Theorem 1 equal 
to {z/2Y /V{u + 1) Jy{z). By using Equivalence (j6]) we arrive at 



n>' + i)JAA fci ^z-" 



fc=0 

where the coefficients hk^i^) are determined from Eq. (ITOl) and are given by 

fc,[fc/2],[fc/3],...,l k 



k,[k/2\,[k/3\,...,l k 



ni,n2,n3,...,nfe=0 i=l 

From this result we see that holu) = 1, hi{u) = l/{v + 1), and h2{i') 
(z/ + 3)/2(z/ + l)2(zy + 2), while for A; = 3, Eq. ^ yields 

M.) = I ( £i^V- 2! £1^ ^<^;±iL + r(.+i) 



r(z/ + 2)/ r(z/ + 2) 2! -1(1^ + 3) 3!-r(z/ + 4) 

z/2 + 8i/ + 19 
3!-(z/ + l)3(z/ + 2)(z/ + 3) 



(45) 



Since ^1/2(2) is related to sinz, the power series expansion for 1/Jij2{z) 
is identical to the one derived in Ref. [3] for cosecant, whose coefficients were 
expressed in terms of the cosecant numbers denoted by c^. These numbers 
were later found to be related to the Riemann zeta function. Hence, we 
obtain 

hk{ll2)=2''ck = 2{2''~2)^-^ . (46) 

Similarly, because J-i/2{z) is related to cos 2;, the expansion for secant, which 
is derived in terms of special numbers called the secant numbers or dk in Ref. 
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[5], is related to z/ = —1/2 in Equivalence (142 p . It was then found that 
the secant numbers could be expressed as the difference of specific values of 
Hurwitz zeta function. Thus, we arrive at 

h,{-l/2)=2^^d, = -i^ (C {2k + 1, 1/4) -a2k+ 1, 3/4)) . (47) 

It was also observed in Ref . [3] that the absolute convergence of the power 
series series expansions for both cosecant and secant were determined by 
the distance from the origin to the first zeros of these functions. In the 
case of cosecant the expansion is absolutely convergent for \z\ < tt, while 
in the case of secant it is absolutely convergent for 1^1 < 7r/2. Therefore, 
the expansion given in Equivalence fH3l) will be absolutely convergent for \z\ 
less than the magnitude of the first zero for Ju{z). For these values of z 
we can replace the equivalence symbol by an equals sign, thereby producing 
an equation. Consequently, we can use the "equation form" of Equivalence 
(143|) to demonstrate that the coefficients hki^i') can also be evaluated by 
recurrence relations. Because there is an infinite number of values of z where 
the "equation form" is valid, a recurrence relation can be be obtained simply 
by multiplying Equivalence (H5]l by the cnovergent power series expansion for 
Ju{z) and setting the resulting product equivalent to unity. Since z is still 
fairly arbitrary within the radius of absolute convergence, we can equate like 
powers of z of the resulting equation. This yields 

"'('^) - Zj (fc-j)!r(*-, + . + i) ■ <***' 

Other recurrence relations can be developed by using integral results in- 
volving Bessel functions. For example, one can express No. 3.768(9) in Ref. 
[23] as 

f\xx^-\l-xr-'-^^^^ = V^ay^~^Ti^^)cos(^) . (49) 

We have already seen that there exists a finite radius of absolute convergence 
for z, where Equivalence ( BTj) becomes an equation. This means that there 
is a region in the complex plane where the expansion for the reciprocal of 
the Bessel function can be introduced into the above result without the ne- 
cessity to replace the equals sign by an equivalence symbol. Furthermore, by 
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introducing the power series expansion for cosine into Eq. (1491) we arrive at 
io ^ ^ f:^ £;^ 24^ (2/c-2j)! '^^ ' ' 



k=0 j=0 



v/ir(/x) f. (-1)'= /ax2. 



22M-ir(/i + l/2)^ {2k)\ \2) 

The integral on the Ihs of the above equivalence is merely the integral repre- 
sentation for the beta function. Hence, we can introduce the gamma function 
product for the beta function. Since a is fairly arbitrary, we can equate like 
powers of a on both sides of the equation, thereby obtaining 



.^P 24i(2A;-2j)! T{2k -2] + 2v + 1) ^^ ' 

^ (-1)^ 1 



22k+2v (2A;)! r(z/+l) 



(51) 



where fi — 1/2 has been replaced by u. 

If we put k = 4 into either recurrence relation and introduce the values of 
/ifc(z/) for /c = l to 3 given earlier, then we find that 

, . _ z/^ + 17z/3 + 117i/^ + 379t/ + 422 
'^''^"4!(z. + l)M^ + 2)M^ + 3)(z/ + 4) • ^^^^ 

This agrees with the value given in Table [H which displays the values of 
/ifc(i/) up to A; = 7. These values have been evaluated by introducing the 
recurrence relation into the Sum routine in Mathematica [12]. In the next 
section we shall develop a programming methodology where summing the 
contributions from the partitions in Eq. (jH]) will be just as expedient as 
using the recurrence relations. 

From the table it can be seen that the /ifc(z/) possess common properties. 
For example, there is always a factor of fc!(l + u)'' in the denominator. In 
fact, if the denominator of a function f{x) is denoted by DN{f{x)), then we 
have is given by 

DN{hkiiy)) =k\{iy + 1)'' {v + 2f'l''^ ■ ■ ■ {u + k ~ if/C'-')^ (u + k) , (53) 

where [x] denotes the greatest integer less than or equal to x. This result also 
follows from Eq. 044p when we examine the upper limits in the summations. 
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1 

1 

u+3 



3\iu+iy{u+2){u+3) 
u'^+17u^ + n7iy'^+379u+422 
4!-(iy+l)4(i/+2)2{,y+3)(i^+4) 

5!-(!/+l)5(i/+2)2(jy+3)(i/+4)(i/+5) 

iy^+42i/^+811i^'^+9412i/^+71155i^"^+349786i/^+1043637i^^+1674616i/+1091052 

6!-(i'+l)6(!^+2)3(iy+3)2(;/+4)(!y+5)(i^+6) 

^^+55i^^+1417i^'^+22535i^^+243311i^-'"'+1827401i/-*+9292435i^^+29539597i/2+51572980i^+36978156 

7!-(i/+l)7(i/+2)3(i.+3)2(z/+4)(i^+5)(i^+6)(i/+7) 



Table 1: Coefficients for the power series expansion of the reciprocal of the 
Bessel function of order u. 

On the other hand, the highest order term in u in the numerator is always 
k less than the highest order term in the denominator. Therefore, for large 
\iy\, hk{iy) ^v-^/k\. 

By applying d'Alembert's ratio test as described on p. 24 of Ref. [21] to 
the series on the rhs of Equivalence fH3|) . we find that for all k the series 
is only absolutely convergent when \z\ < 2^J\hk{l')/hk+l{l')\. Therefore, if 
there is a supremum for this ratio that applies to all values of k, then it 
represents the radius of absolute convergence. Fig. |2] presents a graph of 
hkiy) / hk+iiy) versus z/ for — 2 < z/ < 10 and various values of k. The figure 
shows that the larger k is, the greater the value of the ratio hkii') / hk^i{i>) . 
Therefore, the limit as A; — )■ cxd represents represents the supremum for the 
ratio when z/>— 2 with the nearest singularities to the origin occurring at 



.o,. = ±2 1im J-M^ . (54) 

The figure also shows that as k increases, hi^{v) / hk^i{v) becomes a bet- 
ter approximation to \mik^oohk{i')/hk+i{i') over an ever-increasing range of 
values for v. Eg., we see that hr,[i/) / helu)) is an accurate approximation 
for —2 < u < 4.5, while hglu) / hglu) is an even better approximation for 
-2<z/<7. 

According to p. 372 of Ref. [25], for z/ > — 1, the ffist zero for Ju{z) is 
real, but for u < —1, excluding negative integers, or u complex, the ffist 
zero is complex. Since the ffist zero represents the singularity of 1/Ju{z) 
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Figure 2: hk{u) / hk^i{u) versus v ioi k = 1, 5, 8, 12 and 18 
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V 


BesselJZero 


±2^h^,{y)lh^^{y) 





2.404 825 557695 7 


±2.404 825 557695 4 


1 


3.831705 97 


±3.831705 96 


2 


5.135 622 3 


±5.135 6221 


5 


8.7714 


8.7713 


10 


14.47 


14.46 


-1/3 


1.866 350 858 873 4 


1.866 350 858 873 8 


-4/5 


0.936 806 664 511 


0.936 806 664 510 


-3/2 


2.798 386 045 783 878 


±1.199 678 640 257655 2 


-1/3 + i 


Unable to evaluate 


±2.076 341434 394476 
±1.556 637 759 994 043 z 


3/2 -z 


Unable to evaluate 


±4.529 756 943 967303 
Tl.293 935107111323z 



Table 2: Evaluation of the first zero of J,y(a;) using BesselJZero in Mathe- 
matica and Eq. (15^ with A; =17. 



nearest to the origin, we expect the rhs to be real only for v > —1. Table [2] 
presents a sample of the values for the first zero of Ju{x) for various values 
of z/ obtained by using the special numerical routine known as BesselJZero 
in Mathematica and by setting fc = 17 in Eq. ( 15^ . For z/ > — 1 we see that 
the values obtained via both approaches agree. In particular, for — 1 < z/< 1 
the values obtained from Eq. ( 154|) were found to be accurate to at least 10 
decimal places. As expected, as u becomes greater than unity, the accuracy 
of Eq. (15^ begins to wane compared with the values from BesselJZero. The 
values presented in the second and third columns for these values of z/ only 
agree to the first decimal place. Mathematica is able to evaluate the zeros 
to arbitrary precision, as evidenced by the values for u equal to -1 and - 
3/2 in the second column of the table. On the other hand, if one wishes to 
obtain more accurate values of the zeros via Eq. ( 15^ . then one will have to 
determine the hk{i') for much greater values of k. 

For z/< — 1, BesselJZero [z/,1] gives the first zero of Ju{z), but it is the first 
zero situated on the positive real axis, not necessarily the closest to the origin. 
To see this more clearly, for z/ = —3/2 Mathematica gives as its first zero 
z = 2.798 386 ■ ■ ■ , but this is not the first zero in the complex plane. According 
to Eq. ( 15^ . the closest zero is given by 2 = 1.199 687 ■ ■ -i as displayed in the 
adjacent column of the table. The latter value can be confirmed as a zero by 
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introducing it into the BesselJ routine in Mathematica, whereupon one finds 
that it yields the tiny value of (7.35 ■ • • + 7.35 ■■■i)x 10"^^. 

It also appears that for 5Rz/ > —2, Eq. (1541) still gives the nearest zero 
to the origin of the complex plane. For example, the remaining values in 
the table represent the zeros obtained for complex orders of Bessel functions. 
Zeros for such functions cannot be obtained via BesselJZero since the routine 
can only handle real numbers. However, Eq. f lM|) can yield them. When the 
value obtained from Eq. (15^ for J_i/3+j(2;) is introduced into the BesselJ 
routine in Mathematica, a complex value with a magnitude of the order of 
10^^^ is obtained, while a complex value with a magnitude of the order of 
10^^ is obtained when the final value in the third column of the table is 
introduced into the routine. This last example demonstrates again that the 
accuracy of Eq. flMl) wanes as |z/| increases when k is fixed, which we have 
already observed in Fig. |2j 

The reader may well ask if it is possible to adapt the preceding method 
to evaluate the next Bessel zero or even higher order zeros. This does not 
seem to be possible at this stage unless the specific value for zq^, can be 
determined. For example, by taking the logarithm of No. 8.544 in Ref. [23] 
and differentiating, we eventually arrive at 

In the above equation, J^(z) can be written as a series in powers of zjl with 
coefficients a^ = 1/fc! (i/ + 1)^, while the rhs of Equivalence fl5^ can be used to 
replace the reciprocal of the Bessel function. Multiplying both series yields 

S|-^ii:(r-<^). (-) 

where efc(z/) = Ylij=(d'^k-jhk{i')- To carry out the singularity analysis for 
the next zero, we need to remove the A; = term from the rhs of Eq. (I55p . 
Hence, the power series for this term would need to be introduced into the Ihs 
with the above power series, yielding another power series with coefficients 
expressed in terms of ek{v) and powers of zo,i/- If an approximation for the 
latter is introduced, then the coefficients would only be approximate. For 
large values of k they would still be affected by Zq^^ and thus, it would not 
be possible to isolate the next zero. On the other hand, if the exact result 
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for the ^0,1^5 which is known for z/ = ihl/2, then the series would represent the 
power series for the rhs of Eq. (1551) with the sum beginning at A; = 1 or with 
the singularity at the second zero. This situation is discussed immediately 
below Eq. (83) in Ref. [3]. It should be noted that there are typographical 
errors there since the power of l/(27r) in the expression for c^ should be k, 
not k + 1 and the limit in the next line below should refer to c^_(_i/c^ , not 

Cfc+l/Cfc. 

Before investigating how the BRCP algorithm can be implemented to 
evaluate the coefficients Dk and Ek in Theorem 1, we now prove an interesting 
corollary to the theorem, but before we can do this, the following lemma is 
required: 

Lemma 1. For a complex, regularisation of the binomial series yields 






Remark. This lemma represents the generalisation of the regularisation 
of the geometric series, which reduces to the latter for a = l. 

Proof. The proof of Lemma 1 appears immediately below Proposition 1 
in Ch. 4 of Ref. [7]. Although the proof is concerned with real values of a, 
there is no reason why a cannot be complex. This is due to the fact that the 
proof involves differentiating and integrating the integral representation for 
the beta function, viz. 

/I ifc+a— 1 

This integral is not only defined for real values of a, but also for complex 
values of a. Its convergence is limited by the values for the real part of 
a, not its imaginary part. Since differentiation and integration only affect 
the real part of a, the proof of Proposition 1 in Ref. [7] can be extended to 
complex values of a. Moreover, the result in Lemma 1 can be simplified by 
replacing the equals sign with an equivalence symbol. Then we have the one 
equivalence statement, which is valid for all values of z and a. 

One final remark is in order. The reader should note that a similar no- 
tation to the standard notation for generalised hypergeometric functions has 
been introduced by expressing the binomial series as iJ-q. The reason for 
the slight variation in notation is that the generalised hypergeometric func- 
tion notation of pFq{ai, . . . , a^; /3i, . . . , /3g; z) is only valid when the series is 
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absolutely convergent, i.e. when p < q oi for \z\ < 1 when p = q [2S]- The 
conditions in the lemma are more general. That is, the binomial series has 
been taken into the regions of the complex plane where it becomes either con- 
ditionally convergent or divergent. Hence, it has been necessary to alter the 
notation to avoid any confusion with the standard notation. This completes 
the proof of the lemma. 

We are now in a position to generalise the derivation of the coefficients 
in Theorem 1 by considering the following corollary. 

Corollary 1 to Theorem 1. Given the same conditions on the pseudo- 
composite functions ga° f and haO f as in Theorem 1, there exists a power 
series expansion for the quotient of the pseudo-composite functions raised to 
an arbitrary power p which is given by 

For k>l, the generalised coefficients for the Dk of Theorem 1 or Dk{p) in 
the above result are given by 

k,[k/2Uk/3],...,l k ip.-.ni 

Dk{p)= Yl (-/^)^^o-^n^l^ ' (60) 






and 



k,[k/2],[k/3],...,i ^ (_n ( W- 

Dk{-p) = j: N\ Do(p)-^-^ n ni ' (^'^ 

ni,n2,n3,...,nfc=0 i=l *' 

where {p)n denotes the Pochhammer notation for r{N+p)/r{p). In addition, 
for p = p + v the coefficients satisfy the following recurrence relation: 

k 
D,{p) = J2d,{p)D,.,{u) . (62) 

j=0 

Remark 1. By putting p = — 1 in Eq. fl60l) we see immediately that the 
Di:{—1) become the E^ in Theorem 1 given by Eq. f fT2|) . Furthermore, the 
Dfc(l) represent another representation for the D^ given by either Eq. ([9]) or 
f ITU]) . the latter being valid when Po = 0. 
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Remark 2. Note that when p equals a non-negative integer j, Eq. ( 160 p 
simphfies dramatically because {—p)n vanishes for N>j, reflecting the fact 
that we are dealing with a finite polynomial of order j. For all other values 
of J the coefficients Dk{p) represent polynomials in p of order k. 

Proof. Since the pseudo-composite functions ga° f and haO f are subject 
to the conditions in Theorem 1, we know that there exists a power series 
expansion where 

" ■' k=0 

We also know that there will be a region in the complex plane where the 
equivalence symbol can be replaced by an equals sign. Separating the zeroth 
order term in above result yields 

J, oo 

" ^ k=l 

We can treat the series on the rhs as the variable in the regularised value of 
the binomial series. Then according to Lemma 1 we have 

"■ ■' k=0 ^ ^' j=l 

Now the above equivalence is isomorphic to Equivalence fl2Tl) . which means 
that it is in the form where the partition method can be applied. As stated 
in the proof to Theorem 1, the values to be assigned to elements i in the 
partitions are given by coefficients of y* in the inner series. Therefore, in this 
case each element i is assigned a value of —Di/DQ. To evaluate the contribu- 
tion from each partition, we need to multiply the product of all the assigned 
values by a factor consisting of the multinomial factor N\/ni\ n-i'. ■ ■ ■ n^! and 
r(A^ — p)/A^! r(— p). The latter term arises from the fact that k in the coeffi- 
cient of the outer series plays the role of A^ in the partition method. There- 
fore, by introducing the Pochhammer notation of (p)^ for r(iV-|-p)/r(p), we 
find that the contribution due to each partition is 

C(ni,n2,...,n,)=(-lf (-p^Do^^n^ , (66) 
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As expected, for p = — 1 the above result reduces to Eq. (15^ . Furthermore, 
by summing over all partitions summing to k, we obtain the total coefficient 
in the resulting power series expansion, which is given by Eq. fl60|) . 
As a result of establishing Equivalence f l59|) . we have 

"■ ■' k=0 

We also know that there will be a region in the complex plane where the 
equivalence symbol can be replaced by an equals sign in Equivalence fl59|) . 
For these values of y we can invert the equivalence, thereby obtaining 

'haof\p 11 



Theorem 1 can again be applied to this result by treating the rhs as the 
regularised value of the geometric series with the variable equal to the series 
in the denominator. In this case h{z) = 1, q^ = (—1)*^, a = 1 and pk = 
Dk[p)/DQ{p) for A;>0, while po = 0. The coefficients of the resulting power 
series expansion can be determined by Eq. f lTO|) . Moreover, the resulting 
power series expansion will be equal to the power series expansion on the rhs 
of Equivalence f l67|) since both have the same regularised value. Hence, we 
can equate like powers of both expansions, which yields 

D,{-p) = D,{pr E iV!(-l)"n 1? • (69) 

ni,n2,n3,.--:"fc=0 i=l * 

By taking the factor of (—1)^ inside the product and 1/Dq{p) outside of it, 
we obtain Eq. f l6T|) . 

Since p = /i+z/, we have 

' 9a° f y _ f 9a° f Yf 9a° f y /^qx 

-haOf) \haOf) \haOf) 

We also note that there will be a region of complex plane for y in Equiva- 
lence dSH]) where we can replace the equivalence symbol by an equals sign. 
Therefore, for these values of y, Eq. (!70|) yields 

oo oo oo oo fc 

E^^(/^)^' = E^^(^)^'E^'^(^)2/' = E^'E^^(^)^'^-^('") • c^i) 

fc=0 A;=0 fc=0 fc=0 3=0 
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Since there is an infinite number values of y for which the above equation 
holds, we can equate like powers of y. As a consequence, we obtain Eq. (1^ . 
which completes the proof of the corollary. 

To make the preceding material more concrete, let us consider a couple of 
examples. In order to simplify these examples, we now regard the sum over 
partitions as a discrete operator, which will be denoted by -^p[-]. We shall 
refer to this form for the sum over partitions summing to k as the partition 
operator. That is, the partition operator is defined as 

fc,[fc/2],[fc/3],...,l 

Lp,k[]= Yl ■ (72) 

ni,n2,n3,...,rH:=Q 

For the situation where the partition operator acts on unity, it yields p{k) or 
the number of partitions summing to k. Hence, Lp^A:[l] = p{k)- The number 
of partitions summing to k can also be obtained from the following generating 
function 118i: 



nz) = X{^^Y.p^k)z' ■ (73) 

m=l A:=0 

Note the appearance of the equivalence symbol in the above result since the 
Ihs can become divergent. This is because the Ihs has been treated as an 
infinite product of regularised values of the geometric series in obtaining the 
power series on the rhs. This will become clearer in Sec. 6, where we shall 
also extend the above result. 

By applying Theorem 1 to the simple case where po = 0, pk = b^, h{z) = 1, 
a = l, gfc = l and y = z, we find that the quotient of the composite functions 
becomes 



2bz 2 

fc=0 



h7 1 °° 

9{f{z)) = 1 + ^3^ = 1 + - J2^2bzr' . (74) 

Hence, via Eq. ( TTOl) we arrive at 



k 
^P,k 



^fc 



where, as before, N^ = ^j^i^^i- The k subscript has been introduced here 
for the first time as it will become apparent that we shall need to sum the 
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n.i to different limits sfiortly. If we clioose Pk = {~^^Y ^"^^ 1k = ("1)*^ instead, 
tlien following the same procedure we obtain 



'P,k 



\Nk 



''^'m. 



(76) 



«=1 



where k >2. The above is an interesting result where the partitions with 
an even number of elements are cancelled by those with an odd number of 
elements according to the frequencies of the elements. 

Eq. fl76p is not the only instance where the sum over all partitions van- 
ishes. For example, consider the application of Theorem 1 to the function 
f{z) = exp(a ln(l + z)) or f{z) = (1 + z)"". Here, the coefficients of the inner 
series are given by po = and pk = {—l)'^~^^/k for A; > 1, while the coefficients 
of the outer series are given by qk = l/k\. Then from Theorem 1 we obtain 



D, 



i-lfLp,t 



_l)N,^N, 



n 



«"■» riA 



(77) 



In Ref. [7] it is shown that f{z) represents the regularised value of the bino- 
mial series. That is, for all values of a and z, we have 



ES^(-)'-(i+-r 



fc=0 



rr-a)A;! 



(78) 



For l-z] < 1, the series is absolutely convergent and we can replace the equiva- 
lence symbol by an equals sign. Since the Dk are the coefficients of the power 
series expansion in z, they are equal to the coefficients in the above result. 
Then we find that 



^P,k 



-D^^a^^ 



n 



Hi 



T{k 



T(-a)k\ 



(79) 



The results in Theorem 1 are actually more general than the above result. 
As a consequence, for po = we arrive at 



I?fr 



^P,k 



g^a^^iVfclD- 



N, 



n 



riA 



(80) 
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and 



Ek — Lpk 



■i)-'iv.m-*n?S 



j=l 



l"-^ Tljl 



(81) 



When a = 1, we have f{z) = 1 + z, which, in turn, means that D^ = for 
k>l. Then we find that 



Lph 



-if^ n 



j=i 



2"-* rii ! J 







(82) 



for k>l. On the other hand, when a = —1, we have f{z) = 1/(1 + z), which 
represents the regularised value for the geometric series. Since the coefficients 
of the latter series are equal to (— l)'^, Eq. fITTj) reduces to 



^P,k 



n 



1=1 



1 



z^'-^ m ! J 



1 . 



(83) 



More importantly, the above results can first be generalised by letting a = l, 
where / is an arbitrary integer. Then f{z) = (1+^)', whose coefficients courtesy 
of the binomial theorem equal (^) for fc < / and vanish for the remaining values 
of /. As a result, Eq. ( TfTl) yields 



^P,k 



-If' n - 



1 = 1 



z"' riil 







;-i)iD 



k > I , 
k<l . 



(84) 



If —a is replaced by a in Eq. fl79|) . then the equation reduces to 

k 



k\L 



P,k 



a 



Nk 



n 

i=l 



«"» n,-! 



a fc 



(85) 



According to Chs. 24 and 18 of Refs. [25] and [27] respectively, the Pochham- 
mer polynomials can be written as 



(«). = (-i)'E(-i)''5- 



U) 



a-' 



(86) 



i=o 



:{i) 



where S"^ are known as the Stirling numbers of the first kind and satisfy 



Aj) 



:{i-i) 






U) 



(87) 
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We can proceed further with Eq. (153]) by introducing a new operator that 
only considers a fixed number of elements in the partitions. By setting this 
number equal to j , we can define the operator for a fixed number of elements 
as 



fc,[fc/2],[A;/3],...,l 
ni,n2,n3,...,ni;=0 

Moreover, the above operator is related to the partition operator by 

k 



(88) 






P,k[ 



(89) 



As mentioned in Sec. 2, the number of partitions of k with exactly j parts is 

k 

This means that 



denoted by 



LU^] 



while the recurrence relation given by Eq. ([T]) becomes 



(90) 



(91) 



Furthermore, introducing the new operator into Eq. fl85|) with the rhs re- 
placed with the aid of Eq. fl86|) yields 



^p,k 



n 



i=l 



«"» n,-!J 



■l)^'^' .0) 



k\ 



5 



(92) 



From the recurrence relation given by Eq. ( 1871) . we obtain 

_fc+l ^ k ^ k 

i"-^ rij ! 



fc+i , 

i=l 



L'p~k\T\-^]+kL'pk\f]-^ 

4 = 1 * i = l 



(93) 



Specific results for the Stirling numbers of the first kind when j is rel- 
atively small have been derived in Refs. [1] and [21]. For example, these 
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references give SJ, —{—l)''T{k)Hi{k), where Hi[k) = X]j=i l/j- Therefore, 
we find that 



'P,k 



k ^ 

[11 irii^A 



k 



Hi{k) . 



(94) 



Values for the Stirhng numbers of the first kind are presented for j close to 
k in the appendix of I 
This, in turn, leads to 



k in the appendix of Ref. [2]. For example, when j = k — l, Sj^ = ~{'^- 



k 



«"» n,-! 



k\ 



l)^ fk 
2 



(95) 



Another fundamental result can be obtained by applying Theorem 1 to 
exp(— x). By writing the function as l/exp(a;), we see that the coefficients of 
the inner series, viz. p^, are equal to l/k\ for A;> 1, while Po = 0. Meanwhile, 
the outer series is given by the geometric series so that qu = (—1)'^- The 
coefficients of the power series for exp(— x) are (— l)^//c!, which are also 
equal to the D^. Therefore, according to Theorem 1 we have 



Tk 

^p,k 



mn 



i=l 



?! ' » Ujl 



k\ 



(96) 



In the above result the i = k term in the product is simply the result on the 
rhs. Therefore, it can be simplified to 



fc-i 



'P,k 



mu 



1 



\ni 



i=l 



rii 



-ly 



k\ 



(97) 



Note in the above result that even though n^ = 0, the constraint in the 
partition operator still applies to k. 

In Ref. [3] we found that the cosecant numbers denoted by Ck were the 
coefficients generated when the partition method for a power series expansion 
was applied to scscs. This means that the method was basically applied to 



scscs 



1 - s2/3! + sV5! - sV7! + 



(98) 
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By applying Theorem 1 to this example, in which h{z) = z, we see that y = s'^, 
a = l, pk = {—l)^~^^/{2k + 1)! with Po = and g^ = 1 since the outer series 
corresponds to the geometric series. Therefore, from Theorem 1 we obtain 



scsc s 



= Y,CkS^'' , 



(99) 



A:=0 



where according to Eq. (TIU]) . 



-Ifck 



Lph 



N,\\{ 



i=l 



{2i + l)\J n,\ 



(100) 



and Nk = Eti^i- It should be noted that nti(-l)*"'"^"' = (-1)''"^^', 
although we shall retain the phase factor of (— !)"■• in order to observe a 
remarkable correspondence arising from the inversion of Equivalence (|99|) . 
It was also found in Ref. [3] that the radius of absolute convergence for 
the power series expansion in Equivalence (l99l) was vr, while the cosecant 
numbers were seen to be rapidly decreasing positive fractions given by Eq. 
( H6|l . Hence, Eq. (IIOOI) represents a means of determining even integer values 
of the Riemann zeta function. 

Now we invert the power series expansion in Equivalence (p9|) and apply 
Theorem 1 again, but in the case we have pk = — c^ with pq = 0. The resulting 
power series expansion was found to yield the standard Taylor/Maclaurin 
power series for sin(s)/s or YlT=oi~^)'' ^'^'^ / C^^ + ^)'' which is convergent for 
all values of s, despite the fact that Equivalence f l99|) has a radius of absolute 
convergence equal to vr. This confirms the earlier remark concerning the fact 
that the resulting power series expansion obtained via Theorem 1 can turn 
out to be convergent even though the inner or the outer series may in fact 
be divergent. Furthermore, from Eq. ( ITOj) we obtain 



:-ii 



{2k + l)\ 



^P,k 



NklU^-c 



i=l 



1 

nA 



(101) 



It should also be mentioned that in Ref. [3] numerous recurrence relations 
were derived for the cosecant numbers. One of these is 



fc-i 

E 

i=o 



(_l)fc-i-i 
{2k-2j + iy. ^^ 



Ck 



(102) 
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where Cq = 1. If we introduce Eq. fllOOp into the Eq. (11021) . then we obtain 
the interesting result of 



fc-i 

E 

i=o 



i2k-2j + l)\ 



L 



PJ 



^P,k 



-i)^^iv,!n 



j=i 



1 



(2z + l)! 



rii -| 

nA 



(22 + 1)! 



Hi 



(103) 



Similar results to the above can be obtained by considering the other recur- 
rence relations. Furthermore, in the same reference it was found that just 
as the Bernoulli numbers give rise to Bernoulli polynomials, the cosecant 
numbers give rise to their own polynomials, which are related to the former. 
In particular, the value at unity was found to be given by 



Ckil) = i-ir2''Lp, 



2k 



2k 



-i)^-iV2,!n 



i=l 



1 



1 ! 



rii 



1 

nA 



(104) 



^2A: 



where A^2fc = X]j=i ^i i^ ^^e above result. Consequently, we see the reason 
for the introduction of the fc-subscript to N in Theorem 1. The value of the 
cosecant polynomials at unity was found to equal 



Cfc(l) 



Ck 



)l-2k 



(105) 



Hence, by introducing Eq. fllOOp into the above result and then equating it 
to Eq. (I104p . we arrive at 



^P,2k 



2k 



1 



i=l 



I 



(22fc 



L 



P,k 



(^ + 1)!^ ^ 
1 



{-if^Nklll 



i=l 



i2i + l)\J n. 



(106) 



Consequently, we do not need to consider all the partitions up to 2k to de- 
termine the sum on the Ihs, which is a significant reduction in computational 
effort. 

In Ref. p] another infinite set of related numbers denoted by dk and 
known as the secant numbers were obtained when the partition method for 
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a power series expansion was applied to sec s. Specifically, the method was 
applied to 



sees 



1 - sV2! + sV4! - sV6! + ■ 



(107) 



By applying Theorem 1 to the above result we have y = s"^, po = 0, a = 1 
and qk = l as before, but on this occasion, pk = {—l)''~^^/{2k)\. The resulting 
power series expansion, which can be expressed as 



J] 4 



sec s = > dkS 

k=0 



2k 



(108) 



where the coefficients dk from Eq. ( ITOj) are given by 



-D'dk 



L 



p,k 



NklU 



i=l 



(2z) 



rii 



(109) 



and Nk is the same sum over the frequencies as before. Equivalence O107p 
was found to possess a narrower radius of absolute convergence compared 
with the power series expansion for cosecant, viz. 7r/2 as opposed to tt, while 
the dk or secant numbers were found to be not as rapidly decreasing fractions 
as their cosecant counterparts. In addition, instead of being related to the 
Riemann zeta function as the cosecant numbers are, they were found to be 
related to the Hurwitz zeta function by 



)2fc+2 



(110) 



\2fc+l 



^'^ = ^^^fcTT (C(2^ + 1' 1/4) - C(2A; + 1, 3/4) 

The bracketed expression can also be written as X]fci(~l)"'^^/(2j 

To invert the analysis, either we can apply Theorem 1 to the power series 

expansion in Equivalence (llOSp or we can go directly to Equivalence (II ip . In 

the latter case we replace Di by rfj in Eq. flT^ . while the Ihs of Equivalence 

( ITTj) equals coss, which we replace by its power series expansion. Then we 

find that the coefficients E^ in Eq. ( !T2|) equal (—1)^^/(2^)!. Hence, we arrive 

at 



-11 



(2fc)! 



^P,k 



Nkifii-d^r^ 



;iii) 
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The secant numbers were also found to obey recurrence relations, al- 
though not as many as their cosecant counterparts. Nevertheless, an ana- 
logue of Eq. fll02p was obtained, which is given by 



k-l 



^ (2* - 2i)! "' 



dk 5 



(112) 



with do = l. Introducing Eq. fillip into the above result yields 



V 1 r 

^j=0 {2k-2j)\ ^PJ 



-iP,k 



-1)"^^- Nk\ nil 



1 

(20! 



(113) 



This result is virtually identical to Eq. (I103p except that there are no "+l's" 
in the denominators of the above equation. 

More sophisticated results involving both the secant and cosecant num- 
bers can also be derived. From No. 1.518(2) of Ref. [23] we have 



lnsec(7rz) = N^ 



fc=i 



(22fc _ 1) 

k{2k)\ 



-)2fc-l 



B2k\{T^Z 



i2fc 



(114) 



which is absolutely convergent for \z\ < 1/2. In Ref. ^ it is shown that the 
cosecant numbers are related to the Bernoulli numbers by 

c^ = ^-^ (2^' - 2) B,, . (115) 

Moreover, we can write the Ihs of Eq. (I114p as 

oo 

In sec(7rz) = In f 1 + N^ dk 
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We now apply Theorem 1 to the rhs of the above result. This means that we 
expand the logarithm in terms of its Maclaurin series expansion, in which case 
qk = (— l)'^^^/fc, while y = z^ and pk = dk'n'^''. We then equate the resulting 
power series expansion to like powers of z'^ or y in Eq. (I114p . By substituting 
Eq. (I115p we replace the Bernoulli numbers by the secant numbers, thereby 
obtaining 
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Once again, we see the partition operator acting on another strange argument 
to yield an interesting finite quantity for all values of k. 

In Ref. |1] the partition method for a power series expansion is applied 
to the reciprocal of the logarithmic function ln(l + z). There a power series 
expansion is obtained in terms of special coefficients Ak, which are referred 
to as the reciprocal logarithm numbers. On p. 138 of Ref. {2H] these numbers 
are referred to as the Gregory or Cauchy numbers when their modulus is 
taken. That is, the following result is obtained: 



ln(l + z) 



E-4 

fc=0 



kZ' 



k-1 



(118) 



The reciprocal logarithm numbers are found to be oscillating fractions, which 
are more slowly converging to zero than either the cosecant or secant num- 
bers. Moreover, they are given by 
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Expressing ln(l+-2;) in terms of its Maclaurin series as in the previous example, 
which is absolutely convergent only for |-2| < 1, we are in a position to apply 
Theorem 1 to z/\n{l+z). In this case, h{z) = l/z and f{z) = ln{l+z). Then 
the coefficients of the inner series, viz. p^, are equal to {—l)^^^/{k + 1) for 
k > 0, while for k = 0, po = 0. The resulting denominator can be regarded 
as the regularised value of the geometric series, which means that the q^ are 
equal to (—1)'^ as in the preceding examples. Again, a = l. Hence, according 
to Eq. ( TTOl) . the reciprocal logarithm numbers can be written as 
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The inverse of this result is obtained by putting Dq = 1 and D^ = A^ in Eq. 
( !T2|) . while the E^ equal the coefficients in the Maclaurin series for ln(l+z), 
i.e. (-1)^'+V(A; + 1). Then we find that 
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As an aside, in Ref. [T], Euler's constant is derived in terms of an infinite 
series involving the reciprocal logarithm numbers, where it is also referred 
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to as Hurst's formula. Since the publication of Ref. pQ, it has been revealed 
that the formula was independently discovered by Kluyver [5S]. By using 
Eq. fll20p . we can express Euler's constant as 



7 



Yl h ^P'k 



k=l 



;-i)^iv!n 



i=l 



t + l 



1 

riil 



(122) 



Alternatively, we can introduce Eq. fll2ip into Hurst's formula, which yields 
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Moreover, Euler's constant is not the only result found in Ref. pQ that can 
be expressed as an infinite sum over the reciprocal logarithm numbers. For 
example. In 2 can be expressed as a similar sum to Hurst's formula. Therefore, 
with the aid of Eq. (I120p we find that 
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With the aid of Corollary 1 to Theorem 1 we can generalise the preceding 
examples to where the generating functions are raised to an arbitrary power 
p. For example, the quotient in Eq. ( FM|) raised to an arbitrary power p 
becomes 
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where we have used the regularised value for the binomial series in Lemma 
1. On the other hand, according to Equivalence ([S3D, the Ihs of the above 
result can be written as 
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Hence, Dq = 1 and Dk = 2^ ^b'^ for k > 1, while y = z. From Eq. (1601) we 
obtain 
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Equating like powers of z on the rhs's of the preceding equivalences yields 
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If Equivalence fl99|) is taken to the arbitrary power of p, then the D^ of 
Equivalence fl63|) become the cosecant numbers, i.e. Dj = Cj. If we denote 
-Dfc(p) by Cp^}^, then we find that according to Eq. fIBUl) . these generalised 
cosecant numbers are given by 
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Alternatively, the original equation can be expressed as 
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Now the Dfc in Equivalence ( l63l) are equal to (— 1)'^/(2A; + 1)!, while p has 
changed sign. Thus, the generalised cosecant numbers can be written as 
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This form for the generalised cosecant numbers has been recently been em- 
ployed in evaluating a finite sum of inverse powers of cosines given as Sm,v = 
(7r/2m)2'' YIZi cos(A;7r/m)-2'^ in Ref. |4]. 

In a similar manner we can generalise the secant numbers. By taking the 
p-th power of Equivalence (11081) . we see that the Dk in Equivalence ( l63l) are 
equal to di. By denoting the generalised secant numbers as dp^k, we find via 
Eq. ([60D that 
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while taking the p-th power of Eq. (11071) yields 
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The above result appears as Eq. (290) in Ref. [5]. 

To generalise the reciprocal logarithm numbers, we take the p-th power 
of Equivalence f lllQp . Then the Dk in Equivalence fl63|) are equal to Ak. 
Denoting the generalised reciprocal numbers by Ak{—p), we find via Eq. f l60|) 
that they are given by 

r ^ 4"'i 

A,{-p) = Lp^u {-If {-p)n n -^ ■ (134) 

1=1 

Alternatively, the generalised reciprocal logarithm numbers can be deter- 
mined by taking the p-th power of the Maclaurin series for ln{l + z) , which 
is 

ln(l + z) = Y, ^ ^'^ • (135) 
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As explained in Ref. |2], the above result is absolutely convergent for ^z<—l, 
in which case the equivalence symbol can be replaced by an equals sign. With 
regard to Equivalence fl^ the D^ are now equal to (— l)^"*"^//?, while y = z. 
Consequently, Eq. (160|1 yields 

k 
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Eq. flT36ll appears as Eq. (118) in Ref. [2]. 

At the end of Sec. 2 an alternative algorithm was given for accessing the 
partitions via the tree diagram in Fig. [1] As a result, we can present a new 
formulation of the partition method for a power series expansion. Before 
doing so, however, we need to amend the definition of the partition operator 
Lpfc. The amendment is necessary so that we can sum along each of the 
branches emanating from the seed number. For example, the entire sub-tree 
from {1,5} represents the tree diagram for partitions summing to 5, but to 
obtain those summing to 6, we need to increment ni in all the partitions 
summing to 5 by one. This, of course, affects all the contributions to the 
coefficients. In addition, when we consider the partitions emanating from 
{2,4}, we need to ensure that no partitions with unity will appear, while 
there should be no ones or twos for the partitions emanating from {3,3} and 
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so on. Therefore, we define the restricted partition operator as follows: 
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There are major differences between the above operator and the partition 
operator as defined by Eq. fl72|) . The first is that the sum begins at n^ rather 
than at n\. This is due to the fact that we need to exclude elements less than 
i when summing along the branches emanating from the seed number. The 
next difference is that Ui begins at unity rather than zero as for the other 
elements, which accounts for the fact already a one element of i has occurred 
in moving from the seed number. Furthermore, as a result of separating the 
element «, the maximum element in the resulting partitions becomes k — i. 
Hence, the elements range from i to k — i, which not only affects the number 
of summations, but also their upper limits. In addition, whilst the number of 
summations is restricted in the constraint, the value remains invariant, viz. 
k. 

We are now in a position to implement the algorithm described at the end 
of Sec. 2. Basically, this entails expressing a result like Eq. (TTUj) in terms of 
the partition operator on the Ihs and the sum of restricted partition operators 
on the rhs. Therefore, we arrive at 
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where N, 



i,k- 



Ek- 
j= 



rin 



and Ni^k = N. For po 7^ in Theorem 1, we use 



Eq. IQ instead, which amounts to replacing qiapk and Qn ,,_■ a^^-'^-J Nj^k-r 
in the above result by a F^''\apo) and a^J''-— ^ F^^-'^-^ (apo) respectively. 

To complete this section, we consider the situation where the quotient of 
the pseudo-composite functions in Theorem 1 yields a function r{y), which 
is infinitely differentiable. This results in the following corollary. 
Corollary 2 to Theorem 1. If the functions f{z) and g{z) obey the same 
conditions as in Theorem 1 and the quotient of the pseudo-composite func- 
tions ga° f and ha o f yields an infinitely differentiable function, r{y), then 
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for po 7^0, 
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while for po = 0, 
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Furthermore, if r{0)^0, then inversion of the quotient yields 
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where the E^ are given by Eq. flT2|) . 

Proof. From the proof of Theorem 1, we know that the ratio oi Qa o f 
over ha o f is equivalent to the power series YlT=o E^kV^, where the coefficients 
Dk are given by either Eq. ([9]) for po non-zero, or Eq. ( TTOl) when po = 0- Since 
the ratio of the pseudo-composite functions yields an infinitely differentiable 
function r{y) according to the corollary, the ratio can also be expressed as a 
Taylor/Maclaurin series given by 



9a ° f 
haO f 



E^'"(°)^ 



fc=0 



(142) 



where the superscript (k) now denotes the k times differentiation of r{y) 
w.r.t. y. From Theorem 1 we also know that the above quotient can be 
expressed as 
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where the coefficients Dk are given by either Eq. ([9]) when po 7^ or Eq. ( TTOll 
for po = 0. Since the regularised value is unique as described in Refs. [2j and 
[7j-|10], the rhs's of the two preceding equivalence statements are equal to 
one another. Since y is arbitrary in the resulting equation, we can equate 
like powers or the coefficients of both power series, which yields Eq. (I139P or 
014Up depending upon the the value of po- 
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If the quotient of the pseudo-composite functions is inverted, then it will 
equal l/r{y). If r(0)7^0, then the function l/r{y) can also be expressed as a 
Taylor/Maclaurin series since r{y) is infinitely differentiable. Therefore, we 
have 
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From Theorem 1 we know that the above quotient also represents the reg- 
ularised value of a power series in y whose coefficients are equal to Ek/Do, 
while the E^ are given by Eq. flT2l) . Moreover, since r(0) 7^0, Dq does not 
vanish. Since both power series have the same regularised value, they are 
equal to one another for the same reason as in the first part of the proof. 
Again, as y is arbitrary, we can equate like powers of y, which results in Eq. 
(11411) . This completes the proof of the corollary. 

So far, we have described the partition method for a power series expan- 
sion in terms of a novel discrete operator, which has been referred to as the 
partition operator and is denoted by -Lp,fc[-]. At this stage the operator has 
been used to derive general results for the coefficients of the power series 
expansions obtained via Theorem 1. However, it is not always possible to 
derive general results dependent only upon k for the coefficients D^ and Ek 
in Theorem 1. Often the coefficients become polynomials dependent upon 
another variable or even a function. For these cases we need to develop a pro- 
gram that can automatically evaluate specific coefficients, which is addressed 
in the following section. 

4 Programming the Partition Method for a 
Power Series Expansion 

At the end of the previous section we indicated the need for developing a 
computer program to enable the evaluation of the coefficients via the par- 
tition method for a power series expansion. As indicated earlier, such an 
approach can be developed by employing the BRCP algorithm of Sec. 2, 
but before doing so here, some remarks are necessary. Because the partition 
method for a power series expansion relies on evaluating the contribution due 
to each partition and the number of partitions p{k) according to the Hardy- 
Rademacher-Ramanujan formula is 0{exp{^T^y2k/3)/4:k^/3) jIHl HH [30] , any 
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computer program based on partitions as its input will ultimately become 
very slow. In fact, since all the partitions summing to the order of each power 
are involved, such a program represents a brute-force approach to deriving 
power series expansions. Nevertheless, determining power series expansions 
for orders up to 40 (p(50) = 37338) should be achievable with most number- 
crunching computers around today. So, at least for intermediate values of 
the order k, programming the partition method is still of great benefit, par- 
ticularly for intractable functions where it represents the only method we 
have of deriving a power series expansion. 

With regard to very high orders it should be noted that the partition 
method does not actually use the partitions themselves. What the method re- 
quires is each element appearing in each partition and their frequency, which 
is referred to as the multiplicity representation in Sec. 2. This information 
can be stored in external arrays which can be called upon when one wishes 
to determine the series expansion for different situations. Therefore, there is 
no need to repeat the process of generating the partitions when dealing with 
different problems. In addition, the contributions due to many of the parti- 
tions will often be negligible even by today's computing standards. In those 
cases the calculation of the coefficients can be simplified yielding extremely 
accurate approximations. In other cases it is possible to sum the contribu- 
tions in classes or groups, thereby avoiding the necessity of processing each 
partition separately. This issue will be addressed later in this work. Finally, 
by developing a programming approach to the partition method, we will be 
in a position to consider different problems in the theory of partitions such as 
the evaluation of partitions with specific elements including those with dis- 
crete elements, doubly-restricted partitions and the transposes of partitions. 
The significant issue here is that such problems only require small changes 
to the BRCP algorithm, whereas separate programs are required when other 
codes such as those presented in Sec. 2 are used to generate partitions. E.g., 
in order to determine the partitions with a fixed number of elements in them, 
Knuth presents another algorithm based on the 18-th century dissertation by 
CF. Hindenburg on p. 38 of Ref. [IB]. As we shall see in Sec. 5, this problem 
can be solved by inserting a few lines into the BRCP algorithm. 

In the previous section the partition method for a power series expansion 
was described in terms of the partition operator, Lp^k[-]- There many general 
results involving this operator were given without requiring to evaluate the 
sum over partitions for specific values of k. Consequently, the partition oper- 
ator can be regarded as an intricate abstract operator when compared with 
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more well-known operators such as the differential operator. However, whilst 
taking the derivative of a function is also viewed as an abstract operation, we 
at least have an understanding of the process because we can always calculate 
the limit of Newton's difference quotient provided, of course, it exists. As a 
result of this understanding, general shorthand rules such as dx^/dx = kx^~^ 
have evolved. Yet, the opposite situation applies to the partition operator — 
we have a few general results in the previous section and in Refs. [T]- [3], but 
we do not even have a means of applying the operator outside of those cases 
to evaluate the first few coefficients of the expansions given in Theorem 1. 
Therefore, we require an approach that will allow us to apply the partition 
operator for any value of k to any situation that obeys the conditions in 
Theorem 1, even if it is no longer feasible to evaluate the coefficients for very 
large values of k. 

Now that we have indicated why it is necessary to program the partition 
method for a power series expansion, we turn to the issue of the program- 
ming languages required for the task. The first point to be noted is that if we 
choose a standard high-level programming language like C/C++ or Fortran, 
then our results for the coefficients will inevitably become decimal numbers 
when they could be rational. Moreover, they will invariably be rounded off 
or worse still, may only equal zero if significantly smaller than the preci- 
sion allowed by the computing system. In addition, the coefficients need not 
be numerical as exemplified by the examples appearing after Corollary 1 to 
Theorem 1. All this means is that we require a mathematical software pack- 
age such as Mathematica to retain either the rationality of the coefficients 
or when applicable, their symbolic form. However, programming the BRCP 
algorithm in Mathematica with its bi-variate recursion is also formidable. In- 
stead, the issue can be overcome by using the material in Sec. 2. Therefore, 
the best option is to combine the strengths of both C/C++ and Mathemat- 
ica. Basically, this means that the initial program is to be written in C/C++ 
so that the coefficients can be printed out in a symbolic form. Then these 
forms can be introduced into Mathematica where we can use either the in- 
teger arithmetic routines to evaluate the coefficients, thereby avoiding the 
round-off that occurs with floating point numbers or its symbolic routines to 
reduce all the terms generated by the C/C++ code into simple mathematical 
statements such as polynomials. 

The appendix presents the C/C++ program called partmeth, which 
outputs in symbolic form the coefficients Dj. and E^ given in Theorem 1. 
Here, we are only concerned with the case of po = or Eq. (jH]) for the Dk, 
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while the E^ are given by Eq. ( 1T^ . The case of po 7^ or Eq. (jS]) is left 
as an exercise for the reader. If we compare the code with the final code 
in Sec. 2, then we see that the overall structure remains the same. That is, 
there is a main section with the same two function prototypes termgen and 
idx. In fact, idx or the BRCP algorithm has not been altered at all, but 
termgen and the main routine have been changed to produce the symbolic 
forms for the coefficients in the partition method for a power series expansion. 
Besides evaluating the execution time, main carries out the calculation of the 
coefficients in one for loop, which is limited by the variable dim, representing 
the maximum value of k or the coefficient of the highest order term which the 
user must input. Within this for loop there are two calls to idx, one of which 
applies to the calculation of the D^ and the other to the calculation of the Ej.. 
Therefore, it is termgen that is doing the heavy work in the programme. 
In fact, we shall see in the next section that by modifying termgen we can 
determine many of the properties of partitions, which often require separate 
programs. 

Within termgen we see that the Dk and E^, which are represented by the 
variables DS[k,n] and ES[k,n] respectively, are evaluated depending upon 
the value of the variable inv-case. If it equals zero, then the Dk are evalu- 
ated, while if it equals unity, then the Ek are evaluated. In evaluating the 
latter there is also an extra complication due to the phase factor of (—1)^ in 
Eq. fll2p . Consequently, for this case termgen must determine the number 
of distinct elements in each partition. When dim = 4, partmeth prints out 
the first four values of the Ek and Dk in symbolic form. E.g., the k = A values 
that it prints out are : 

DS[4,n_]:= p[4,n] q[l] a + p[l,n] p[3,n] q[2] a^(2) 2! 

+ p[l,n]^(2) p[2,n] q[3] a^(3) 3!/2! + p[l,n]^(4) q[4] a^(4) 

+ p[2,n]^(2) q[2] a^(2) 

ES[4,n_]:= -DS[0,0]^(-2) DS[4,n] + DS[0,0]^(-3) DS[l,n] DS[3,n] 2! 

- DS[0,0]^(-4) DS[l,n]^(2) DS[2,n] 3!/2! + DS[0,0]^(-5) DS[l,n]^(4) 

+ DS[0,0]^(-3) DS[2,n]^(2) 
From these results we see that each coefficient is composed of five distinct 
terms corresponding to the fact that the number partitions summing to 4, 
i.e. p(4), is equal to 5. These results allow for the situation where the p[k, n] 
may be dependent upon another variable, viz. n, even though it may not be 
necessary. 

The first code presented in the appendix is suitable for values up to 
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and around k = 20. In fact, all the values of D/^ and Ej. for k < 20 are 
computed within one CPU second. For k > 20, however, the expressions 
become unwieldy and thus, it is better to evaluate them separately so that 
each can be introduced directly into Mathematica. This amounts to removing 
the for loop in main and computing only for the value of k or dim. The code 
for computing the D^, which is called mathpm, appears immediately after 
partmeth in the appendix. Furthermore, in order that the coefficients can 
be introduced directly in Mathematica, only three terms appear on each line 
of output, while a plus sign now appears as the last character on each line 
except, of course, on the final line. 

Let us now consider the evaluation of D^q via mathpm. Since p(30) = 
5604, this is the number of distinct terms when mathpm prints out DS[30, n]. 
Even though the output file for DS[30, n] is very large, it can still be imported 
into Mathematica. If we now set pfc = (— l)^/(2/c + 1)!, qk = l{—l)'' and a = l, 
which represent the inner and outer series for the cosecant numbers, then it 
takes 0.15 CPU sec to evaluate c[30] or C30 in integer form on the same Sony 
VAIO laptop mentioned in Sec. 2. In this instance the numerator is given by 
a 60 digit integer, while the denominator is given by a 90 digit integer. In 
decimal notation the value of C30 is approximately 2.965 x 10~^°. If we use 
Eq. (llOOp to evaluate C30, then we find that it takes almost zero CPU sec to 
evaluate the same result. On the other hand, if we set pk = {—l)^/(2k)\, i.e. 
the situation for the secant numbers d^, then we find that it takes 0.14 CPU 
sec to evaluate d^o in integer form on the Sony VAIO laptop. In this case the 
numerator and denominator are respectively 67 and 78 digit numbers, while 
in decimal form d^o is approximately equal to 2.176 x 10^^^. Unfortunately, if 
Eq. flllOp is implemented in Mathematica, then we only obtain approximate 
values in decimal form for the secant numbers. Hence, we need to implement 
a recurrence relation such as Eq. flll2p in order to obtain them in integer 
form. When this is done, it is found that Mathematica takes 6548 CPU sec 
to compute d^o- 

If we set pk = {—l)^/{k + 1), which represents the situation for the recip- 
rocal logarithm numbers A^, then we find that it takes only 0.1 CPU sec to 
determine A^q in integer form. In this instance the numerator and denomi- 
nator are 35 and 38 digit numbers yielding an approximate decimal value of 
1.474 X 10~^, the slowest converging of the numbers considered so far. The 
reciprocal logarithm numbers can be evaluated by either relating them to the 
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Stirling numbers of the first kind [T] via 
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or by the recurrence relation of 

If the first form is implemented in Mathematica, then it takes 0.1 CPU sec to 
compute ^30, while with the second form it takes 5719 CPU sec. Therefore, 
we see that the evaluating the coefficients of power series expansions via the 
partition method for a power series expansion can be vastly superior to us- 
ing recurrence relations and is almost on a par with the cases where intrinsic 
forms have already been implemented within a mathematical software pack- 
age. Furthermore, by altering the relations for the coefficients of the inner 
and outer series in addition to a, we obtain results for other mathematical 
quantities with different power series expansions. 

As discussed previously, the coefficients of the inner and outer series do 
not need to be numbers as has been the case so far. If we set qk = (p) — k/k\, 
Pk = (— 1)'^^^/(2A;+1)! and a=l, then DS [30, n] yields the generalised cosecant 
number Cp^so as given by Eqs. (11291) and (11311) . Then we find that it only takes 
0.36 CPU sec to compute the resulting polynomial, which is thirtieth order 
in p. By invoking the Simplify routine in Mathematica the polynomial can 
be arranged in increasing order within another 0.36 CPU sec. Before this 
calculation was performed, the results produced by DS[5,n] and DS [8, n] had 
been found to agree with the generalised cosecant numbers, Cp^s and Cp^g, 
obtained in Ref. ^. 

As stated earlier, the calculation of the coefficients via the forms generated 
by either of the first two programs in the appendix can be continued beyond 
the thirtieth order, but eventually problems arise due to the combinatorial 
explosion occurring in the number of partitions. It has already been stated 
that there are 190 569 272 partitions summing to 100, which means that this 
number of terms will be present in DS[100,n]. If mathpm is run for this case, 
then it takes around 600 CPU sec to compute DS[100,n]. Whilst this is not an 
overly long time of computation compared with the earlier results obtained 
via recurrence relations, it produces a file whose size is over 16 GB. Files of 
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this size are going to present a problem when imported into mathematical 
software packages. For example, it appears that Mathematica 8.0.1 is only 
able to import files with 2Gb of data. One method of circumventing this 
problem would be to divide the file into smaller files so that could be handled 
by the different processors on a supercomputer. Then the results generated 
by each processor could combined to yield the final answer. 

Another method of overcoming this problem is to introduce the values 
for pk, qk and a first and then evaluate a specific number or limit of terms 
via mathpm. Once the limit point is reached, the values where this occurs 
would need to be stored. In terms of Fig. [T]this amounts to storing the values 
of both arguments in idx. Then the partial value of the coefficient could be 
evaluated and stored, while all the terms outputted in running mathpm can 
be either deleted or overwritten in a re-run of mathpm. In the re-run of 
mathpm the code would not print out any terms until idx reaches the values 
of the arguments of idx stored from the first run. Then the code would either 
continue to print out the next limit of terms stopping at two new values of 
idx or would terminate on reaching the central partition. Then the terms 
stored in the second run could be evaluated and combined with the result 
obtained from the first run. If the central partition has not been reached in 
the second run, then the process can be continued until the central partition 
is eventually reached. Of course, the disadvantage in this approach is that 
we have lost the ability to evaluate a new coefficient by altering the pk, qk 
and a as we were able to do with DS[30,n] above. This second method of 
solving the problem of very large data files produced by running mathpmis 
contingent on whether we can stop and re-start the program at specific points 
in the tree diagrams for the partitions. Hence, we need to be able to adapt 
the BRCP algorithm so that specific partitions can be determined, which 
represents the topic of the following section. 

5 Specific Types of Partitions 

In this section we aim to investigate how the BRCP algorithm presented in 
Sec. 2 can be modified to determine specific types of partitions. By specific 
partitions we mean such general problems in partitions as the determination 
of: (1) partitions with a fixed number of elements, (2) doubly-restricted 
partitions, (3) discrete partitions, (4) conjugate partitions and (5) partitions 
with specific elements in them. Solving such problems invariably means 
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creating a different algorithm or program for each problem as can be seen 
in Refs. [TB], [TS] and [22]. However, we shall see here that such problems 
and others can be solved with relatively minor modifications to the BRCP 
algorithm presented in Sec. 2, once again highlighting its versatility. As we 
shall be modifying partgen in Sec. 2 when solving these problems, we shall 
generate the partitions in the compact multiplicity representation, although 
it should be mentioned that the various algorithms presented in this section 
would only require minor modification to termgen to generate partitions in 
the standard representation. 

Of all the problems mentioned in the previous paragraph perhaps the 
simplest one to consider is the determination of those partitions with a spe- 
cific element or elements in them. E.g., suppose we wish to determine all 
those partitions summing to 15 with the element of {5} in them? From the 
material presented in Sec. 2, it is obvious that the total number of partitions 
for this problem must equal the number of partitions summing to 10 or p{10). 
Moreover, the partitions generated by the new code should yield the same 
partitions obtained by running the various codes presented in Sec. 2 except 
that each partition generated by the new code will have an additional ele- 
ment of {5} in them. We shall see that this is indeed the case, although the 
order in which the partitions are generated is different from those discussed 
in Sec. 2. 

In order to modify partgen in Sec. 2 so that it generates partitions with 
a specific element in them, we first need to alter the main prototype of the 
program. This is required to enable users to input the specific values of the 
elements that they wish to appear in the partitions printed out by the new 
program. Consequently, main becomes 

int main( int argc , char *argv[] ) 

{ 

int i ; 

if(argc! = 3) printf (" usage : specpart <partition sum> 

<sp_val> \n" ) ; 
else{ 

tot=atoi ( argv [ 1 ] ) ; 

sp_val=atoi ( argv [2 ] ) ; 

part = (int *) malloc ( tot * sizeof ( int ) ) ; 

if (part ^= NULL) printf (" unable to allocate 
array\n" ) ; 
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else { 



} 



for ( i=0;i<tot ; i++) part[i]=0; 
idx ( tot , 1 ) ; 
free ( part ) ; 



} 

printf("\n"); 

return (0) ; 

} 

Here we see that the new program called specpart has a global variable 
called sp-val, which represents the specific element that is to appear at least 
once in each partition printed out by the program. 

As discussed in Sec. 2, idx(tot,l) scans over all the partitions summing 
to tot, while termgen, which is called in idx, is responsible for generating or 
printing out partitions. Thus, in order to determine those partitions with a 
specific element or elements in them, we need to modify termgen. That is, 
we still need to scan over all partitions by calling idx(tot,l). In fact, aside 
from making minor modifications to main, we shall find that to solve all the 
problems mentioned in the introduction to this section, we need only modify 
termgen. 

In Sec. 2 the termgen function prototype in partgen was responsible 
for printing out all the partitions summing to tot in the multiplicity repre- 
sentation. This was achieved by processing the array part, which stored the 
frequencies of the elements in each partition. That is, the variable freq was 
used to represent the frequency of the element i + 1 in the partition with i 
ranging from to tot — 1 . Now that we wish to determine those partitions 
with a specific element or elements in them, we need to restrict the partitions 
that are printed out by termgen. This is accomplished simply by introduc- 
ing a local variable called freq^spval, which evaluates the frequency of sp^val 
in each partition. If this value is non-zero, then we know that there is at least 
one occurrence of sp-val in the partition and the partition is then printed out 
in the same manner as in Sec. 2. li freq^spval is zero, then the partition is 
ignored. Therefore, the termgen function prototype for specpart becomes 

void termgen () 

{ 

int freq , i , freq_spval ; 
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freq_spval=part [sp_val —1]; 
if ( freq_spval ) { 

printf ("%ld : " ,term ++); 
for (i=0;i<tot ; i++){ 
freq=part [ i ] ; 
if(freq) printf ("%i (%i 



freq , i +1^ 



printf ("\n' 
} 



} 



} 

In Sec. 2 tlie partitions summing to 5, which amounted to 7, were printed 
out by running partgen. The output produced by running specpart with 
tot and sp^val set equal to 11 and 6, respectively, is: 
1: 5(1) 1(6) 
2: 3(1) 1(2) 1(6) 
3: 2(1) 1(3) 1(6) 
4: 1(1) 2(2) 1(6) 
5: 1(1) 1(4) 1(6) 
6: 1(2) 1(3) 1(6) 
7: 1(5) 1(6) 

Hence, we observe that the number of partitions is once again 7 or p(5). If 
we remove one six from each partition, then we obtain the same partitions 
as those generated in Sec. 2 by partgen except that the order in which they 
appear is now different. If we define the specific element partition operator 
LsEP,k,j[-] as 



^SEP,k,j 



k,...,[k/{j-i)],[k/jUk/(j+i)]...,i 



E 



(147) 



ni,...,nj_i=0,nj=l,n 



■i+i' 



,,rafc=0 



^^=l,^7^j •"i^'^-J 



then it follows that 



^SEP,k,j 



'P,k-j 



p{k - j] 



(148) 



As a result of the preceding analysis, it becomes a simple matter to con- 
sider partitions with more than one specific element in them. All we need 
to do is introduce more values in main and then create local variables like 
freq^spval to represent the number of occurrences or frequencies for each of 
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these values. As before, each frequency is set equal to part[sp^val — i] in 
termgen. Next, we modify the if statement to include all the frequencies. 
For example, if we wish to determine all the partitions with two specific 
elements, viz. sp-val and sp-val2, then the if statement becomes 

if (freq_spval &fe freq_spval2 ) { etc. 

On the other hand, we may want those partitions with at least one occurrence 
of either sp-val or sp-val2. In this instance all we need to do is replace the 
logical AND operator by the logical OR operator or && in the if statement. 

By making a few minor modifications, mostly to termgen, we have been 
able to solve several problems involving specific partitions. Now we consider 
evaluating the partitions with a fixed number of elements in them, which 
has already been considered when we modified the partition operator into 
the form given by Eq. flHHj) . Previously, it was stated that to solve this 
problem, often a completely different approach is employed. For example, 
Knuth presents a different algorithm on p. 38 of Ref. jTS] compared with the 
algorithm used to generate partitions in reverse lexicographic form. In the 
case of the BRCP algorithm, however, the number of elements in a partition 
is determined by the number of branches along its path prior to termination 
in the tree diagram. E.g., the partitions with only two elements in them 
in Fig. [1] are obtained by searching for the terminating tuples appearing in 
the third column, viz. (0,5), (0,4) and (0,3). In other words, these are the 
terminating tuples two branches away from the seed number situated in the 
first column. When we include the elements from the second column, they 
yield the partitions {1,5}, {2,4} and {3,3}. Since partitions with a fixed 
number of elements can be determined from the tree diagram, it means that 
only minor modifications to the program in Sec. 2 are again required in order 
to generate these partitions. 

As in the previous example, we need to modify main so that the user can 
input the number of parts, which will be represented by the global variable 
numparts. Then we can proceed to the modification of termgen, which is 
presented below. There we see that a new local variable called sumparts has 
been introduced. This variable determines the number of elements in each 
partition. When this number equals numparts, the partition is printed out. 
Otherwise, the partition is discarded. 

void termgen () 

{ 

int freq ,i ,sumparts = 0; 
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/* sumparts is the number of parts or elements in a 
p artition */ 



for (i=0;i<tot ; i++){ 

sumparts= sumparts+part [ i ] ; 

} 
if(sumparts = numparts){ 
printf("%ld: " ,term ++); 
for (i=0;i<tot ; i++){ 

freq=part [ i ] ; 

if(freq) pr intf ("%i (%i ) ",freq,i+l); 

} 
printf ("\n" ); 

} 

} 

When the code is run for the case where the partitions sum to 10 and 
possess 5 elements, i.e. for toi = 10 and numparts = 5, the following output is 
produced: 



1 


4(1) 1(6) 


2 


3(1) 1(2) 1(5) 


3 


3(1) 1(3) 1(4) 


4 


2(1) 2(2) 1(4) 


5 


2(1) 1(2) 2(3) 


6 


1(1) 3(2) 1(3) 


7 


5(2) 


T 


herefore, we see that the number of partitions in this instance is 7, which 



can also be expressed as either 



10 

5 



according to Sec. 2 or as Lp]^q[1] according 



to Eq. (|90|). 

It should also be mentioned that if the condition, 

if ( sumparts^=numparts) { etc., 

in the above program is replaced by 

if ( sumparts<=numparts) { etc., 

then the resulting code generates all those partitions summing to tot with 
at most numparts elements. E.g., the number of partitions summing to 10 
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in Sec. 2, this is also equivalent to 



However, the partitions generated 



with at most 5 elements, which according to Sec. 2 is equal to 30. As stated 

15 
5 

by the code given above for tot and numparts equal to 15 and 5 respectively, 
are different from those generated by the version of the program with the 
modified if statement. 

Doubly-restricted partitions are those partitions where all the elements 
are greater than a particular value and less than another value. Since this is 
a combination of two separate conditions, first we need to be able to modify 
partgen so that it generates the partitions where the elements are either 
greater than or lower than a specified value. Therefore, let us consider the 
situation where all the elements in the partitions are less than or equal to 
a value, which will be represented by the global variable largest^elt. As in 
the previous examples, this value must be introduced into main. Then we 
proceed to the modification of termgen. 

To modify termgen so that only partitions with elements less than or 
equal to largest- elt are generated, we need to introduce an extra for loop. 
This extra loop is required so that if a partition is encountered where an 
element is greater than largest- elt, it is discarded via a goto statement as 
demonstrated by the modified version of termgen given below. Although 
goto statements are generally frowned upon by programmers, it is being 
used here to abandon processing in a nested structure of two loops. In fact, 
the code behaves much like partgen in Sec. 2 when all the elements in the 
partitions are less than or equal to largest_elt. However, when an element 
is greater than largest_elt, the goto statement discards the partition by 
diverting to the end statement label. 

void termgen () 

{ 

int f , i ; 

for(i = largest_elt ; i<dim ; i++){ 

f=part [ i ] ; 

if(f) 

goto end ; 

} 
printf("%ld: " ,term ++); 

for ( i =0; i<dim ; i++){ 

f=part [ i ] ; 
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if(f) printf("%i(%i) " ,f ,i+l); 

} 
printf("\n"); 

end : ; 

} 

When the above code is run for partitions summing to 14 in which the ele- 
ments are greater than 3, the following output is produced: 
1: 1(14) 
2: 1(4) 1(10) 
3: 2(4) 1(6) 
4: 1(4) 2(5) 
5: 1(5) 1(9) 
6: 1(6) 1(8) 
7: 2(7) 

Hence, we see that there only seven partitions summing to 14, in which all 
the elements are greater than 3. 

As a result of the above code, it is now a simple matter to consider the case 
where all the elements are greater than or equal to another value specified 
by the user. In this instance we simply replace largest- elt by another global 
variable called smallest^ elt and alter the condition in the first for loop of 
the previous version of termgen. That is, the first for loop in the preceding 
version of termgen simply becomes 

for ( i =0; i<smallest _elt — 1; i++){ etc. 

Note also that if the upper value in the for loop had been smallest- elt instead 
of smallest-elt — 1, then all the elements generated by the code would only 
be greater than the value specified by the user. 

For doubly-restricted partitions, where the elements are greater or equal 
to value and less than or equal to another (larger) value, all we need to do 
is incorporate two for loops that divert to the end statement label. For 
example, termgen for this situation would become 

void termgen 

{ 

int f , i ; 

for ( i =0; i<smallestpart — 1; i++){ 

f=part [ i ] ; 

if ff) 
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;oto end ; 



-){ 



} 



for(i = largestpart ;i <dini ; i 
f=part [ i ] ; 
if(f) 

goto end ; 

printf("%ld: " ,term ++) 
for ( i =0; i<dim ; i-H- ){ 

f=part [ i ] ; 

if(f) printf("%i(%i 

} 

\n" ) ; 



f,i+l): 



printf 
end : ; 

} 

When the above code is run for partitions summing to 13 with the ele- 
ments greater than or equal to 3 and less than or equal to 9, the following 
output is produced: 
1: 2(3) 1(7) 
2: 3(3) 1(4) 
3: 1(3) 1(4) 1(6) 
4: 1(3) 2(5) 
5: 1(4) 1(9) 
6: 2(4) 1(5) 
7: 1(5) 1(8) 
8: 1(6) 1(7) 

Thus, we see that there are 8 partitions with all elements lying in the interval 
[3,9]. 

In Ch. 3 of Ref. [21] Andrews defines restricted partitions differently from 
the earlier definition given by Eq. (11371) . There, the partitions represent those 
in which the elements are less than a value, say eLmax, while the number 
of elements is less than or equal to another value, which we take again to be 
numparts as in the preceding examples. Although there is now a condition 
pertaining to the number of parts, evaluating these partitions is again similar 
to the doubly-restricted case studied above. First, we must introduce eLmax 
into main in addition to numparts. Then we need to insert an extra for loop 
into termgen so that it can make use of the different condition. The new 



loop appears first since if it is true, we immediately by-pass any action to 
process the current partition. Therefore, this modified version of termgen 
becomes 

void termgen () 

{ 

int freq , i , sumparts = 0; 

for( i=0;i<tot ; i++){ 

if(i> el_max— 1 &fe part[i]>0) goto end; 

} 
/* ( 1 ) sumparts is the number of elements in a p artition 

(2) all elements are now less than or equal to eLmax 

*/ 

for( i=0;i<tot ; i++){ 

sumparts= sumparts+part [ i ] ; 

} 
if (sumparts <= numparts){ 
printf ("%ld: " ,term++); 
for (i=0;i<tot ; i++){ 

freq=part [ i ] ; 

if (freq) pr int f ("%i (%i ) ",freq,i+l); 

} 
printf("\n"); 

} 

end : ; 

} 

When the code is run with tot, numparts and el^max set equal to 10, 3 and 

5 respectively, the following output is produced: 

1: 1(1) 1(4) 1(5) 

2: 1(2) 1(3) 1(5) 

3: 1(2) 2(4) 

4: 2(3) 1(4) 

5: 2(5) 

Hence, we see that there are 5 partitions summing to 10 with at most 3 

elements and all elements less than or equal to 5. 

The number of partitions summing to k with at most M parts and each 
element less than or equal to A^ is represented as pG{N,M,k) in Ref. 
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The subscript G has been introduced here so that the reader is not con- 
fused with similar notation in the next section. From the above example 
we have pg(5,3, 10) = 5. li k > MN, then pg{N, M,k) vanishes, while 
Pg{N,M,NM) = 1. These numbers also appear as the coefficients in the 
generating function for Gaussian polynomials, which are given by 

^ n _ a^+' ) ^^ 

G(iV,M;g) = n \^_ ,/ =l + J2pG{N,M,k)q'' . (149) 

j=i ^ ^ ^ k=i 

Hence, Gaussian polynomials are polynomials in q of degree NM. Moreover, 
to avoid confusion with the restricted partitions studied earlier, we shall refer 
to the above partitions as Gaussian partitions. As a consequence, we define 
the Gaussian partition operator LGP,k,N,M[-] as 



LGP,k,N,M 



M,Min{[k/2],M},--;Min{[k/i],M} 

E ■ (150) 

ni,n2...,n]v=0 



As a result, we have i^GP,fc,Af,M[l] =Pg{N, M, k). 

We now turn our attention to a more complicated example — the problem 
of determining discrete or distinct partitions. By discrete partitions, we mean 
those partitions in which the elements appear at most once, if at all. Like 
all the preceding examples, they too represent a subset or class of the set 
of integer partitions. As we shall see in the next section, such partitions 
figure prominently in the theory of partitions. Because of this, we define the 
discrete partition operator, L£)pfc[-], by 



1,1 1 

^DP± 



E • (151) 



n\,n2...,nk=0 

Therefore, the only difference between this operator and the partition op- 
erator introduced in Sec. 3 is that the upper limits of the summations are 
restricted to unity, whereas for the latter the upper limits were set to [k/i] 
for each element i. 

A computer program that generates discrete partitions is presented in 
its entirety in the appendix. As stated previously, the number of partitions 
summing to 100 or p(lOO) equals 190 569 272, which is the reason why it is 
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Figure 3: The number of discrete partitions to the total number of partitions 

versus k 



cumbersome to evaluate the 100-th coefficient via the partition method for 
a power series expansion. On the other hand, if we run the program in the 
appendix called dispart to determine the discrete partitions summing to 
100, then we find that it only takes 5 CPU seconds to generate the 444 793 
partitions, which represent 0.2 per cent of p(lOO). Moreover, Fig. [3] displays 
the ratio of the number of discrete partitions or L£)p^k[^] to the number of 
standard partitions versus k for k < 50. From the figure we see that this 
ratio decreases monotonically for k> 10, reaching a value of 0.0179 ■ ■ ■ when 
k = 50. 

The final example we shall consider in this section is the evaluation of 
the conjugate partition as a partition is generated. According to Ref. [TS] . 
the conjugate partition is obtained by transposing the rows and columns of 
the corresponding Ferrers diagram for the original partition. In a Ferrers 
diagram a partition is represented by an array of dots in which the first 
element, say ai, is allocated a row of ai dots, while the next element (02) 
is represented by another row of 02 dots immediately below the first row of 
dots. The process of allocating rows of dots for each element in a partition is 
continued up to the final element in the partition. The conjugate partition or 
a^ of the partition {ai, a2, ■ ■ ■ , a^} is obtained by transposing the rows and 



columns of the Ferrers diagram for the partition. For example, the Ferrers 
diagram for the partition {1, 1,2,3,3,4} has two rows with one dot, a row 
with two dots followed by two rows with 3 dots and finally a row with 4 dots. 
The transpose is obtained by counting the dots in the vertical columns of 
the Ferrers diagram. Hence, a^ for our example is found to be {6,4,3, 1}. 
Note also that the conjugate partition does not necessarily possess the same 
number of elements as the original partition. 

A code that determines the conjugate partition for each partition gener- 
ated by partgen is presented in the appendix. As with all the other examples 
in this section, it is termgen that has been modified. The major complica- 
tion in this code compared with the others is that one is required to create 
and allocate space for a two-dimensional array ferrers of size tot x tot in addi- 
tion to creating a one-dimensional array of pointers called rptr to it. In both 
cases the arrays are of integer type. In reality they should be of type char 
since a true Ferrers diagram consists of dots. That is, instead of allocating 
dots the program allocates ones in creating a Ferrers diagram. Nevertheless, 
by adding the ones vertically, one obtains the conjugate partition. For the 
purists it should be a simple matter to alter the termgen to count dots 
rather than ones. 

When the program called transp for the partitions summing to 4, the 
following output is produced: 
Partition 1 is: 1(4) and its conjugate is: 4(1) 
Partition 2 is: 1(1) 1(3) and its conjugate is: 1(2) 2(1) 
Partition 3 is: 2(1) 1(2) and its conjugate is: 1(3) 1(1) 
Partition 4 is: 4(1) and its conjugate is: 1(4) 
Partition 5 is: 2(2) and its conjugate is: 2(2) 

As can be seen from the output, the conjugate partitions are printed out in 
a different order to those produced by the BRCP algorithm. For example, 
the conjugate partition to Partition 2 or {1,3} is printed out as 1(2) 2(1) 
or {2,1,1}, whereas the form produced by the BRCP algorithm, viz. Par- 
tition 3, is 2(1) 1(2) or {1,1,2}. Throughout this work we have not been 
concerned with the order of the elements in the partitions. When the order 
is important, partitions become compositions ^T\. For example, there are 
three compositions of the partition {1,1,2} are: (112), (121) and (211). We 
shall not be concerned with compositions here. 

There are some interesting properties concerning conjugate partitions. 
For example, there is at least one partition whose conjugate is itself, although 
it can turn out to be a different composition. Such self-conjugate partitions 
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arise when the partitions sum to a square of an integer k since the partition 
k{k) in the multiphcity representation is a self-conjugate. On the other hand, 
if the partitions sum to k{k + l)/2, then the partition {1,2,. . . ,k} is a self- 
conjugate, although again the composition is different. When partitions sum 
to 2k for k > 2, the partition given by {k — 2)(l)l(2)l(fc) is also a self- 
conjugate, while if they sum to 2k+l, where k > 7, then the partition given 
by {k — 6)(l)l(2)l(3)l(/c — 2) is another self-conjugate. In fact, according to 
Ref. |30], the number of self-conjugate partitions is the same as the number 
of partitions with distinct odd elements. 

There are other problems that can be solved by modifying termgen given 
in Sec. 2. One such problem is determining those partitions with only odd 
elements in them. As a consequence, we can also define an odd-element 
partition (OEP) operator as 

2fc+l,[(2fc+l)/3],...,l 

LoEP,2k+l[-] = Yl ' (1^2) 

ni,n3,...,n2k+i=0 

Efci(2i+l)n2i+i=2fe+l 

which is only valid for partitions summing {2k+l), while for those summing 
to 2k, the operator becomes 

2fc,[2A:/3],...,l 

LoEP,2k[-]= Yl ■ (^^3) 

ni,n3,...,n2fe=0 

To obtain partitions with only odd elements in them, all we need to do is 
insert the following for loop at the beginning of termgen just after the type 
declarations: 

for (i=0;i<tot ; i++){ 

freq=part [ i ] ; 

if (i % 2 &fe freq > 0) goto end; 

} 
printf ("%ld: " ,term++); 

It has already been mentioned that the number of partitions with discrete 
odd elements equals the number of self-conjugate partitions. To obtain the 
former, the upper limits in the definitions for LoEP,2k+i[-] and LoEP,2k[-] must 
be set equal to unity. In addition, if we wish to generate partitions with 
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discrete odd elements in them, all we need to do is introduce the above for 
loop at the beginning of termgen in the program dispart that is presented 
in the appendix. With regard to self-conjugate partitions, transp in the 
appendix would need to be modified so that the original partition is stored 
in a temporary array before it undergoes conjugation. Then a test would 
need to be introduced to see if both partitions are identical to one another. 
If they are, then the partition is printed out. Otherwise, it is discarded. This 
problem is left as an exercise for the reader. 

In a similar manner we can define an even-element partition (EEP) oper- 
ator that only applies to those partitions summing to 2k with even elements 
in them. This is defined as 



k,[k/2],...,l 



^EEP,2k\-\ 



E 



(154) 



n2,n4,...,n2fc=0 



The number of partitions generated by this operator is equal to the number 
of partitions summing to k. Hence, we find that 



LEE,2km — -^P,fc[l] 



(155) 



In Ref. [3] the cosecant numbers or Ck are first derived in terms of the parti- 
tions summing to 2k with even elements in them before the form in terms of 
the partition operator given by Eq. ( llOOp is derived. Consequently, we arrive 
at the following identity: 





Lp^k 


[ '^ ^'=-ni,(2. + l)!j n^\ 






= LEEP,2k 




where A^^ -- 


= Eli 


n2i- 







(156) 



In this section various programs have been presented for obtaining specific 
types/classes of partitions, which can be regarded as subsets of the total num- 
ber of partitions summing to a particular value. Consequently, the various 
operators given in this section represent restricted forms of the partition op- 
erator Lp^kl']- In the previous section we showed how the partition method 
for a power series expansion could be developed into a computer program 
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scanning the entire set of partitions. From the material presented in this 
section, it should, therefore, be possible to determine the contributions that 
the specific partitions contribute to the partition method for a power series 
expansion. In particular, in the next section we shall see that the partition 
function p{k) can be obtained via the partition method for a power series 
expansion involving a restricted set of the partitions summing to k, although 
this restricted set is more difficult to evaluate than many of the examples 
considered in this section. We shall not only present the program that can 
generate this restricted set of partitions, but also present the program that 
outputs the partition function in symbolic form so that it can be handled by 
Mathematica. As for the more important programs discussed in this section, 
the new program will appear in the appendix. Finally, it should be noted 
that by possessing the capacity to adapt the partition method for a power 
series expansion to handle specific partitions, we may be able to determine 
which partitions make the largest contribution to the coefficients in Theorem 
1. As a result, accurate approximations to the coefficients can be considered, 
which may avoid the combinatorial explosion that occurs for large orders as 
described at the beginning of the previous section. 

6 Generating Function for P(k) 

An important topic in the theory of partitions are the generating functions 
whose power series expansions possess coefficients that are dependent on the 
properties of partitions. One of the greatest achievements in this context is 
the derivation of the asymptotic formula for the number of partitions p{k). 
The first step that led to this formula was the derivation of a remarkable 
formula for the product P{z) in Equivalence ( 175]1 by Dedekind. As explained 
in Ref. |T8], this can be derived by the application of standard analytic 
techniques, namely Poisson's summation formula, to the logarithm of P{z). 
Then by studying the behaviour of Dedekind's formula for lnP(exp(— t) with 
3fJt > 0, Hardy and Ramanujan [32] were able with amazing insight |31] to 
deduce the asymptotic behaviour of the partition function p{k) for large k. 
Eventually, the asymptotic behaviour of p{k) was completely evaluated by 
Rademacher [33], which culminated in the now famous Hardy- Rademacher- 
Ramanujan formula mentioned in the introduction to Sec. 4. 

Although we shall not reach such lofty heights, we shall, nevertheless, turn 
our attention in the remainder of this work to how the partition method for 



85 



a power series expansion can be applied in the analysis of the various gener- 
ating functions that occur in the theory of partitions and their extensions or 
generalisations. We begin in this section by applying the partition method 
for a power series expansion to P{z), but before we can embark upon this 
task, we need to determine when the power series expansion or generating 
function is convergent or in another words, for what values of z Equivalence 
f l75]) becomes an equation. According to Knuth [TB], it was Euler who noticed 
that the coefficient of 2;" in the infinite product oi {1 + z + z"^ + z^ + ■ ■ ■ + 
z^ + ■■■){! + z'^ + z'^ + ■■■ + z'^'' + ■■■){! + z^ + z^ + ■■■ + z^'' + ■■■)■■■ is the 
number of non-negative solutions to k + 2k + 3k + ■ ■ ■ = n or the partition 
function p{n) and that 1 + z"^ + z^™ + ■ ■ ■ equals 1/(1 — z™"). As a result, he 
arrived at Equivalence (1751) except that the equivalence symbol was replaced 
by an equals sign, which is not entirely correct as can be seen by the following 
theorem. 

Theorem 2. The equivalence statement relating P{z) or nm=i ^/i^ " 
z^) to the generating function given by the power series expansion with 
coefficients equal to p{k), viz. Equivalence (1751) . is absolutely convergent for 
1^1 < 1, in which case the equivalence symbol can be replaced by an equals 
sign. On the other hand, it is divergent for all other values of z. Then P{z) 
represents the regularised value of the series on the rhs. For \z\ = 1, the 
generating function is singular. 

Proof. The reason why an equivalence symbol has had to be introduced 
is due to Euler's second observation concerned with the geometric series. 
Replacing the series by its limit value of 1/(1 — z'"^) is strictly not valid for all 
values of z as described in Refs. [2], [7j, [8] and [11]. There it is shown that 
the standard geometric series, i.e. YlT=o ^^^ ^^ divergent for 5R z> 1, absolutely 
convergent for |2;| < 1 and conditionally convergent for |-2| > 1 and ^z <1. 
Regardless of the type of convergence, the limit value of the series is found to 
equal 1/(1 — z). For ^z>l, however, summing the series yields an infinity. 
In this case, if the infinity is removed after the series is summed, which is 
the essence of regularisation, then the remaining finite part is found to equal 
1/(1 — z) again. Hence, for z > 1, the regularised value of the geometric 
series is equal to 1/(1 — z). Moreover, along the line ^z = 1, the limit of 
the series is undefined or indeterminate, while at the point z = l, where the 
line is tangent to the unit disk of absolute convergence, it is singular. This 
type of behaviour is expected since ^z = l, represents the border between 
the convergence to the left and divergence to the right. Since the limit value 
is the same on both sides after regularisation, we set the regularised value to 
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l/{l-z) along 3f? 2 = 1. 

With regard to Equivalence (1731) we are dealing with an infinite product 
of geometric series involving different powers of z in the limit value. Never- 
theless, we can use the above knowledge of the geometric series to determine 
where the series on the rhs of Equivalence fl73|) is convergent and where it 
is divergent. For the product on the Ihs to equal the series on the rhs of 
Equivalence (175]) . we must have for all positive integer values of /, 3fJz' < 1. 
For / = 1 we end up with the standard geometric series, but for 1 = 2, the 
series will now only be convergent for ^ z"^ < 1 or —1 < ^ z < 1. Thus, the 
range of values of z has changed, which means that the convergence of the 
series of the series on the rhs of Equivalence (175]) will be affected by each 
value of / or each series in the product. 

Let us examine the third series in the product, whose limit is 1/(1 — 2;'^). 
In order to analyse this version of the geometric series, we write the limit 
value as 

(157) 



Z' 



3 (l-;2)(l-2e2W3)(l_^-2W3) 



Decomposing the rhs into partial fractions, we see that this version of geomet- 
ric series is actually the sum of three geometric series, each with a different 
limit. The first yields the standard geometric series discussed above. The 
second series has a limit of 1/(1 — zexp(227r/3)). In this case we replace z 
by zexp{2i7T/3) and continue with the same analysis. Then the second series 
is convergent for 3fJ (2;exp(2z7r/3)) < 1 or y < (2 — x)/\/3 when z = x + iy, 
while it is divergent for y > (2 — x)/^/3. That is, the line 'Siz = 1 sepa- 
rating the regions or planes of convergence and divergence has been rotated 
by 27r/3 in a clockwise direction. The "left side" of the line representing 
where the series is convergent is now given by y < (2 — x)/\/3. On the other 
hand, it is divergent for y > [2 — x)/a/3 in which case the limit becomes 
the regularised value of the second series. The third series, whose limit is 
1/(1 — 2exp(— 2i7r/3), represents the opposite of the previous series. That 
is, the singularity at 2; = 1 in the standard geometric series has now been 
rotated by 27r/3 in an anti-clockwise direction. Hence, the third series is 
convergent for 3fJ (zexp(— 2i7r/3)) < 1 or y > —{x + 2)/-\/3 when z = x + iy. 
This represents the "left" side, while the "right" side or where it is divergent 
is given by y < — (s + 2)/v^. For these values of z the limit represents the 
regularised value of the series. 

It is the intersection of the "left" sides for the three series that represents 
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the region of the complex plane for which the third series in the product 
or YlT=o^^^' ^^ convergent. Outside this region the series is divergent. The 
intersection of the tangent lines yields an equilateral triangle, where the mid- 
points of the edges coincide with the three singular points of the component 
geometric series on the circle \z\ = 1. Moreover, the unit disk of absolute 
convergence is circumscribed by this triangle. Those parts of the triangle not 
in the unit disk of absolute convergence represent the regions of the complex 
plane for which the third series in the product is conditionally convergent. In 
total they are significantly less than the corresponding region of the complex 
plane for the second series in the product, which we found is given by the 
region outside the unit disk of absolute convergence in the plane —l<^z< 
1. 

If we consider the fourth series in the product, i.e. X^fclo '^^'^' then decom- 
posing its limit into partial fractions yields four distinct geometric series with 
the singularities situated again on the unit circle, but at ±1 and exp(±i7r/2). 
If we draw tangent lines through each of these singularities, then we find the 
common region or the intersection of their "left" sides is now a square cir- 
cumscribing the unit disk of absolute convergence. The singularities in the 
component geometric series appear at the mid-points of the square's sides. 
Outside the square, the fourth series in the product is divergent. Then the 
regularised value of the series is obtained by combining the limits of the 
component series. The regions inside the square, but outside the unit disk, 
represent the values of z for which the fourth series in the product is condi- 
tionally convergent. As expected, these regions in total are less than either 
of the regions of conditional convergence for the second and third series in 
the product. 

If we continue this analysis to the /-th series in the product, i.e. for 
^fcLo^"^' then we find that the intersection of tangent lines yields an l- 
sided polygon that circumscribes the unit of disk of absolute convergence. 
For the values of z outside the polygon the series will be divergent, while 
for those values of z within the polygon, but outside of the unit disk, the 
series will be conditionally convergent. Furthermore, as higher values of / are 
considered, the number of tangent lines not only increases, but also the total 
region of conditional convergence contracts. In the limit as Z — )■ cxd we will 
be left with the unit disk as the sole region where the series is (absolutely) 
convergent, while outside the disk the series is divergent. Since the product 
in Equivalence fl73|) includes all values of /, the / = oo limit by virtue of the 
fact that it possesses the smallest region of convergence in the complex plane 



determines the values of z, where the series on the rhs of Equivalence (I73p 
is convergent. This means that the series or generating function, which we 
shall call from here on the partition number series, is equal to the product 
on the Ihs only for z situated in the unit disk. Outside the unit disk, the Ihs 
represents the regularised value of the series and an equivalence symbol must 
be used instead of an equals sign. Finally, the circle |2;| = 1 represents a ring 
of singularity separating the divergent values of the partition number series 
from the absolutely convergent values. In the case of the standard geometric 
series, there is only one point where the absolutely convergent region is sep- 
arated from the divergent region, namely the singularity situated at z = 1. 
Elsewhere, the line ^z = 1 separates conditionally convergent values from 
divergent values. As mentioned above, the limit for the geometric series is 
indeterminate along the line, but is assigned a regularised value of 1/(1 — z). 
However, for z = l the regularised value yields infinity. Hence, it can be seen 
that there is a difference between separating absolutely convergent values 
from divergent values and separating conditionally convergent values from 
divergent values. This completes the proof of the theorem. 

For |2;| <1, we can invert Equivalence fl75]) . thereby obtaining 

1 °° 1 

Since the above product produces a power series that is vahd for all values 
of z, we can write it as 

oo oo 

{z; zU = l[{l~z')=l + Y, q{k) z^ . (159) 

fc=l k=l 

The leftmost expression is a special case of the q-Pochhammer symbol [55] . 
which is defined as 

ra-l 

{a-z)n = \[{l-az^) . (160) 

fc=0 

According to Knuth [18] , it was Euler, who was the first to discover that 
much cancellation occurs when multiplying the various terms in the infinite 
product given in Eq. (I160p . Specifically, he found that 

oo oo 

J](l-z'") = l + 5^(-l)^(z(3^''-'^)/2 + z(=^'='+'=)/2) . (161) 

m=l fc=l 
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Therefore, comparing the above result with rhs of Eq. fll59p we that the q{k) 
are frequently equal to zero and when non-zero are either equal to 1 or -1. 
The values of k for which the q{k) vanish are known today as the pentagonal 
numbers [36] • They are themselves a particular case of a broader class of 
numbers known as the figurate or figural numbers [37| [38] . By applying 
Theorem 1 to Eq. fll58p . we see that the coefficients of the outer series, viz. 
Qk, are equal to (—1)'^, while those for the inner series or the pk are equal to 
the partition function p{k) for k > 1 and zero for ^ = 0. Therefore, with the 
aid of Eq. (fTOj) we arrive at 



q{k) 



^P,k 



-l)^'=iV,!j] 



p{i) 



i=l 



rijl 



(162) 



As a consequence of the fact that the q{k) are non-zero when k= (3j^ ± j)/2, 
Eq. (1162P can also be written as 



^P,k 



■i^'^iVfcin 



p{i)"' 



i=l 



riA 




k = {3f±j)/2, 
otherwise . 



(163) 



Although they do not give the actual number of discrete partitions, we 
shall refer to the q{k) as the discrete partition numbers. Shortly, we shall 
see how these coefficients are related to the number of discrete partitions. 
They do, however, have an interesting connection with the partition function, 
which follows when both power series on the rhs's of Eqs. (I158P and (11590 
are multiplied by one another. Then we find that 



$^2;'=$^p(j)g(fc 
i=o 



(164) 



fc=0 



where p{0) = q{0) = 1. Again, since z is fairly arbitrary, like powers of z can 
be equated on both sides of the above equation. For A; > 1, we obtain the 
following recurrence relation: 



^p{j)q{k-j) = 



(165) 



i=o 



This is simply Euler's recurrence relation for the partition function, which is 
given by Eq. (20) in Ref. P^. Occasionally, it is referred to as MacMahon's 
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recurrence relation as in Eq. (20) of Ref . [32] • Because most of the discrete 
partition numbers or q{k) vanish, it means that only a few of the previous 
values of the partition function are required to evaluate the latest value of 
the partition function. 

Eq. (11621) is an interesting result, but perhaps, not very practical for deter- 
mining the discrete partition numbers, when it is realised that the partition 
function or p{k) grows exponentially. We can, however, use the partition 
method for a power series expansion to derive another result for the discrete 
partition numbers. First, assuming that |-z|<l, we write 1/P{z) as 



oo oo 



-p|3 = exp(5^1n(l-.'"))=exp(-5^5^(.-)7j) • (166) 

^ ^ m=l m=l j=l 

The Taylor/Maclaurin series for the logarithm has been introduced in this 
result, since it is absolutely convergent for |;z| < 1 [I]. The double sum can 
be expressed as a single sum by realising that the coefficients of z can be 
expressed as a sum over divisors or factors of the power [31]. Then we arrive 
at 

fc=i 

where jj = Xldii ^/^ ^^'-^ ^ represents a divisor of j. That is, the sum is only 
over the divisors of j. Some values of the ■jj are: 71 = 1, 72 = 3/2, 73 = 4/3, 
74 = 7/4, and 75 = 6/5. More explicitly, we find that 76 = 1/6 + 1/3+1/2 + 1 = 
2, while for the case of j = p"^, where p is a prime number, the sum yields 

(1 - 1/P'"+^) ....^ 

%- = (1 _ 1/^) • (168) 

Eq. (11681) can be derived simply by using the limit for the geometric series. 
Furthermore, Eq. (I167p also means that 

^ 1 / 3^2 4^3 7z^ , \^ . . 

P(z) = l + J2-(^z + — + — + — + . .■ + ^,z^ + ...j . (169) 

fc=i 

Therefore, we have obtained alternative representations for the generating 
functions of both P{z) and its inverse. 
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Now we are in a position to apply Theorem 1 to Eq. fll67p . whereupon 
we see that the coefficients of the inner series pk equal 7^ for k >1, while 
Po = 0. On the other hand, the coefficients of the outer series qk are equal to 
(— l)'^//c!. Hence, from Eq. f lTO|) the discrete partition numbers are given by 



qik) 



^P,k 



\Ni 



n 



ni^ 



(170) 



Hence, we have derived an alternative version of Eq. fll62p . Moreover, when- 
ever /c = (3j^ ± j)/2, for i an integer, this result can also be written as 



-P,(3j2±i)/2 



4)^(3.^..)/. J] :h_ 

1=1 



(171) 



For all other values of fc, the sum over all partitions in Eq. (I170p vanishes. 
Therefore, to evaluate ^'(6), we require all the 7^ ranging from j = 1 to 6, 
which have already been given above. Summing over the eleven partitions 
summing to 6 in Eq. (I170p yields 



g(6) = -2 + 6/5 + 21/8 - 7/8 + 16/18 - 2 + 4/18 
- 27/48 + 9/16 - 3/48 + 1/6! = , 



(172) 



which is indeed the value of this discrete partition number. By using these 
results, the reader can readily verify that g(0) = 1, g(l) = —1, g(2) = —1, 
g(5) = l and g(3) = g(4) = 0. Moreover, since Dq in Theorem 1 is non-zero in 
this case due to the fact that g(0) = 1, we can use Eq. flT2|) to determine the 
coefficients of the inverted power series expansion or the generating function 
for P{k). Hence, we find that 



p{k) 



^P,k 



\N 



iv-n 



q{i) 



i=l 



riA 



(173) 



This result, which represents the inverse of Eq. (11620 . incorporates much re- 
dundancy since the q{i) are only non-zero when for specific values of i. Con- 
sequently, both the sum over the partitions and the product are only non-zero 
for those values of i, which are of the form of (3j^— j)/2 or (3j^+j)/2, where 
j is an integer ranging from 1 to jm = [(1 + vl + 24A;)/6]. In addition, when 
the q{i) are non-zero, they are only equal to unity in magnitude. Therefore, 
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it is the factor of A^! that is responsible for the exponential increase in the 
partition function as k increases, although this factor will often be countered 
by the l/rij terms in the denominator of the product. For example, if we 
wish to determine p(6), there will not be any sums over n^, 77.4 and uq in the 
above equation since g(3), g(4) and g(6) vanish, while in the product, 71,3!, 
nj and tiqI will equal 0! or unity. Furthermore, since g(l) = g(2) = —1 and 
g(5) = 1, Eq. (^W> yields 



p(6) = 1 



+ (-1)^(-1)^(-1) 5!/4! + (-l)^(-l)2(-l)24!/(2! ■ 2!) 



-l)^3!/3! + (-l)2(-l)2! = ll 



(174) 



Hence, the contributions from the five relevant partitions are positive except 
for the last one, which represents the contribution due to the partition {1,5}. 
On the other hand, if we apply Theorem 1 to Eq. f ll69p . then the only 
difference to the previous evaluation of the discrete partition numbers is that 
the coefficients of the outer series g^ are now equal to l/k\. That is, the 
coefficients of the outer series are still equal to 7^. Hence, we arrive at 



p{k) 



^P,k 



n 



n-i -, 



Lp_, \\^ 



1 



(175) 



Therefore, we have an entirely different means of evaluating the partition 
function with the sum of the reciprocals of the divisors of each element being 
the assigned values in the partition method for a power series expansion. 
Consequently, Eq. (I165p becomes 



i=o 



Eip.-n^i''^'(-i)"-n^ 






n,;! J 



;i76) 



In deriving these new results for the partition function p{k) we have en- 
countered a more complex example involving a specific class of partitions 
than any of those studied in the previous section. In this case we only re- 
quire partitions whose elements are pentagonal numbers or of the form of 
(3j ± l)j/2, where j is any integer lying between zero and jm- In view of the 
importance of the Eq. (I173p . let us consider introducing the modifications to 
the program partgen presented in Sec. 2. As in the examples of the previous 
section, most of these modifications will be made to the function prototype 
termgen. 
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The first modification is that we need to make is introduce the math 
hbrary with the other header files. This is needed so that the floor function 
can be called to evaluate the maximum value of j or im when the value of 
k or tot is typed in by the user. Since jm is called in both the main and 
termgen prototypes, it must be declared as a global variable. Once these 
modifications are carried out, we can concentrate on the changes that are 
required for termgen, which is displayed below: 

void termgen () 

{ 

int freq ,i ,j , jval ; 

for (i=0;i<tot ; i++){ 
j val =0; 
freq=part [ i ] ; 
if (freq >0){ 

for (j=l;j<=j_m; j++){ 

if (i = ((3*j-l)*j-2)/2) jval=j; 
if (i = ((3*j+l)*j-2)/2) jval=j; 

} 
if ((jval==0) &fe (freq>0)) goto end; 

} 
} 

printf("%ld: " ,term ++); 
for (i=0;i<tot ; i++){ 

freq=part [ i ] ; 

if (freq) printf ("%i (%i ) ",freq,i+l); 

} 

printf("\n"); 

end : ; 

} 

Basically, this version of termgen is similar to those that appeared in 
the previous section. That is, before a partition is printed out, testing is 
done in the function prototype to see if each partition belongs or conforms 
to the particular set of partitions, which the user desires. In this case our 
aim is to print out only those partitions whose elements can be written in 
the pentagonal number forms of (3j^ — j)/2 and (3j^ + j)/2, where j is an 
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integer. This is accomplished by introducing another for loop in the func- 
tion prototype, which evaluates the value of j called jval, for the appropriate 
elements. If an element is not of the required form, then jval remains zero. 
Otherwise, it is non-zero. If jval is zero for an element, then the partition 
is examined to see if there any occurrences of the element by checking the 
variable freq. If it is greater than zero, then the entire partition is discarded 
by the goto statement. This test is carried out on all elements in the parti- 
tion. The procedure is then applied to all partitions summing to tot. Only 
those partitions whose elements are of the required form are printed out by 
the code. As an example, when k = 6, the output for this program called 
partfn is: 
1: 1(1) 1(5) 
2: 4(1) 1(2) 
3: 6(1) 
4: 2(1) 2(2) 
5: 3(2) 

The above partitions represent the five that contributed to the calculation 
of p{6) given above. Interestingly, when the code is run for partitions sum- 
ming to 100, it only prints out 42205 partitions which represents about 0.02 
percent of the total number of partitions or p(lOO). Moroever, the num- 
ber of partitions summing to k whose elements are pentagonal numbers or 
Lpentei.kl^] starts off greater than the number of discrete partitions for the 
same value of k, but then drops away when k > 17. Fig. H] displays the ratio 
of Lpente«,fc[l] to the numbcr of partitions summing to k for k < 50. Once 
again, it is monotonically decreasing for k beyond a certain value, which in 
this case is about 8. 

It was also stated at the end of the previous section that the partition 
method for a power series expansion can be adapted to handle situations 
where only a subset of the partitions summing to a particular value are re- 
quired as in the above example. We shall demonstrate this by modifying 
the code in Sec. 4 so that the partition function p{k) can be evaluated via 
Eq. 0173p . The resulting code called pfn is presented in its entirety in the 
appendix. Basically, it expresses the partition function in a symbolic form 
where the final values can be evaluated by introducing the output into Math- 
ematica. As expected, termgen has the same for loop presented above. The 
interesting feature about the code is that the for loop appears twice at the 
beginning of termgen. This is necessary because the first step prints out 
p{k) := before considering the single element partition of {k} in the first 
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Figure 4: The ratio of the number of partitions whose elements are pentag- 
onal numbers to the total number of partitions versus k 

branch of the if statement. Clearly, the single element partition may not 
be of the required form — hence, the first appearance of the for loop. All 
the other partitions are processed by the second branch of the if statement, 
which means that we require the for loop in this part of the program in order 
to determine which are of the required form. 

The code is also capable of expressing the final symbolic form for the 
partition function in two different forms. In the first form a printf statement 
prints the final form for p{k) in terms of the q{i) as they appear in Eq. fll73p . 
This statement has been commented out in favour of the version of the code 
appearing in the appendix. In the second version of the code the value of 
each q{i) is evaluated. That is, the code prints out -1 to the power of jval 
for the elements in a suitable or candidate partition. In order to accomplish 
this, the for loop mentioned in the previous paragraph has had to be intro- 
duced into the latter part of termgen again. Finally, the quantity q[N]a^ 
appearing in the output of the code discussed in Sec. 4 has been replaced 
(— l)^(A^) in pfn. Therefore, if the version of pfn in the appendix is run for 
tot = 6, then the following output is produced: 

p[6]:= ((-l)-(l)) ((-1)^(2)) (-1)A(2) 2! + ((-1)^(1))^(4) ((-1)^(1)) (-1)A(5) 
5!/4! + ((-1)^(1))^(6) (-1)^(6) + ((-1)^(1))^(2) ((-1)^(1))^(2) (-1)^(4) 4!/(2! 
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2!) + ((-1)^(1))^(3) (-1)^(3) 

Time taken to compute the coefficient is 0.000000 seconds. 
This output can be imported into Mathematica whereupon it gives the cor- 
rect value of 11 for p{6). When the code is run for tot equal to 100, it takes 34 
seconds to compute the symbolic form for p(lOO) or Lp_ioo[l] and only 0.27 
seconds to produce the value of 190 569 292 in Mathematica on the same 
Sony VAIO laptop mentioned in previous sections. However, whilst this is 
not too bad, the fastest method of obtaining p{k) or -^p,fc[l] from the various 
programs considered in this work is to comment out termgen and introduce 
the statement 

term++; 

immediately below in partgen of Sec. 2. Then term needs to be initialised 
to zero and a printf statement introduced into main. When the resulting 
program is run for tot equal to 100, it only takes 3 seconds to compute 
p(lOO). Nevertheless, both methods for computing the partition function 
will slow down dramatically as tot continues to increase due to combinatorial 
explosion. This is where either Eq. (I165p or even the Hardy- Ramanujan- 
Rademacher formula should be used. In fact, we may write 



lim Lp^k 

fe— >oo 



^ ^^ exp(VW37r) , (177) 



and 

-Ifmfl'-^ -, -^ exp(v/2i73.) . (178) 



hm Lpfc 



fc— >oo 



iWi-, 



. 1 rii 



AVSk 



In addition, the arrow symbols in the above results can be replaced by equals 
signs. 



7 Generalisation of the Inverse of P(z 



In the previous section we demonstrated how the partition method for a 
power series expansion can be used to derive different forms for the generating 
functions of the product P{z) and its inverse or reciprocal. In the case of 
P{z) two different forms for the coefficients of the generating function or the 
partition function p{k) were obtained in terms of the partition operator acting 
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on different arguments. The first given by Eq. fll73p involved the discrete 
partition numbers q{k), which represent the coefficients of the generating 
function for the inverse of P{z), while the second given by Eq. fll75p involved 
a sum over the inverses of the divisors for each element i in a partition. For 
the inverse of P{z) two different forms for the discrete partition numbers were 
also obtained, but now the partition operator was found to act on either the 
partition function as in the first form given by Eq. fll62p or in the second form 
given by Eq. fll70p the same sum over the inverses of the divisors with an extra 
phase factor of (—1)^*, where Nk represents the number of elements in each 
partition summing to k. This extra phase factor is responsible for ensuring 
that the discrete partition numbers equal either ±1 when fc is a pentagonal 
number or zero, otherwise. Because the phase factor does not appear in Eqs. 
f ll73p and (11750 . the partition function experiences exponential growth as 
can be seen from the Hardy-Ramanujan-Rademacher formula given above. 
Nevertheless, it was possible to derive all these results because the coefficients 
of the powers of z in the product P{z) are simple, namely equal to -1. In 
this section we aim to generalise P{z) to the situation where the coefficients 
of z'' are now equal to C^. We begin by studying the inverse of P{z). As a 
consequence, we arrive at the following theorem. 
Theorem 3. The infinite product defined by 

oo 

H{z) = l[{l + az') , (179) 

can be written as a power series expansion or generating function of the form: 

oo 

H{z) = 1 + J2hkz'' , (180) 



fc=i 



where the coefficients hk can be expressed in terms of the discrete partition 
operator defined by Eq. (I15ip and are given by 



k 

hk = Lop,k 

1=1 



K 

[n^"] • (181) 



Proof. If we multiply out the lowest order terms in z in the product 
given in Eq. (I179p . then we obtain a power series expansion for -/^(-z) in a 
similar manner to the proof of Theorem 1. The zeroth order term in the 
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resulting power series is unity, while the first order term becomes CiZ. We 
then find that the second order term in the power series becomes C2-2^. Once 
we go beyond second order, however, the coefficients become more complex 
to evaluate due to an ever-increasing number of terms appearing in them. 
For example, the third, fourth and fifth order coefficients are respectively 
equal to C3 + C1C2, C4 + C1C3, and C5 + C1C4 + C2C3. In fact, on close 
inspection we find that the coefficient of z^ in the power series depends on 
the number of discrete partitions summing to k. For example, the third and 
fourth order coefficients are composed of only two terms because there are 
only two discrete partitions for k = 3 and k = 4, while the fifth order term 
is composed of three terms due to the three discrete partitions summing to 
5, viz. {5}, {1,4} and {2,3}. Therefore, instead of summing over all the 
partitions summing to k as we did in Theorem 1, we need only sum over 
the subset comprising the discrete partitions summing to k, which we have 
already seen is significantly less than p{k). Hence, the sum over all partitions 
in Eq. f jTOj) simplifies drastically with all the frequencies lying between zero 
and unity, not between zero and [k/i] for each element i. Furthermore, since 
there is no outer series in the expansion of the product, we can drop q^ and 
a^ in Eq. ( ITOj) . In addition, there is no multinomial factor. Consequently, 
we find that according to Eq. (fTOj) the coefficients hk in the power series for 
H{z) are given by 

1,1,1,...,! k 

hk= Y. HCT • (182) 

ni,n2,n3,...,nj.=0 «=1 

The sum over partitions in the above result is simply the discrete partition 
operator as defined by Eq. fllSip . When this is introduced into Eq. f ll82p . we 
arrive at Eq. fllSip . This completes the proof of the theorem. 

It should also be mentioned that we can invert the rhs of Eq. f llSOp and 
apply Theorem 1 again. Then we obtain a power series expansion for 1/H{z). 
As a result, we find that 



Hk 



^P,k 



W^ 






(183) 



while from Eq. f ll2p . we obtain 



hk 



WP± 



He. 



N, 



^P,k 



■i)'^'=iVfc!n 



h: 



Hi 



rii 



(184) 



99 



Since Theorem 3 is quite general, it means conversely that any power 
series expansion can be expressed as an infinite product of the form given 
by the rhs of Eq. fll79p . For example, as discussed on p. Ill of Ref. [18] the 
geometric series can be represented by the following product: 

oo oo 

fc=0 J=0 

In this example C2» = 1, while the other Cj simply vanish. This means 
that the product of ni=i ^T i^ ^^e hk is either zero or unity. In addi- 
tion, since the coefficients of the geometric series are equal to unity, we have 
Ldp^A: [ni=i Q"'] = 1- Since all the terms in the discrete partition operator are 
either equal to zero and unity, this means that there can only be one discrete 
partition where all the elements are of the form 2^ with each value of j lying 
between zero and [log k/ log 2] . For fc = 5 and fc = 7, these discrete partitions 
are respectively {1,4} and {1,2,4}, while for A; = 11 it is {1,2,8}. 

Another example of a well-known power series that can be expressed as 
an infinite product is the exponential function, which can be expressed as 

e^ = (1 + y){l + yV2)(l - |/V3)(1 + 3yV8)(l - 1/V5)(1 + W H'^) 
X (1 -y77) (1 + 27yVl28)(l-8yV81)(l + 912/^7800) ■■■ . (186) 

From Eq. f ll86p we find that Ci = 1, while the other Ci can be determined by 
the following recurrence relation involving the divisors of i: 

(_-\\d+l 

d\i 

This means that whenever i is a prime number greater than 2, say p, we find 
that Cp = —l/p- Eq. f ll87p is left as an exercise for the reader. 

We can derive other interesting results by taking the logarithm of both 
Eqs. flT79D and fITSOD . which yields 

oo oo 

In^l + ^/ifc^^^) =5^1n(l + Ci2') . (188) 

fc=l i=l 

By introducing the Taylor/Maclaurin series expansion for logarithm into the 
rhs, we obtain 



oo oo oo 



ln{l + J2hz')=J2T.-^C'''' ' (189) 



fc=l i=0 i=l 

100 



3 



where we have now assumed that 1^1 < 1. We now apply Theorem 1 to the 
left hand side of the above result. The inner series coefficients are given by 
Pk = hk, while the outer series coefficients are given by qk = {—l)^~^^/k. By 
equating like powers of z in the resulting expression, we arrive at 



^P,k 



■i)"Hiv.-i)!n^l=E^ 

i=l i\k 



^k/i ■ 



(190) 



For the particular case of the geometric series, viz. Eq. f llSSp . hk = 1 and 
C2J = 1, while the other values of the Ck vanish. Introducing these results 
into the above equation produces 



^P,2J 






(191) 



If we put the Ci equal to unity in Eq. (I179P , then we find that 



Y[{l + z'')=l + J2LnP,k 



(192) 



fc=i 



fc=i 



In another words, the coefficients of the power series expansion for the prod- 
uct are given by the number of discrete partitions summing to k. Therefore, 
from Eq. (I18ip we have 



hi 



WP,k 



On the other hand, if we put the Ci = —1, then we obtain 



qik) 



^DP,k 



■1) 



Nk 



(193) 



(194) 



Prior to Eq. (I162p it was mentioned that the q{k) or discrete partition num- 
bers are equal to (— l)-* whenever fc is a pentagonal number [30] or is equal 
to j(3j ± l)/2. Hence, Eq. (I194p can be expressed as 



^DP,k 



-1 



,Nk 



i-iy, forfc = j(3j±l)/2, 
0, otherwise . 



(195) 



The above result tells us that the number of discrete partitions with an odd 
number of elements is equal to the number of discrete partitions with an 
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even number of elements when k is not a pentagonal number. It also tells us 
that the discrete partition numbers q{i) come in pairs of either 1 or -1, the 
former corresponding to a pentagonal number derived from an even number, 
i.e. by setting j equal an even number and the latter to a pentagonal number 
derived by an odd number. A value of 1 in the above result or for j an even 
number, means that the number of discrete partitions with an even number 
of elements is one more than the number of discrete partitions with an odd 
number of elements. A value of -1 or j odd, represents the opposite situation. 
A proof of this result based on Ferrers diagrams appears in Ch. 1 of Ref. [3T] . 
Now if we put Ci = —uj in Eq. (11801) , then the power series for the ensuing 
product becomes 



fc=i 



uoz 



••) 



-U]Z), 



(l-c) 



l + ^g(A;,-a;)z^ 



(196) 



fc=i 



where (— w; -2)00 is again the q-Pochhammer symbol presented in Eq. (11601) 
and 



q{k, -u) = Ldp,i 



-UJ 



lAffc 



(197) 



Thus we see that q{k, —l) = q{k) or the discrete partition numbers. Note also 
that the power of u gives the number of elements in each discrete partition. 
By summing over all partitions each power of u, say w", in q{k,u) gives 
the total number of discrete partitions with n elements in them. A similar 
situation will arise when we study the inverse of Eq. (I196p . Moreover, by 
putting Ci = uj in Eq. (I180p . we obtain 



Q{z,u) = \[{l + ujz^) 
k=i 






^q{k,uj)z^ 



k=l 



where 



q{k,uj) 



^DP,k 



UJ 



Nk 



(198) 



(199) 



The above polynomials will be referred to as the discrete partition polyno- 
mials. From Eq. (11921) we have already seen that the a; = 1 case gives the 
coefficients that represent the total number of discrete partitions summing 
to k, viz. L£ipfc[l], while from the above the discrete partition numbers are 
given by q{k) = Ldp^^K—I)^'']. Hence, we see that the difference between 
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k 


q{k,u) 


p{k,uj) 





1 


1 


1 


u 


u 


2 


u 


u'^ + u 


3 


Uj'^ + UJ 


00^ + 00"^ + 00 


4 


Uj'^ + UJ 


co^ + co^ + 2uP' + uo 


5 


2b? ^bj 


w^ + w^ + 2^3 + 2^2 ^ w 


6 


u^ + 2^2 + w 


w^ + w^ + 2w^ + 3^3 + 3^2 ^ ^ 


7 


u^ + 3^2 + w 


w^ + w*^ + 2^5 + 3w^ + 4^3 + 3^2 ^ ^ 


8 


2u^ + 3w^ + w 


w^ + w^ + 2^6 + 3w^ + 5u;^ + 5u;3 + 4a;2 ^ cj 


9 


3u^ + 4^2 + w 


w^ + w^ + 2w^ + 3^6 + 5u;^ + 6w^ + lu^ + 4^^ 


10 


o;"^ + 4w^ + Alo^ + CO 


^10 + c^9 + 2a;S + 3a;^ + Sw^ + 7u;^ + 9a;^ + 80;^ 
+5tj^ + bo 



Table 3: Discrete partition polynomials q{k,u) and partition function poly- 
nomials p{k,u). 

the discrete partition numbers and the number of discrete partitions for a 
particular value of k is that for the latter the discrete partition operator acts 
on unity when scanning the discrete partitions summing to k, while for the 
former it acts on (—1)^*, where Nk represents the number of elements in 
each discrete partition. That is, the difference is caused again by a phase 
difference in the number of elements in the discrete partitions summing to k. 
The discrete partition polynomials up to A; = 10 are displayed in the 
second column of Table |3l As expected, for u = —1 these results reduce to 
the discrete partition numbers q{k), while for u = l they yield the number of 
discrete partitions. From the table it can be seen that the discrete partition 
polynomials are polynomials of degree n, where n{n + l)/2 < fc < (n + 1) (n + 
2)/2 since the partition with the most discrete elements when fc = ra(ra+l)/2 is 
{1,2,3,. . . ,n}. The lowest order term in u corresponds to the only one-element 
partition, viz. {k}. If we run the program dispart, which is displayed in 
the appendix, then we find that there are 4 discrete partitions summing to 
6, which are {6}, {1,5}, {2,4} and the self-conjugate {1,2,3}. By referring 
to Fig. [H we see that the first partition is just one branch from the seed 
number, the second and third partitions are two branches away and the third 
is three branches away. Hence, we arrive at q{6,uj) = u + 2cj^ + 00^, where the 
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magnitude of the coefficients of u'^ represent the number of distinct partitions 
with k elements. That is, if the number of distinct elements summing to k 
with i elements in them are given by qi{k), then the polynomials can be 
expressed as 

n 

q{k,u) = Y,(l^ikp'' , (200) 

where n(n + l)/2< A; < (n + l)(n + 2)/2. Hence, if we let g(A;, 1) =Ldp^A;[i] or 
the number of distinct partitions summing to k. then we obtain the trivial 
equation of q{k, 1) = X^ILi liW- Moreover, the lower bound on k gives us a 
limit as to the maximum number of elements that can appear in a discrete 
partition, which is given by 



imax = [{V8kTT-l)/2] . (201) 

By fixing the number of elements to i in the discrete partition operator 
so that it becomes 



1,1 

WP,k,i 



E ' (202) 



we arrive at L^ip^ j[l] = qi{k). This means that we need to input a second 
value into the program, which represents the number of elements the user 
desires. If this is set equal to a global variable called numparts, then the only 
changes to be made to dispart are: (1) introduce a local variable sumpart. 
which adds all the values of freq in the ffist for loop and (2) insert the 
following if statement before anything is printed out: 

if(sumparts != numparts) goto end; 

When these modifications are implemented and the resulting code run for 
the discrete partitions summing to 100 with the number of elements, i.e. 
numparts, set equal to 5, one finds that there are 25,337 discrete parti- 
tions beginning with {1,2,3,4,90} and ending with {18,19,20,21,22}. That 
is, 55(100) or -^£)pioo,5[l] is equal to 25 337. According to Eq. (I20ip . the 
maximum number of elements in the discrete partitions summing to 100 is 
13. When the code is run for numparts set equal to 13, 30(= gi3(100)) 
partitions are printed out beginning with {1,2, 3,4,. ..,12, 22} and ending with 
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{1,2, 3,4, 6, 7,. ..,13, 14}. Running the code for higher values of numparts with 
/;; = 100 does not result in any partitions being printed out. Hence, ^^(lOO) =0 
for z>13, which is consistent with Eq. f l20ip . 
From Eq. f ll96p we have 

n (1 - ^'z'') = ^-^^^ = 1 + E ^(^' -^')-'' • (203) 

k=i ^ ' k=i 

The product on the Ihs of this result can also be written as 

oo oo 

\[{l-uj\^'')=\[{l-uz''){l + ujz'') . (204) 

fc=l fc=l 

Introducing the rhs of Eq. (11961) into the rhs of the above equation yields 

oo oo fc 

n (1 - ^'z^'^) =J2''Y1 ^(•^■' ^) ^(^ - ^■' -^) • (205) 

fc=l fc=l j=0 

The equals sign is only valid in Eqs. (I203p and (12051) for |ci;| < 1 and |z| < 1. 
By equating like powers of z on the rhs's of both these equations, we find 
that for k, an odd number equal to 2?7i+l, 

2m+l 

Y, qU, w) g(2m + 1 - J, -a;) = . (206) 

i=o 

On the other hand, for k = 2m, we obtain 

2m 

Yc[iJ,^)c[i'^m-j,-u) = q{m,-u'^) . (207) 

j=0 

When a; = l, Eq. (1207]) reduces to 

2m 

J2 LdPJ [1] LDP,2m-j [(-1)^^--] = LDP,m [("l)^'"] • (208) 

3=0 

From Eq. (I195p we see that the Ihs of the above equation is effectively a sum 
over the pentagonal numbers less than 2m, while the rhs is non-zero if and 
only if m is a pentagonal number. Furthermore, if we let u equal the complex 
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number i, then Eq. (120 7p gives the number of discrete partitions summing to 
m or the number of discrete partitions summing to 2m with even elements. 
The foregoing analysis can also be extended by raising Eq. f ll96p to an 
arbitrary power p and applying Corollary 1 to Theorem 1. Then we obtain 
a power series expansion like Equivalence fl59|) . but in powers of z with the 
coefficients depending upon p. Hence, we arrive at 

oo oo 

Y[(^l^uJz'Y^l + Y^q(^k,-u,p)z'' , (209) 

k=l k=l 

where the coefficients q{k,u,p) can determined from Eq. fl^U]) with Di = 
q{i,co) and are given by 

qik, u, p) = Lp,, [(-1)^'= i-pU, n ^^^1 ■ (210) 

1=1 

It should be emphasised that in the above results p can be any value including 
a complex number. For integer values of p greater than zero the equivalence 
symbol can be replaced by an equals sign. As we shall see shortly, when 
p = — 1 and 0; = 1 in Equivalence fl209p . the q{k, —1, —1) equal the partition 
function or p{k), while if p is equal to a positive integer, say j, then Eq. fl210p 
simplifies drastically due to the fact that for k>j, the factor {—p)n,, vanishes 
for Nk>j. That is, the partitions with more than j elements in them do not 
contribute to the q{k,u,p). Moreover, if u = —l, further redundancy occurs 
in Eq. ( ]210p since the q{i,uj) become the discrete partition numbers or q{i), 
which we have seen are only non-zero when z is a pentagonal number [36] . 

There is also another approach to developing a power series expansion 
to the generating function given on the Ihs of Equivalence fl209p . This too 
involves the partition operator, but rather than consider the generating func- 
tion in Equivalence (I209p . we shall investigate the more general case of Eq. 
fll79p raised to an arbitrary power p as set out in the following theorem. 

Corollary to Theorem 3. The generalised version of the product in 
Equivalence fl209p whereby the coefficients of z'' are set equal to Ck can be 
expressed as a power series or generating form given by 

oo oo 

IKl + CkzY'^l + Y^BMz' , (211) 

fc=i fc=i 
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while the coefficients in the series are found to be 



Bkip) = Lp^k 



n 



Pijrii 



rii 



cr 



(212) 



In the above results (p) denotes (pi,P2, • • • ,Pk)- If 5' = inf ICjI""*^/* >0, then 
for |2;| < S* the equivalence symbol can be replaced by an equals sign. In 
addition, the coefficients Bk{p) satisfy the following relations: 



^P,k 



-if'NklH 



b,{pY^ 



i=l 



riA 



Lph 



\Nk TT \Pi)ni Qm 



riil 



(213) 



and 



Bkin + iy) = J2 Bjip) Bk-j{v) . 

j=0 



(214) 



Remark. The reader should observe that in previous cases involving a 
constant power of p the Pochhammer symbol appeared outside the product 
in the partition operator as in Corollary 1 to Theorem 1. Thus, the p- 
dependence of the coefficients in the resulting power series was only affected 
by the total number of elements in each partition. In the above case each 
element i is assigned a value that is dependent upon pi and consequently, 
{Pi)ni appears inside the product being acted upon by the partition operator. 

Proof. In order to prove the theorem, we shall use Lemma 1 again. Then 
the generating function can be expressed as 



n(i+Cfc^') 



pk 



-pi] 



k=l 



X 






^i-)0(i + E 



-pi] 



1\^ 



f- 



i-Cz^) 



i+E^T"(-c--T 1 + E 



i=i 



i=i 



j! 






"P4)j / ^ ^4^J 



(215) 



According to Lemma 1, each binomial series in the above product is ab- 
solutely convergent in the disk given by \z\ < \Ci\~^^\ Each series is also 
conditionally convergent in a specific region of the complex plane. From the 
proof of Theorem 2, which deals with the pi = —1 and Cj = —1 case, we 
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found that these regions can cancel each other when all the series appear in 
a product as in the above generating function. As a result, we were left with 
the region in which all the series are only absolutely convergent. Therefore, 
if S* = inf |Ci|~^/*>0, then all the series in Equivalence (I215P are absolutely 
convergent whenever |2;| <5' and the equivalence symbol can be replaced by 
an equals sign. Furthermore, when the pi all equal a positive integer, say /c, 
the binomial series become polynomials of degree k. Then the equivalence 
symbol can be replaced by an equals sign for any value of z. 

Expanding the product of all the series in Equivalence (12151) in powers of 
z yields 



k=l ^ 



''''^'cl-{-p,),cAz' 



2! 
+ {-^-^^Cl + {-p^),{-p,),C,C,-{-p,),C^z^ + --- . (216) 

From this result we see that there is one term appearing in the first order 
coefficient in z, while there are two terms appear in the second order co- 
efficient. The third order coefficient is composed of three terms. Had the 
fourth order term been displayed, there would have been five terms in the 
coefficient. In fact, the number of terms in the fc-th order coefficient is the 
number of partitions summing to k or p{k). Therefore, we need to develop 
a means of coding the partitions as was done for the partition method for a 
power series expansion in Theorem 1. 

Since the first order term corresponds to {1}, we assign a value of —Ci 
to each occurrence of a one in a partition. The first term in the second 
order coefficient possesses a factor of C^. Therefore, it must correspond to 
{1,1}, while the other term must correspond to the other partition summing 
to 2, viz. {2}. This term possesses a factor of — C2. Hence, each occurrence 
of a two in a partition yields a factor of — C2. If we continue this process 
indefinitely for the higher orders, then we find that each occurrence of an 
element i yields a factor —Ci. 

This, however, is not the complete story. Accompanying the factor of 
—Ci in the first order term is the factor of (— pi)i, while the first and second 
terms in the second order coefficient possess factors of (— pi)2/2! and (—^2)1, 
respectively. That is, when there is one occurrence of an element i m. a. 
partition, its assigned value must be multiplied by (— pj)i, but if there are 
two occurrences of an element in a partition, then the assigned value must 
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be multiplied by (— pj)2/2!. Therefore, if an element i occurs Ui times in a 
partition, then it contributes a factor of (— pj)„.(— Cj2;*)/nj! to the coefficient. 
This can be checked with the various terms comprising the third order term. 
The total contribution made by a partition is then given by the product over 
all elements, viz. Y[i=ii~Pi)ni{~Ci)'^^ /riil. Finally, the coefficient Bk{p) is 
evaluated by summing over all partitions summing to k. Hence, we arrive at 
Eq. ^rm- 

If we assume that z < S, which is actually not necessary, we can invert 
the equation form of Equivalence (12111) . thereby obtaining 

Because the rhs can be regarded as the regularised value of the geometric 
series, we can apply the method for a power series expansion where the 
coefficients of the inner series, viz. p^ in Theorem 1, are equal to —Bk{p), 
while the coefficients of the outer series qu are once again equal to (— 1)'^. 
Then we have 

oo oo 

'[[[l + Cuz'')-''' ^l + Y.Dkz'' , (218) 

fc=i fc=i 

where according to Eq. ( TTOl) . the coefficients of this expansion are given by 

k 

■ a 

rii 



Du 



^P,k 



{-if^mfi^^] . (219) 



i=l 



We also know from earlier in the proof that the above product can be ex- 
pressed in terms of a power series of the form given on the rhs of Equivalence 
(12111) except that now the pk are replaced by —pk- Because the resulting 
power series expansions possess the same regularised value, they are equal to 
one another. Moreover, since z is arbitrary, we can equate like powers of z, 
thereby yielding Eq. (I213p . 

The final identity is easily proved by noting that 

oo oo 

n (1 + CkzT'"^'' = n(l + C'^'T (1 + CkzT ■ (220) 

1=1 j=l 
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The Ihs represents the regularised value for the series on the rhs of Equiv- 
alence f l21ip with coefficients B^^dj, + u), while the rhs represents the regu- 
larised value of the product of two series, one with coefficients Bk{^i) and the 
other with coefficients Bklv). As the regularised value is the same in both 
situations, we have 



oo oo 



oo 



fc=l j=l k=l 

oo k 

= ^ + Y.^'llBMBk-A^) . (221) 

fc=i i=o 

Since z is arbitrary, we can equate like powers yet again. Therefore, we 
obtain Eq. fl214p . which completes the proof. 

In order to make the foregoing material clearer, we now consider a few 
examples. We have already generalised the discrete partition polynomi- 
als q{k,u) by introducing the power of p into their associated product as 
demonstrated by Equivalence fl2U9p . It was found that the new polynomials 
q{k, uj, p) could be expressed in terms of the partition operator acting on the 
discrete partition polynomials given by Eq. fl210p . Now we apply the corol- 
lary to Theorem 3 to the product with the Ck and pk set equal to u and p 
respectively. Consequently, we find that 

g(A;,a;,p) = Lp,,[(-a;)^'=nt^] . (222) 



So, let us evaluate g(3, w, p) via Eqs. fl210p and fl222p . In the case of Eq. 

we require ^(l,^), q{2,uj) and g(3,w), which have already been evaluated. 

Therefore, we obtain 

g(3, u, p) = -^^ ^' + (-P)2 u;' - (-P)i (c^' + ^) ■ (223) 

After a little algebra, the above becomes 

3 2 

g(3, a;, p) = (^ - ^ + u;3 ^ p2^2 ^ ^^ _ (224) 

By setting fc = 3 in Eq. (12220 . we obtain 

g(3, u, p) = -^^ u' + (-p)i (-p)i u^ - (-p)i u . (225) 
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Again after a little algebra, we end up with Eq. 

see that Eqs. (12231) and (12251) are composed of different quantities. 



although we can 
This 



demonstrates that while Eqs. (12101) and (I222p both possess the same sums 
over partitions, they are in fact different representations for q{k,u,p). 

The next example, which is a result attributed to Euler, appears on p. 56 
of Ref. CHI. This is 



JJ(1 - z'^Y = l-3z + 5z^-7z^ + --- = ^(-1)^(2A; + l)^^^') . (226) 



k=l 



From Equivalence 



we have 



fc=0 



n(i-^') =i+E^(^'-i'3) 



(227) 



fc=i 



fc=i 



Note that the equivalence symbol has been replaced by an equals sign because 
p = 3 in this case. Euler's result shows that only when the power of z is 
another type of figurate number called a triangular number [37| l38j . are the 
coefficients of the power series or generating function non-zero. If we equate 
like powers of z in both power series expansions given above, then we arrive 
at 



^P,k 



i-ir^ ( 



-3)a^. n 



q{i) 



Ui 




if A:=K^), 
otherwise. 



(228) 



On the other hand, putting Cj = — 1 and pj = 3 in Eq. ( I212p yields 



g(fc,-l,3) = L 



P,k 



[n 



(-3), 



Ui 



(229) 



In the above result (— 3)„- is only non-zero for rii equal to 1, 2 and 3, in which 
case it equals -3, 6 and -6, respectively. This is a different type of restricted 
partition from those we have encountered previously since it means that 
partitions in which an element appears more than three times are excluded. 
To generate such partitions all that needs to be done is to scan each partition 
twice by introducing another for loop in termgen of the partition generating 
program presented towards the end of Sec. 2. If in the first scan the rii or the 
variable freq is greater than 3, then a goto statement is required so that the 
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program avoids the next for loop, which is responsible for printing out the 
specific partitions. When such a code is constructed, one will find that out 
of the total of 627 partitions summing to 20, there are only 320 partitions in 
which all the elements occur at most three times. As a result of Eq. fl229p . 
we also arrive at 



Lp^k 



' Mki _ j(-i)^(2j + i), if fc=e+^). 



r-rrlzi^j^U-in^J + iJ, It A:=nj, (230) 

L-'--'- nA \ 0. otherwise. 



Moreover, it should be noted that if p = 1 such that n^ < 1 for all elements 
i in the partition, then the sum over all partitions in Eq. f l212p reduces to 
the discrete partition operator, viz. Lnp^kY], irrespective of the values for the 
coefficients Cj. 

If we put p= 1, take the cube power of the series on the rhs of Eq. fl209p 
and equate like powers of the resulting series with those on the rhs of Eq. 
f l227p . then we obtain 

k ji 

q{k, -1, 3) = 5^ ^ g(fc - ji)g(ji - j2)g(j2) , (231) 

ii=Oi2=o 

while multiplying the p = 1 and p = 2 versions of the series on the rhs of Eq. 
f l209|) yields 

g(fc,-l,3) = J2^ijMk-j,-l,2) = |(-l)^(2z + l), k= Ct;), (232) 
^—: 0, otherwise. 

In the above result g(2, u, 2) can be evaluated by putting p = 2 in Eq. f l210p . 
They can also be determined by equating like powers of z when taking the 
square of the series on the rhs of Equivalence f l209p with the p = 2 series. 
This gives 



g(A;,a;, 2) = N^g(z,a;)g(A; — z,a;) . (233) 

For 0; = —!, Eq. ( I233P reduces to 

k 

q{k,-l,2) = Y,<l{Mk-i) ■ (234) 
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k 


q{k,u,2) 


q{k,uj,3) 





1 


1 


1 


2u 


3a; 


2 


2uj + uj^ 


3u; + 3^2 


3 


2u + 4a;2 


3u; + 9^2 + ^3 


4 


2u + 5u^ + 20? 


3a; + 12a;2 ^ g^s 


5 


2u + 8u;2 + 4^3 


3w + ISw^ + IStu^ + 3w^ 


6 


2u; + 9u;2 + lOw^ ^ ^4 


3a; + 21a;2 ^ 37^3 ^ ^2a;^ 


7 


2u; + 12^2 ^ ^4^3 ^ 4^4 


3a; + 27a;2 + 54a;3 + 33a;^ + Sco^ 


8 


2u; + 13^2 ^ 22^3 ^ g^4 


-3a; + 30a;2 + Sla;^ + 66a;^ + 12a;5 


9 


2uj + 16^2 ^ 30^3 ^ ^g^4 


3a; + 36a;2 + 109a;=^ + 114a;^ + 39a;5 




2u;5 


+a;6 


10 


2u; + 17^2 + 40a;^ + 30a;^ 


3a; + 39a;2 + 144a;3 + 189a;^ + 81a;^ 




Au^ 


+9u^ 



Table 4: Coefficients q{k,u,2) and g(A;,a;,3) in the power series expansions 
of the p = 2 and p = 3 cases of Equivalence f l2U9p . 

Table m displays both q{k,u,2) and g(fc,a;,3) up to A; = 10. 

On p. 23 of Ref. |3T] it is stated that Gauss derived the following result: 



E 

k=0 



{k^+k)/2 



k=l \ 



(1 - Z^fe) 



fc\2 



y2k^ 







n(i-.^)(i+.^) 



(235) 



fc=l 



Once again, the equals sign is only valid for |2;| < 1. From this result we see 
that only the powers of z equal to (j^ + j)/2, where j is any non-negative 
integer, possess a non-zero coefficient. For any power k we obtain contribu- 
tions from all the partitions summing to k with only odd elements in them 
and from the discrete partitions summing to k consisting of only even ele- 
ments. If the number of even elements in these discrete partitions is even, 
then the partition will yield a value of unity. Otherwise, it will yield a value 
of -1. In addition, the partition function possesses mixed partitions com- 
posed of discrete partitions with only even elements and standard partitions 
with only odd elements. The values contributed by the mixed partitions to 
the coefficients of the power series depend upon the number of elements in 
the discrete partitions. For example, if we consider z^^ or ^ = 4 on the Ihs of 
Eq. fl235p . then it will be composed of the contributions due to the partitions 
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summing to 10 with only odd elements. There are 10 of these beginning with 
{1,9} and ending with {5,5}. Hence, these partitions contribute a value of 
10 to the coefficient. On the other hand, there are only 3 discrete partitions 
with even numbers summing to 10. These are {10}, {2,8} and {4,6}. Since 
the last two possess an even number of elements, they each contribute a 
value of unity, while the single element partition gives a value of -1. Overall, 
the discrete partitions summing to 10 contribute a value of unity to the the 
coefficient, which now becomes 11 when the contribution from the standard 
partitions summing to 10 with only odd elements is included. However, when 
the discrete partition is {2}, we need to consider the partitions summing to 8 
with odd elements. There are six of these, beginning with {1,7} and ending 
with {3,5}. Because there is only one element in the discrete partition, these 
mixed partitions contribute a value of -6 to the coefficient, which now drops 
to 5. However, there are still more mixed partitions. We need to consider the 
discrete partition of {4}. In this case we need the partitions summing to 6 
with odd elements. There are four of these, beginning with {1,5} and ending 
with {3,3}. In this instance the mixed partitions contribute a value of -4, 
yielding a value of 1 for the coefficient as indicated above. We do not need 
to consider the contributions where the discrete partitions sum to either 6 
or 8 because in these cases there are only 2 partitions, one of which has two 
elements and one with one element. Hence, they cancel each other yielding 
a value of 0. 

If we introduce Eq. fl209p into the Eq. fl235p . then we find that 

CO OD k 

J2 z''+'/' = 1 + E ^' E ^(-^X^ - ^■' 1' 2) ■ (236) 

k=0 k=0 j=0 

By equating like powers of z, we obtain for i, a positive integer, 

2_^q{j)q{k-j,l,2) = < . (237) 

^—: 0, otherwise. 

Therefore, when k is not a triangular number [37], by combining the above 
result with Eq. fl232p we arrive at 

k 

E ^(■^■) (^(^ - ^■' 1' 2) ± (lik - J, -1, 2)) = . (238) 



114 



8 Other Products 

By using the material of the previous section we are in a position to study 
more advanced products. We begin by introducing the variable u next to 
the power of z in the denominator of P{z). Based on the similar extension 
of the product yielding discrete partitions in the previous section, we expect 
to obtain polynomials as the coefficients of the resulting generating function. 
Therefore, we define the new product as 



oo ^ 

Hn ^ • (239) 

-'■-'■ 1 — UJZ'^) 



P{z,u) = 

k=l 

Obviously, when a; = 1, this reverts to the generating function of the partition 
function given by Equivalence fl73|) , which we have seen becomes an equation 
when |-z| < 1. According to p. 112 of Ref. [18], the product in the above result 
can be written alternatively as 

oo ^ ^2 

while inversion of the rhs of Eq. ( 1196^ yields 

We have seen that P{z) or u; = 1 in Eq. fl239p yields a power series ex- 
pansion whose coefficients are given by the partition function or p{k). This 
expansion is obtained by expanding each term in the generating function 
into the geometric series for each value of i. It is this value in the generating 
function, which is responsible for yielding the specific elements in a parti- 
tion, while the power of z^ in each geometric series represents the frequency 
or number of occurrences of the element in the partition. For example, mul- 
tiplying (-2^)^ in the expansion of 1/(1 — 2;^) by (z^)"^ in the expansion of 
1/(1 — z'^) means that the partition has three twos and four threes in it. By 
introducing u into the generating function as indicated above, we see that the 
overall power of u yields the total number of elements in a partition. In the 
example just mentioned we now obtain [uz"^)^ multiplied by (uz^)"^, which 
yields u'^z^^. The power of 7 on w represents the total number of twos and 
threes in the partition. Therefore, we expect the coefficient of each power of 
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uj in the coefficients of the resulting generating function to indicate the total 
number of partitions where the number of elements equals the power of u. 

By regarding the rhs as the regularised value of the geometric series and 
the g(fc, —oj) as the inner series whose coefficients are pk-, we can apply The- 
orem 1 to the last form for P{z, uj) with y = z. If we express the generating 
function as a power series expansion with coefficients, p{k,u), i.e. 



oo .. oo 

fc=l ^ ^ k=l 



(242) 



then according to Eq. (jTOl) the coefficients of the resulting power series ex- 
pansion are given by 



p{k,u) = Lp^t 



'if^mll 



q[t,-UJ) 



i=l 



nA 



(243) 



Furthermore, if we put p= —1 and Cj = —u in the corollary to Theorem 3, 
then we observe that the p{k, u) become the coefficients Bk{—1) given by Eq. 
(12121) . Thus, we arrive at 



p{k,uj) 



^P.k 



UJ 



Nk 



(244) 



This tells us that the coefficients in the p(/c, u) will be the number of par- 
titions summing to k where each power of u corresponds to the number of 
elements in the partitions. That is. 



p{k,uj) = ^^Piik) 



UJ 



(245) 



j=i 



where from Sec. 2, Pi{k) 



. It has already been stated the sub-partition 



numbers obey the recurrence relation given by Eq. ([T]). In terms of the fixed 
number of elements operator defined by Eq. (IHHj) we also find that 



Pi{k) = L 



P,k 



(246) 



The first few partition function polynomials are found to be: p(0,C(;) = l, 

p{l,uj) = u, p{2,u) = u'^ + u, p{3,uj) = u^ + u'^ + u and p{4:,uj) = u'^ + u^ + 
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IlkP'^uj. In fact, those up to A; = 10 are displayed in the third column of 
Table 121 Since the coefficients of these polynomials represent the number 
of partitions in which the number of elements is given by the power of a;, 
the highest order term of these polynomials is fc, which arises from the k- 
element partition of {1,1,. . . ,lfe}- The other partitions are unable to provide 
an bj^ term because the highest order of all other g(fc, iS) is less than h. The 
partition {1,1,. . . ,lfc_i,2} only produces an ij^~^ term as its highest order 
term because the highest order term of (3'(2,a;) is 1. Hence, deg p{k,u) = k. 
Conversely, the lowest order term in the p{k, u) is the lowest order term 
in q{k,u) stemming from the single element partition {k}, which is unity. 
Therefore, p^^k) =pk-i{k) =pi{k) = 1, p^,_2{k) = 2 and p{k,l) = Y!l^-^pi{k) = 
p{k). Moreover, the total number of partitions with even elements is given 
t'y ^i=i P2i{k), while the total number of partitions with odd elements is 
equal to Yl^i P2i-i{k), where m = k/2 when k is even and m= [k/2] + 1 when 
k is odd. On the other hand, the coefficient P2{k) can be evaluated by noting 
that it represents the product of the two lowest order terms in each partition, 
namely {j,k — j}, where j ranges from 1 to [k/2]. According to Ref. [39] , 
the number of two-element partitions summing to k is given by P2{k) = [k/2], 
while the number of three-element partitions is given by Pslk) = [k'^/12] for 
k > 3. Moreover, a table of the sub-partition numbers for k and i ranging 
from to 11 is presented on p. 46 of Ref. [18]. All these results agree with 
those obtained via Eq. (11811) . confirming that the latter result does yield the 
number of partitions summing to k with i elements in them. 

An interesting property of the partition number polynomials can be de- 
rived by setting z = z^ and u = u'^ in Equivalence fl240p . Then we obtain 

oo ^ oo 

P (-'. -') = n a-lv^) - 1 + E H*. -') --" ■ (247) 

fc=l ^ ^ fc=l 

The quantity on the Ihs of the above equivalence can also be written 

P{z'^,uj^) =P{z,io)P{z,-uj) . (248) 

Introducing Equivalence (12421) into the above equation yields 

oo oo 

P{z,u)Piz,-u)^(l + Y,p{k,u)z'^[l + J2pik,-oo)z'^ . (249) 

fc=i fc=i 

Since the series on the rhs's of Equivalences (I247P and (I249P are derived from 
the identity given by Eq. (I248p . they are equal to one another in accordance 
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with the concept of regularisation [21 El El IS] • Because z is arbitrary, once 
again we can equate hke powers of z. Therefore, we obtain 



2fc+l 



^p(j,a;)p(2A; + l-j,-c^)=0 , 



(250) 



and 



2k 



^p{j,u)p{2k-j,-uj) =p{2k,uj'^) 

j=0 



(251) 



As in the previous section we can generahse the foregoing analysis by 
replacing —u in P{z,u) by C,, which effectively represents the inversion of 
Eq. 01791) . From Theorem 2 we know that 



Hiz) 






(1 _ c,z + Clz^ + . . . ) (1 - C^z^ + Clz' + ■ ■ ■ ) 

(252) 



(1 - C^z'' + C|z« + . . . ) (1 _ C4Z^ + C^« + ■ ■ ■ ) ■ 



where the equivalence symbol can be replaced by an equals sign provided 
|Cj2;*| <1 for all i. Expanding the above yields 



E{z)-^ = l-Ciz+ [Cl - C2)z' + {-Cl + C2C1- C^)z^ + 



(253) 



Each coefficient of z^ in the above power series is composed of contributions 
that can be related to the partitions summing to k as was the case in the proof 
of Theorem 1. The major difference between the above situation and that in 
the proof of Theorem 1 is that there is no multinomial factor associated with 
each contribution made by a partition as we also found in the proof to the 
corollary to Theorem 3. Furthermore, each element i in a partition is now 
assigned a value of —Ci so that the above result becomes 



H{z)-'^Y.^^'' 



(254) 



fe=0 



where Hn = l, and 



Hu 



^p± 



-D^.jjcr 



(255) 



j=i 
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The above result represents the case when the pk in the corollary to Theorem 
3 are set equal to -1. 

If the Ci are set equal to unity in Equivalence ( 12521) , then we find that 



n 



\ (1 + ^) 



1 + ^ Lp,fe 



-1) 



Nu 



(256) 



k=l 



where the equivalence symbol can be replaced by an equals sign for |z| < 1. 
Therefore, the coefficients on the power series expansion on the rhs represent 
the difference between the number of even- and odd-element partitions sum- 
ming to k. In addition, the above equivalence is analogous to putting the Ci 
equal to unity in Eq. fll79p and applying Theorem 1 to its inverted form. In 
this case the coefficients of the inner series are given by pk = -^i:)P,fc[l]; the 
number of discrete partitions, while the coefficients of the outer series are 
given by qt = {—^)^- Then the coefficients of the power series expansion 
on the rhs of Equivalence fl256p can be expressed in terms of the number of 
discrete partitions as 



^p,fc 



-1) 



Ni. 



^P,k 



N,\l[LDP,[ir/n,\ 



i=l 



If we multiply the product on the Ihs of Equivalence 
both the numerator and denominator, then we find that 



(257) 



by (1 — z* 



m 



n 



1 



\ (1 + ^0 



1=0 



— z 



2i+l 



) 



(258) 



When the product on the rhs of Eq. (12581) appeared in powers of z"^ rather 
than 2;^*+^, we saw that the resulting power series possessed coefficients which 
were equal to the number of discrete partitions summing to k. In the above 
result all the even powers are now missing. This means that the coefficients 
of the resulting power series will be the number of discrete partitions with 
only odd elements in them. That is. 



n(i--"^')=i+E(-i)'^«^^.' 



i=0 



A;=0 



(259) 



where ODP denotes that only partitions with odd elements are to be con- 
sidered in the sum over partitions. That is, n2i = for all values of i. The 
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phase factor of (—1)'^ in the series expansion arises from the fact that only 
an even number of odd discrete elements yields an even power of z, while 
only an odd number of discrete elements yields an odd power of z. That is, 

[fc/2] 



{-lfLoDP,k\l]=LDP^k\\{{-lT''-'\ . (260) 

4 = 1 

In a similar fashion we arrive at 

CXD OO 

\[{l + Z^')=l + Y,LEDP,2k\lY'' , (261) 



4 = fc = l 



where only even elements are to be considered in the even discrete partition 
operator, i.e n2i+i = for all i. Alternatively, we can replace the even discrete 
partition operator by the discrete partition operator since 



^EDP,2k 



[I] =LDP,k[l] . (262) 



Moreover, by equating like powers of z in the power series on both rhs's of 
Equivalences f l256p and fl259p . we arrive at 



^P,k 



{-If" = {-lfLoDP,k 1 . (263) 



From this result we see that when k is even, the number of even partitions 
is greater than the number of odd partitions, while for odd values of k, 
the opposite applies. Multiplying both sides by (—1)^ results in taking the 
absolute value or modulus of the Ihs. Thus, the above statement tells us 
that absolute value of the difference between the number of even and odd 
partitions is equal to the number of discrete partitions with only odd elements 
in them or the number of partitions with distinct odd parts, a result first 
proved by Euler according to p. 14 of Ref. [3T] . 

As a consequence of the previous section, it is a relatively simple exer- 
cise to produce a code that evaluates the difference between the number of 
partitions with an even number of elements and those with an odd number 
of elements. Two new global variables are required, one for evaluating the 
difference as each partition is scanned and another that is either equal to 1 
or -1 depending on whether there is an even number of elements or an odd 
number of elements. Once the second value is determined, it needs to be 
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Clut[4]= 
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Figure 5: The difference between even and odd partitions to tlie total number 
of partitions summing to k versus k 



Lp,k 


'{-lf\ 


> 


Lp, 


'i-lfj 



added to the first global variable in the main function prototype. The sec- 
ond global variable must be evaluated in the termgen function prototype 
after the for loop has been altered to calculate the total number elements in 
the partition, which is determined by summing all components of part. As a 
consequence, one finds after running the code for several values of k that 



(264) 



for k>j. Fig. O presents the graph of the ratio of the absolute value of the 
difference between odd and even partitions summing to k to the total number 
of partitions or p{k) for /c<50. Whilst the absolute value of the Lp^k[{—^)^''] 
increases with k, we see that in relation to the total number of partitions the 
ratio decreases monotonically, once k exceeds 15. 

From Eqs. fll92p and fll94p we see that the infinite products, Ilfclill + 
z'') and Y[T=i (-'- ~ ^^) ^^^ ^^ expressed in terms of a series in successive 
powers of z, where the coefficients are equal to Lpip^ [1] and Lpip^ [(—1)^'=] 
respectively. The first product can also be written as 



fc=i fc=i ^ 



2fc-l^ 



(265) 
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Eq. (12651) is easily obtained by manipulating the rlis after multiplying it by 
(l — z"^^) / (l ~ ^'^'^) ■ ■'■f ^^ P^^ ^i ~ ~^ i^ Equivalence f l252p . then according 
to Equivalence fl254p and Eq. fl255p we arrive at 



n(i+^M-i+Ei 



OEP,k 



fc=l 



fc=l 



(266) 



where LoEP,k[-] represents the odd element partition operator, which we have 
seen has two different forms given by Eqs. fll52p and (11530 depending upon 
whether k is an even or odd number. Because the generating function on 
the Ihs also yields a power series whose coefficients represent the number of 
discrete partitions summing to A;, we see immediately that 



^DP,k 



L 



OEP,k 



(267) 



Hence, the number of discrete partitions is equal to the number of partitions 
composed only of odd elements, another result attributed to Euler according 
top. 5of Ref. taij. 

By putting Q = -1 in Eq. flT79|) . we obtained Eq. flTMl) . If we put 
Ci = —1 and z = z^, then the product on the Ihs of Eq. (I194p yields a series 
expansion in powers of z^, but now the coefficients represent the number 
of distinct partitions summing to 2k with only even elements operating on 
(— l)^2fe^ Alternatively, this is equivalent to the number of distinct partitions 

-1)^*. Therefore, we can write 



summing to k operating on 

oo 

nil 



fc=i 



2fc 



IH-Ei 

fc=l 

oo 

1+Ei 



EDP,2k 



DP,k 



;-i) 



\N2k 



Nk 



^2fc 



.2fc 



(268) 



k=l 



In Eq. (I268P LEDP,2k[-] denotes the even discrete partition operator, which 
acts on the number of discrete partitions summing to 2k where the elements 
are only even integers. This is opposed to the odd discrete partition operator, 
LoDP,k[']^ where the elements are odd numbers and discrete, but can sum to 
both even and odd integers. Since the infinite product of (1 - 
product of two separate infinite products involving {1 — z'') and 
have 



z'^'') is the 
'1 + z''), we 



n(i-^") = (i+Ei 



WP,k 



-1) 



Nk 



i+Ei 



DP,k 



(269) 



k=l 



k=l 



k=l 
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Equating like powers of z in Eq. 

2fe 



with those in Eq. 



yields 



j=0 



N, 



WP,2k-j 



(270) 



and 



2fc+l 

i=o 



\N, 



WP,2k+l~j 



. 



(271) 



Since L£)p^fc[(— 1)^*=] and LDP,k[i] are equal to q{k) and q{k,l) respectively, 
Eqs. (12701) and (12711) can also be written as 



fe-i 



q{k) (l - q{k, 1)) = J2 {(lUH^k - J, 1) + q{2k - j)q{j, 1)) , (272) 



j=0 



and 



J2 {qUWk + 1 - J, 1) + g(2A; + 1 - j)q{j, 1)) = . 



(273) 



i=o 



In these results it should be borne in mind that g(0) =q{0, 1) = 1. Isolating 
the j = terms in the above equations yields 



fc-i 



g(2/c, 1) + qi2k) = q{k) (l - q{k, 1)) - ^ (g(j) qi2k - j, 1) 

+ g(2A;-j)g(j,l)) , (274) 



and 



q{2k + 1,1) + g(2A; + 1) = J^ (g(j) q{2k + 1 - j, 1) 

+ g(2A; + l-j)g(j,l)) • 



(275) 



Eqs. (12741) and (12751) represent the recurrence relations for determining the 
number of discrete partitions or g(fc, 1). Like the Euler/MacMahon recur- 
rence relation given by Eq. (I165p . they utilise the discrete partition numbers 
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or q{k) and consequently, most of the terms in the sums vanish when the 
summation index j is not equal to a pentagonal number. 

It should also be noted that the analysis resulting in Eq. f ll69p can be 
adapted to provide another representation for the partition number polyno- 
mials or p{k,u). First, we re-write the generalised product in Equivalence 
as 



oo oo 



P(z,a;) = exp(5^5^a;^z-Vj) , (276) 

m=l jr'=l 

where now it is assumed that |a;2;| < 1. Consequently, the modified version of 
Eq. fll69p becomes 

oo _ 

P{z, u) = 1 + ^ — (uz + {u + uj'^/2) z^ + {uj + wVS) z^ 

k=l 

+ {u + u?l1 + a;74) ^^ + ■ ■ ■ + ^j{uj)z^ + ■ ■ ■ )\ (277) 

where '^j[uj) = J2d\ji'^/J)'^''^^ ^^^ ^ represents a divisor of j as before. The 
first few divisor polynomials are: 70(1^) = 1, 71(0;) = w, 72(0;) = uj + ^"^ l'2.i 
and 73(0;) = (jj + cu^/S. As expected, for w = 1 the 7j(co') reduce to the 7^ 
below Eq. (I167p . This means that Eq. (I170p can be generalised to 

q{k,-uj) = LpA{-lf>'\{^-^] . (278) 

1=1 

Furthermore, the 7j(c(;) are polynomials in u whose highest and lowest orders 
are respectively j and unity. That is, deg 7j(w) = j. We shall refer to these 
unusual polynomials as the divisor polynomials. By applying Theorem 1 to 
Eq. (I277P with the coefficients of the inner and outer series set equal to '~ik{<^) 
and l/k\ respectively, we arrive at 

p(A;,a;) = Lp,Jn^^j . (279) 

4 = 1 

Since deg 7^(0;) = i and nj=i '^*"' =uj^, the highest order term in the p(fc, u) 
is k. On the other hand, since the lowest order term in the 7j(a;) is w, 

124 



Y[i=i ^^' = u)'^'' and Nk ranges from unity to k, the lowest order term in the 
p{k, u}) is unity, again confirming that the partition number polynomials are 
polynomials in u with deg p{k,u) = k. 

The sixth and last program presented in the appendix is called dispfnpoly. 
It prints out both q{k, —u) and p{k, u) in symbolic form so that it can be 
handled by Mathematica [12]. To run this program, the user must specify 
the order k of the polynomials. The program is different from the other 
programs in the appendix since it is not required to determine the factorial 
of the total number of distinct parts, i.e. iV^!. When the global variable 
polytype is equal to unity in the for loop in main, the program determines 
the discrete partition polynomial for the order specified by the user. When 
it becomes equal to 2, the program determines the corresponding partition 
function polynomial. E.g., for k = 6, the following output is generated: 
Q[6,-w_]:= DP[6,w] (-1) + DP[l,w] DP[5,w] (-1)^(2) + 
DP[l,w]^(2) DP[4,w] (-l)^(3)/2! + DP[l,w]^(3) DP[3,w] (-l)^(4)/3! + 
DP[l,w]^(4) DP[2,w] (-l)^(5)/4! + 

DP[l,w]^(6) (-l)^(6)/6! + DP[l,w]^(2) DP[2,w]^ (2) (-l)^(4)/(2! 2!) + 
DP[l,w] DP[2,w] DP[3,w] (-1)''(3) + 
DP[2,w] DP[4,w] (-1)^(2) +DP[2,w]^(3) (-l)^(3)/3! +DP[3,w]^(2) (-l)^(2)/2! 

P[6,w_]:= DP[6,w] + DP[l,w] DP[5,w] + 

DP[l,w]^(2) DP[4,w] /2! + DP[l,w]^(3) DP[3,w] /3! + DP[l,w]^(4) DP[2,w] 
/4! + 

DP[l,w]^(6) /6! + DP[l,w]^(2) DP[2,w]^(2) /(2! 2!) + DP[l,w] DP[2,w] 
DP[3,w] + 

DP[2,w] DP[4,w] + DP[2,w]^(3) /3! + DP[3,w]^(2) /2! 
Time taken to compute the coefficients is 0.000000 seconds 
The terms on the rhs denoted by DP[fc,w] represent the divisor polynomials 
of order k. These can be obtained by typing in the following line in Mathe- 
matica: 

DP[k_,w_]:= SumKd/d,{d,Divisors[k]}]. 

From Eqs. f l278p and (I279p we see that the discrete partition polyno- 
mials or rather the q{k, —u) are almost identical to the partition function 
polynomials or p{k, u) when they are both expressed in terms of the divisor 
polynomials, the only difference being a phase factor that appears in the 
former. This factor is positive when the number of elements in a partition is 
even, but is negative when there is an odd number of elements in a partition. 
Hence, the only difference between the two sets of polynomials occurs for 
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odd partitions or those with an odd number of elements. 

If we multiply Eq. f ll96p with Eq. (I242p , then by equating like powers on 



(280) 



both sides of the resulting equation we obtain 

k 
j=0 



Alternatively, Eq. fl280p can be expressed as 

k j / \„- k-j 



j=0 



i=l 












(281) 



These results represent the generalisation of the Euler/MacMahon recurrence 
relation given by Eq. (I165p . 

We can also generalise the product in Equivalence (12420 by introducing 
the arbitrary power of p as we did in the discrete partition case of Equivalence 
(I209p . Thus, the generalised product becomes 



Pp{z,u) 



oo ^ 

k=l ^ 



- uz'')p 



i + X^p(^,w,p) 



(282) 



fc=i 



In this case the coefficients p{k,u,p) in the generating function can be de- 
termined by setting the Di equal to p{i,uj) in Eq. (!60|) . Then one finds that 



p(k,u,p) 



^P,k 



'-^r^ {-p)n, n 



p{i,u) 



1=1 



nA 



(283) 



As a consequence of the preceding analysis, we are now in a position to 
study more advanced products. In particular, let us consider the following 
quotient: 



°° /'i a k\ 

P{z, (3u, au) = Q{z, -/3u)P{^, «^) = TT i ~ J 



(284) 



fc=i 



In deriving a power series expansion or generating function for the above 
product we expect the power of u in the coefficients to yield the total number 
of elements in the partitions summing to the power of z. Furthermore, the 
power of P in the coefficients should represent the number of elements due 
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k 


QPk{u;,a,f3) 





1 


1 


(a - P)u 


2 


{a - /3){uj + aoo'') 


3 


{a - ^){uj + {a - /3)uj^ + a^uj^) 


4 


{a - /3)(a; + (2a - P)u^ + a{a - l)u^ + a^u^) 


5 


{a - /3){uj + 2{a - /3)uj^ + 2a{a - /3)io^ + a\a - (3)io^ + aV) 


6 


{a -(3)luj + (3a - 2(3)uj^ + {3a^ - 3a/3 + /3^)uj^ + 2a\a - (3)uj^ 




+ a3(a - /3)u;^ + a^w^) 


7 


(a -(3)iuj + 3(a - (3)lo^ + (4a - /3)(a - /3)uj^ + (3a - /3)(a - /3)aw^ 




+ 2a3(a - /3)uj^ + a\a - (3)uj^ + a^w^) 


8 


(a - /3)(w + (4a - 3/3)^2 ^ (5^ _ 2/3) (« _ /3)^3 ^ (5^2 _ g^^ ^ 2/32)aa;^ 




+ a2(3a - /3)(a - /3)u;5 + 2a^(a - /3)u;^ + a%a - /3)a;^ + a^t^«) 



Table 5: Coefficients QPk{u}, P, a) in the power series expansion for the gen- 
erating function given by Eq. (I280p . 

to the discrete partitions, while the power of a should give the number of 
elements due to the standard partitions. By adopting the same approach 
as for the other infinite products that have already been presented in this 
section, we can express Eq. (I284p as 

00 00 

P{z, l3u, au) = ^2 li^^ ~l3uj)z^ ^ p{k, auj)z'' 

k=0 k=0 

00 

= ^z'=QPfc(c^,/3,a) (285) 

fc=0 

where 

k 

QPk{u,P,a) = J2(lij^-Mpik-J,aoo) . (286) 

j=0 

From Eq. (I286p we see that QPo{u, /3, a) = 1. Furthermore, since P{z, au, au) 
equals unity, it follows that QPk{uj,a,a)=0 for k>0. 

Table O presents the coefficients QPkioj, f3,a) up till k = 8. As can be 
seen from the table they are polynomials of 0{k) in u. The power of u 
in these polynomials gives the number of elements in the final partitions. 
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which combine the elements from standard partitions with those from discrete 
partitions. As expected, the polynomials vanish when a = f3 since the product 
P{z, (3u, au) equals unity in this case. Furthermore, the highest power of a 
is k, which is also the highest power of u. This corresponds to the fact that 
the power of a represents the number of elements in the standard partitions. 
Therefore, the greatest number of elements in the final partitions will be 
due to the partition {1,1,. ..,1^} with no elements coming from a discrete 
partition. The highest power of (3, however, is considerably lower since it 
is determined by the partition with the most number of discrete elements 
summing to k. In this instance the power of a will be zero. E.g., for k = 8, 
the highest power of /3 is three, which is in accordance with Eq. fl20ip . When 
a = 0, the polynomials reduce to the polynomials arising from the generating 
function for discrete partitions, i.e. q{k,—Pu), while for /3 = 0, they reduce 
to the partition function polynomials or p{k,au). In addition, for a = the 
coefficients in the resulting polynomials give the number of discrete partitions 
where the number of elements is equal to the power of /3. For /3 = the 
coefficients of the resulting polynomials become the sub-partition numbers 
or Pi{k), where i represents the power of a. 

The interesting terms in the polynomials displayed in Table [5] are the 
cross-terms involving a and /9, which represent the mixture of the discrete 
partitions and standard partitions with the total number of the elements 
equal to the power of co. For example, in QP^^Uja, (3) the coefficient of u^ 
has a term equal to —2a(3, which tells us that one element in the partition 
has come from the generating function for standard partitions and the other 
has come from the generating function for discrete partitions. There are two 
instances where this can occur: either the one has come from the discrete 
part or numerator on the rhs of Eq. fl280p and the two from the denominator 
or vice- versa. On the other hand, there is only one instance of the partition 
{1,2} emanating only from either the numerator or the denominator. Thus, 
the coefficients of o;^ in QP^^Uja, (3) for only standard and discrete parti- 
tions are respectively a^ and /3^. In this instance the partition maintains its 
discreteness when accepting an element from the discrete partitions and one 
from the standard partitions, but this will not always be the case. In addi- 
tion, the power of a can be much higher than the power of /3 reflecting the 
fact that the greatest number of elements in a discrete partition summing to 
a particular value is significantly less than the greatest number of elements 
in a standard partition summing to the same value. 
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According to p. 23 of Ref. [21], Gauss derived the following result: 



2D-l)'-'' = nfe^- (287) 



fc=l A;=l 



The rhs of the above result is a special case of P(2;, /3a;, aw), namely P(2;, 1, — 1) 
If the values for au and /Sw are introduced into the rhs of Eq. fl285p . then we 
can equate like powers of z with the Ihs of the Eq. (12871) . Consequently, for 
i equal to a positive integer we arrive at 

gP.(l, 1, -1) = X: q{j) P{k - J, -1) = 1^^"^)^' ^ = ''' . (288) 

■^ U, otherwise . 

Eq. (12881) can be checked with the results appearing in Table [51 

We can also use the preceding analysis to derive a power series expansion 
or generating functon for the product of two specific forms of P{z,x,y) in- 
volving the three parameters, u, x, and y, and the variable, z. This product 
was first studied by Heine. According to p. 55 of Ref. [TS] he found that 

11 il-coz'^) (l-ujxyz^)-^^ (z:z), (ujz:z\^, ^''''^'^ ' ^^''^^ 



fc=i (l-^^'') {^-^^yz^) ^ {z;z)^, {ujz;z) 



If we set a = 1/x, b = 1/y, c = uz and q = z with \c/ab\ < 1 and |g| < 1, 
which are the conditions for guaranteeing absolute convergence, then the 
above result can be expressed as a q-hypergeometric series. This is perhaps 
the more familiar form for the product, where it is written as 

^ (g; q), {b; g). / c x ^ -j^ (l - {c/a)q'') (l - {c/b)q') 
2^^{q;q),{c;q),\ab) 11 {l - cq'^) {l-{c/ab)q^) ' ^ ^ 



This result appears as Corollary 2.4 on p. 20 of Ref. 

The Ihs of Eq. (I289p represents the product of P(2;, ux, co) and P{z, uy, uoxy). 
If we denote the product of P{z, x, y) and P(z, s, t) by P2(-2, x, y, s, t) and in- 
troduce Eq. (12851) . then we find that 

(1 — uoxz^) (1 — ujyz^) 



P2{z,ujx,uj,ujy,ujxy) = W—-- 

^J- (1 — coz'^) (1 — uxyz'^j 

OO 

= 5^;2'=i7Pfe(a;,x,i/) , (291) 



fc=0 
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k 


HPk{u,x,y) 




1 


1 

ix-l)iy-l)u 


2 


ix-l){y-l)uj{l + il + xy)u) 


3 


{x - l)(y - l)u;(l + {x-l){y- l)co + (1 + xy + x^y^)u^) 


4 


{x-l){y-l)uj{l + {2~x-y + 2xy)io + {x - l){y - 1)(1 + xy)io^ 




+ (1 +xy + x^y^ + xV)^^) 


5 


(x - l)(y - 1V(1 + 2{x - l){y - l)uj + 2(x - l)(y - 1)(1 + xy)u^ 




+ {x -l){y- 1)(1 + xy + x'^y^ + x^y^)u^ + {1 + xy + x^y^ + xV 




+xY)oo^) 


6 


{x - l)(y - l)uj{l + {3-2x-2y + 3xy)uj + {3 - 3x - 3y + x"^ + y^ 




+21xy - 3xy'^ - 3x'^y + 3x^y^)uj^ + {x - l){y - 1)(2 + 3xy + 2x^y^)u^ 




+ {x -l){y - 1)(1 + xy + x^y^ + xV)^'^ + (1 + xy + xV + a^V 




+x^y^ + x5i/5)w5) 



Table 6: Coefficients of tlie polynomials HPk{u, a, /3) arising from the three- 
parameter one variable generating function given in Eq. (12911) . 

where 

k 

HPkioo, x,y) = J2 W(^' a;, 1) QPk-j{oo, y, xy) . (292) 

3=0 

The equals sign appears here because of the conditions given below Eq. (I289p . 
Table [6] presents the coefficients HPk{uj,x,y) up to fc = 6. They have 
been obtained by implementing Eq. (I292p in Mathematica. From the table 
it can be seen that the HPk{u,x,y) are polynomials in u of degree k. The 
coefficient of the leading order term is 

cr= ^\[fjy^\ ^-i)iy-i) , (293) 

while that for the penultimate leading order term is found to be 

Cl'. = ^ ., \ ^ (X - my - 1? . (294) 

(1 - xy) 

In the above equation k>2. The lowest order term in u for these polynomials 
is linear and its coefficient is equal to 

Cr = (x-l)(y-l) . (295) 
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As expected, the polynomials are zero when either x= 1 or y = 1 since in these 
cases P{z, ux, u, uy, coxy) is equal to unity. In addition, they are symmetrical 
in X and y in the sense that any power of u with a term of ax^y\ where i^j, 
in its coefficient will also possess the term of ax^y^ with the same power of 

UJ. 

9 Conclusion 

Originally, this work set out to devise a programming methodology on the 
partition method for a power series expansion, which has been used recently 
to solve important and intractable problems in applied mathematics [I]-[l]. 
In this method the coefficients of the resulting power series for a function 
are obtained by summing the contributions made by each partition summing 
to the order k. These contributions are evaluated by: (1) assigning values 
Pi to each element i in a partition, (2) multiplying by a multinomial factor 
composed of the factorial of the total number of elements in the partition, N^l 
divided by the factorial of each element's frequency, nj and (3) multiplying 
by the coefficient of an outer series for the total number of elements in the 
partition, viz. q^i^. In order to apply the method, it means that one needs to 
know the composition of all the partitions summing to k, which includes the 
frequencies or numbers of occurrences of all the elements in each partition. 
Therefore, an algorithm is required that is capable of scanning all these 
partitions, the number of which increase exponentially with k. Whilst Sec. 
2 discusses various methods of generating partitions, it turns out that the 
novel bi-variate recursive central partition or BRCP algorithm is the most 
suitable method for implementation in the partition method for a power 
series expansion because it is based on the graphical representation of the 
partitions in the form of a non-binary tree diagram as depicted for A; = 6 in 
Fig. [1] As a consequence, the BRCP algorithm is able to print out partitions 
in the multiplicity representation more efficiently than the other algorithms 
discussed in Sec. 2, while the multiplicity representation turns out to be the 
minimum amount of information required for carrying out the method for a 
power series expansion. 

The theory behind the partition method for a power series expansion is 
presented in Sec. 3 as Theorem 1, which shows how power series expansions 
can be derived from a quotient of pseudo-composite functions. It also rep- 
resents the lynchpin of this work. Next the regularisation of the binomial 
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series is presented in Lemma 1. With this result a corollary to Theorem 1 
is developed, whereby the partition method is adapted to the situation in 
which the quotient of the pseudo-composite functions can be taken to an 
arbitrary power. As a result of Theorem 1 and its corollary, we observe that 
the process of evaluating the contributions made by each partition can be 
viewed as a discrete operation, giving rise to a partition operator denoted 
by Lpfe[-]. While Lp,A:[l] =p{k) or the number of partitions summing to k, 
varying the argument inside the operator yields completely different identi- 
ties. Moreover, the partition operator can be modified so that it only applies 
to specific types of partitions such as discrete or odd/even partitions, again 
resulting in further new and fascinating identities when the arguments are 
altered. 

Because the number of partitions increases exponentially, it becomes 
rather onerous to apply the partition method for a power series expansion 
when the order k is greater than 10. This problem is overcome by modifying 
the BRCP algorithm to calculate the contribution due to each partition in 
symbolic form. Sec. 4 presents two such programs. The first partmeth cal- 
culates all the coefficients Dj^ and Ej. in Theorem 1 up to and including the 
value of k specified by the user. Unfortunately, for k > 20 the output files 
generated by this program become too large. Then we require a code that 
only evaluates the coefficients Dj. and E^ for a particular value of k, which 
is accomplished by the second code mathpm. For much larger values of k, 
say for A;>100, storing the coefficients is no longer a viable option. In these 
cases the output files need to be divided into smaller files. Then each file can 
be evaluated separately in Mathematica and the result stored. Once a result 
is stored, the file can be deleted and another file can be imported into the 
software package. Once it is evaluated and the result stored, it too can be 
deleted. Finally, the stored results can be summed to yield the value for the 
coefficient. 

As a result of the success in developing a programming methodology for 
the partition method for a power series expansion, we have been able to create 
programs that can determine various types of integer partitions such as those 
with either a fixed number of elements or specific elements, doubly-restricted 
partitions and discrete/distinct partitions. Normally, different programming 
approaches are required to solve each of these problems, but as explained 
in Sec. 5, they can all be solved by introducing minor modifications to the 
program partgen in Sec. 2. In the process new operators such as the discrete 
partition operator -Ldp,A:['] and the odd- and even-element partition opera- 
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tors, LoEP,kV\ and LEEP,k[-] are defined. In particular, the number of discrete 
partitions summing to k or q{k, 1) is equal to L£)pfc[l]. Another interesting 
application in this section is the development of the program transp, which 
determines conjugate partitions by means of Ferrers diagrams. These are 
created by the dynamic memory allocation of two-dimensional arrays in the 
C/C++ code appearing in the appendix. 

In Sees. 6-8 the operator approach of Sees. 3 and 5 is employed in the 
derivation of new power series expansions or generating functions for nu- 
merous infinite products of increasing complexity that arise in the theory of 
partitions. Sec. 6 begins by studying the product P{z) defined by Equiva- 
lence (1731) . which produces a generating function or power series expansion 
whose coefficients are the partition function or p{k). Theorem 2 shows that 
the generating function of this important product is absolutely convergent 
for |2r| < 1 and divergent elsewhere. That is, unlike the geometric series, 
there is no region in the complex plane where the generating function is con- 
ditionally convergent. Instead, there is a ring of singularity separating the 
absolutely convergent unit disk from the rest of the divergent complex plane. 
As a result of Theorem 2, Equivalence (1751) is only an equation when \z\ < 1. 
For these values of z, the generating function can be inverted and Theorem 
1 can then be applied to the ensuing result. The coefficients of the resulting 
power series expansion or generating function are referred to as the discrete 
partition numbers q{k), which equal (—1)-' when k is a pentagonal number 
or equal to (3j^ ±j)/2 and zero, otherwise. We also find that they can be 
expressed in terms of the partition operator with each element i assigned to 
the partition function value oi p{i) as given by Eq. fll62p . On the other hand, 
when Theorem 1 is applied to the inverse of P{z), we find that the discrete 
partition numbers can be expressed in terms of the partition operator acting 
with each element i assigned to a value of p(z) as given by Eq. fll73p . 

Sec. 6 continues with the derivation of alternative representations for the 
generating functions of P{z) and its inverse, which are given by Eqs. (11691) 
and (11671) . respectively. In these results the coefficients for each power j 
of the inner series are expressed in terms of a sum over the divisors d of j 
divided by j and are denoted by jj. When Theorem 1 is applied to these 
new forms for the products, the coefficients of the generating functions now 
have the partition operator acting with the elements i are assigned to values 
of 7j as in Eqs. (I170p and (11750 . The difference between these results is the 
appearance of the phase factor (—1)^'= in the case of the discrete partition 
numbers. 
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Because a result like Eq. (11731) gives the partition numbers in terms of 
the partition operator acting with the elements i assigned to the discrete 
partition numbers or q{i), it becomes necessary to develop a program that 
excludes those partitions in which the q{i) vanish. This represents a com- 
pletely different type or class of partition studied in Sec. 5. Thus, Sec. 6 
presents the program called partfn, which describes the modifications that 
need to be made to partgen in Sec. 2 so that only those partitions in which 
all the elements are pentagonal numbers are printed out. 

Although the programs in Sec. 4 were developed with the partition oper- 
ator acting on all partitions, it was stated at the end of Sec. 5 that they could 
be adapted to handle situations where only a subset of the total number of 
partitions is required. That is, the programming methodology in Sec. 4 is 
not restricted to the partition operator, but as a result of the material in 
Sec. 5, it can be adapted to handle situations in which the coefficients are 
expressed in terms of different operators or specific types of partitions. Such 
a situation occurs with Eq. (11731) . Consequently, Sec. 6 concludes by pre- 
senting the program pfn, which determines the contributions made by the 
partitions whose elements yield non-zero discrete partition numbers. This 
program expresses the partition function or p{k) in two symbolic forms, both 
of which must be imported into Mathematica to obtain the actual values for 
p{k). In the first form the partition function is expressed directly in terms of 
the discrete partition numbers or q{i). This means that in order to obtain the 
values for the partition function, another module for calculating the discrete 
partition numbers must be created in Mathematica. In the second form the 
non-zero values of the q{i) are replaced by their symbolic form of (-l)^j for 
z = (3j^ ± j)/2. When the second form is entered as input into Mathematica, 
it immediately computes the partition function in integer form. 

Sec. 7 begins with the presentation of Theorem 3, which generalises the 
infinite product of 1/P{z). Instead of the coefficients of the powers of z be- 
ing equal to -1, they are now assumed to be equal to general values Ck as in 
Eq. f ll79p . The theorem is proved by adapting the proof of Theorem 1 and 
means that the coefficients of the generating series for this product are given 
in terms of the discrete partition operator acting with each element i assigned 
the value of Ci as in Eq. fllSip . Conversely, the theorem implies that any 
power series expansion for a function can be expressed as an infinite product. 
After some elementary examples are studied, viz. the geometric series and 
the exponential power series, the Ci are set equal to u, whereupon we find 
that the coefficients of the resulting generating function become the polyno- 
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mials q{k,u), which are referred to as the discrete partition polynomials. As 
expected, for u; = — 1 they reduce to the discrete partition numbers, while for 
cu = 1, they yield the number of discrete partitions summing to k. In fact, 
the coefficients of the discrete partition polynomials represent the number of 
discrete partitions where the number of elements is equal to the power of u. 
Then a brief description appears on how the program dispart presented in 
Sec. 5 can be adapted to evaluate these numbers, which can also be written 

as LDP,k,i[^]- 

The infinite product of Theorem 3 is further generalised by the intro- 
duction of an arbitrary power pk as given by Eq. ( 121 ip in the corollary to 
Theorem 3. In this case the coefficients of the resulting generating function 
are expressed in terms of the partition operator acting with the elements i 
assigned a value of —Ci multiplied by (— Pi)n,, where rii is the number of 
occurrence of rii and {p)k denotes the Pochhammer symbol. With the aid of 
other results appearing in the corollary, we are able to study the generating 
functions when 1/P{z) is squared or cubed. For these cases pi has been set 
equal to the uniform value of p and the discrete partition polynomials q{k, u) 
are extended to become q{k,u,p) as defined by Eq. (12221) . Several identities 
involving the q{k,oj,p) are also derived. 

Sec. 8 is also devoted to deriving generating functions for other forms of 

infinite products using the material from the previous two sections. First of 

all, the coefficients of the powers of z in P{z) are set equal to — w instead 

of -1. Consequently, the coefficients of the generating function obtained 

after applying Theorem 1 to the modified version of P{z) or P{z,u) become 

polynomials in u, which are given by the partition operator acting with the 

elements i assigned to the value of u as in Eq. ( I244p . As expected, these 

polynomials denoted by p{k, u) reduce to the partition function p{k) for 

u = l. Their coefficients are referred to as the sub-partition numbers Pi{k), 

k 
but are also equal to , where the latter notation was introduced in Sec. 2 

m 

to denote the number of partitions summing to k with m elements. 

Sec. 8 continues with further generalisations of P{z), where the coeffi- 
cients of the powers of z are set equal to Ck- This represents not only the 
inversion of Theorem 2, but is also a special case of the corollary to Theorem 
3. The specific case of the Ck equalling unity has the interesting property 
that the coefficients of the generating function represent the difference be- 
tween the number of even- and odd-element partitions summing to k. The 
absolute values of these coefficients are shown to equal the number of dis- 

135 



Crete partitions with only odd elements, a result first obtained by Euler. 
Then recurrence relations are developed for the number of discrete partitions 
summing to 2k and 2k + 1 given by Eqs. fl274p and fl275p . respectively. Like 
the Euler/MacMahon recurrence relation or Eq. (I165p . which evaluates the 
values oi p{k), the new recurrence relations also require the discrete parti- 
tion numbers. This means that not many of the previous values of q{k, 1) 
are required to determine the highest successive value. 

The approach that resulted in the alternative representations for P{z) 
and its inverse given by Eqs. fll67p and fll69p is then apphed to P{z,u). 
Replacing the 7^ are the divisor polynomials 7i(co'), where each coefficient 
represents a divisor d oii divided by i while the power of u equals the inverse 
or reciprocal of this value. As a result, the discrete partition and partition 
function polynomials can be expressed in terms of the partition operator 
acting with the elements i assigned to —•ji^u) and 7i(co'), respectively. Next, 
the program dispfnpoly is presented. The purpose of this code is to express 
both q{k, —u) and p{k, u) in symbolic form so that they can be imported 
into Mathematica, where the Divisors routine can be exploited to yield the 
final forms as polynomials in u of degree k. 

Sec. 8 concludes by applying the the preceding material to derive the 
generating functions for more advanced infinite products. First, the product 
Pp{z,oj) is studied, in which an arbitrary power of p is applied to P{z,ijj). 
The generating function for the new product has coefficients p{k, u, p), which 
are expressed in terms of the partition operator acting on elements i assigned 
values of p{i,u), but now the multinomial factor is altered as described in 
the proof of Corollary 1 to Theorem 1. The next example is the product 
of Q{z, —(3u) as defined by Eq. (11980 and P{z,auj). The coefficients of the 
generating function become special polynomials QPk{uj, (3,a) as defined by 
Eq. (12861) and are tabulated in Table O In these coeffcients the power of u 
represents the number of elements in the final partitions, which are composed 
of the elements from both discrete and standard partitions, while the powers 
of a and (3 represent the numbers of elements emanating from the standard 
and discrete partitions, respectively. The final example in the section is the 
derivation of the generating function for the famous three-parameter plus 
one variable product first studied by Heine and given by Eq. (12910 . This 
product, which represents the product of P{z,u!x,u!) and P{z,uy,uxy), is 
found to possess a generating function whose coefficients are polynomials 
denoted by H Pk^uj , x , y) . Like the QPk{uj, (3,a), they are of degree k in u. 
General expressions for some of the coefficients are derived, while the first 
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seven polynomials are presented in Table El 
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11 Appendix 

The first code presented in this appendix is the program partmeth discussed 
in detail at the beginning of Sec. 4. This code employs the BRCP alogrithm 
of Sec. 2 to generate the coefficients arising from the partition method for 
a power series expansion according to Theorem 1 in symbolic form so that 
the values of the coefficients can be evaluated either in integer form or as 
algebraic expressions in Mathematica. 

/* This code deals with the application of the p artition 
method for a power series expansion to the pseudo- 
composite function g ( af ( z ) ) . Here it is assumed that 
g(z)= h(z)(l+ q_l z + q_2 z^2+ . . . + q^k z^k +...), where 
h(z) can he any function , hut is usually equal to unity 
or some factor multiplied hy a non— integer power of z. 
In addition , f(z) is assumed to he a power series 
expansion of the form (p-0 + p_l y+ p_2 y^2 + ...+ 
P-k y^k + ...) with y=z^ alpha . This code is only valid 
for the case of p-0=0. The co efficients DS[k,n] and 
ES[k,n] are those in Theorem 1 and are computed in a 
format suitahle for processing in Mathematica . */ 

T^include <stdio.h> 

T^include <memory.h> 

T^include <stdlib.h> 

T^include <math . h> 

T^include <time.h> 

int dim , sum,* part , inv_case ; 
time_t init_time , end_time ; 

/* // inv-case equals unity, then the code will evaluate 
the coefficients of the inverse power series, i.e. 
the power series for h( af ( z) ) / g ( af ( z ) ) . */ 

void termgen(int p) 

{ 

int freq , i , num_p arts =0,1 ,spacing=0,num_dis_parts=0,dis_parts=0; 
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double sign , dnum_parts ; 

/* num_parts is the total # of parts in the partition , 

while num_dis_parts is the number of distinct parts 

with greater than unity frequency snd is required 
for the multinomial factor. */ 

if(j^^um) printf ((inv_case >0)?"ES[%i , n_] := ":"DS[%i,n. 
if ( inv_case==0 &fe p!=sum) printf ("+"); 
for(i=l; i<=dim; i++){ 
freq=part [ i ] ; 
if ( freq ) { 

dis_parts++; 

if ( inv_case==0){ 

printf("p[%i ,n]" ,i); 

if(freq>l) pr int f ( " "(% i ) ",freq); 

else printf (" "); 

} 
num_dis_parts += (freq>l); 

num_parts=num_parts+freq ; 

} 

} 
if ( inv_case >0){ 

dnuni_parts = (double) num_parts ; 
sign=pow( — 1.0 , dnum_parts ) ; 
printf ((sign >0.0)?"+" : "-" ); 

} 



if ( inv_case >0){ 

printf("DS[0,0]"(-%i) " , num_parts + 1); 

} 
else { 

print f ("q[% i ] a" , num_parts ) ; 

printf (( num_parts >1)? ""(%i) " : " " , num_parts ) ; 

} 

if ( inv_case >0){ 
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:")"); 



for ( i =1; i<=dim ; i++){ 
freq=part [ i ] ; 

if(freq==l) printf ("DS[%i ,n] ",i); 
else if(freq>l) printf ("DS[%i , n] ^ (% i ) ",i,freq); 

} 

} 
if ( num_parts>l &fe dis_parts >1){ 

pr int f ("%i ! " , num_parts ) ; 
if(num_dis_parts){ 

printf ((num_dis_parts>l)? "/(" : "/"); 
for(i=l; i<=dim; i++){ 
freq=part [ i ] ; 
if (freq >1){ 

if ( spacing++) printf (" "); 
printf ("%i!" ,freq); 

} 

} 
if ( num_dis_parts >1) printf 

} 
printf (" 

} 

printf("\n"); 

} 

void idx(int p , int q) 

{ 

part [p] + + ; 

termgen (p ) ; 

part [p]--; 

p -= q; 

while (p >= q) { 

part [q] + + ; 

idx(p ,q); 

part [q++] ; 

} 
} 

int niain(int argc , char *argv 
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" "); 



{ 

int i ; 

double delt a_t ; 

time(&init _t ime ) ; 

if(argc != 2) printf (" usage : ./partmeth <#partitions >\n" 

else { 

dim=atoi ( argv [ 1 ] ) ; 
for(sum=l; sum<=dini ; sum++){ 

part=(int *) malloc ( (dim + l)* sizeof ( int ) ) ; 

inv_case =0; 

for(i=l; i<=dim; i++) part[i] = 0; 

idx (sum , 1 ) ; 

free ( part ) ; 

printf("\n"); 

inv_case =1; 

part = (int *) malloc ( (dim + 1)* sizeof ( int ) ) ; 

for ( i =1; i<=dim ; i++) part[i] = 0; 

idx (sum , 1 ) ; 



free ( part ) ; 
printf("\n"); 



} 

} 

printf("\n"); 

time(&;end_time ) ; 

delta_t= difftime( end_time ,init_time); 

printf (" Computation time is %f seconds\n" , delt a_t ) ; 

return (0) ; 

} 

The above code represents the implementation of Theorem 1 for the case 
where po vanishes. Once the order of coefficients becomes sufficiently large, 
viz. for k > 20, the code needs to be adapted so that only specific values 
for one of the two types of coefficient are evaluated in symbolic form. This 
requires: (1) separating the inverse case or the Ek from the Dk, and (2) 
removing the first for loop in main so that only i = dim is computed. The 
modified code called mathpm, which determines only the D^. in symbolic 
form, is presented below. 
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/* The code mathpm determines the c o efficients of the 
power series for the pseudo— composite function 
g(af(z)), where g(z)= h(z) (1+ q^l z + q-2 z^2+ ... 
+ q-k z^k +...) and h(z) can he any function . The 
function f(z) must he expressed as (p-0 + p-1 y+ 
p_2 y^2 + ... + p_k y^k + ...) where y=z^ alpha and 
p_0 =0. */ 

^include <stdio.h> 
T^include <memory.h> 
T^include <stdlib.h> 
9^include <time.h> 

int dim ,* part ; 

long unsigned int term = l; 

time_t init_time , end_time ; 

void termgen(int p) 

{ 

int f , i , num_parts=0, dis_part_cnt=0,l , num_dis_parts=0; 

/* num_parts is the total # of parts in the partition , 

while num_dis_parts is the numher of distinct parts. 

The latter is required for the multinomial factor. *^ 

if(j^=dim) printf("DS[%i ,n_]:= ",p); 
else { 

printf ("+ " ); 

if (term %3 ^ 0) printf (" \n" ) ; 

term++; 

} 
for(i=l; i<=dim; i++){ 

f=part [ i ] ; 

if(f){ 

printf("p[%i]" ,i); 

if(f>l) printf(""(%i) " , f); 

else printf(" "); 

num_parts += 1 = f ; 

num_dis_parts += (f>l); 
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} 

} 

print f ("q[% i ] a" , num_parts ) ; 

printf ( ( num_parts >1)? ""(%i) " : " " , num_parts ) ; 
if( num_parts > 1 ){ 

printf ("%i ! " , nuni_parts ) ; 
if(num_dis_parts){ 

printf ((num_dis_parts>l)? "/(" ■ "/"); 
for(i=l; i<=dim; i++){ 
f=part [ i ] ; 
if(f>l){ 

if ( dis_part _cnt++) printf(" 
/* dis-part-cnt counts the # of distinct parts with 
greater than unity frequency It is required to 
insert blank spaces in the denominator of the 
multinomial factor */ 

printf ("%i!" , f ); 

} 

} 

if ( nuni_dis_parts >1) printf(" 

} 

printf(" "); 

} 

} 

void idx(int p , int q) 

{ 

part [p] + + ; 

termgen (p ) ; 

part [p]--; 

p -= q; 

while (p >= q) { 

part [q] + + ; 

idx(p ,q); 

part [q++] — ; 

} 
} 
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"^"); 



int main( int argc , char *argv[]) 

{ 

int i ; 

double delt a_t ; 

FILE *ptr ; 

char filename [10] = " times" ; 

time(&init _t ime ) ; 

if (argc != 2) printf (" usage : . /mathpm <7^partitions >\n" ); 

else { 

dim=atoi ( argv [ 1 ] ) ; 

part=(int *) malloc ( (dim + 1)* sizeof ( int ) ) ; 

if ( part==NULL) printf (" unable to allocate array\n\n"); 

else { 

free ( part ) ; 

} 

} 
printf ("\n" ); 

time(&end_time ) ; 

delta_t= difftime( end _t ime , init_time ) ; 

ptr=fopen (filename , " a" ) ; 

fpr int f ( ptr , "Time to compute p[%i,n] is %f seconds\n", 

dim , delt a_t ) ; 

f close ( ptr ) ; 

In Sec. 5 we discussed the problem of generating discrete partitions or 
partitions in which the elements only appear once. Since discrete partitions 
are important in the theory of partitions as can be seen from Sees. 6-8, the 
entire code called dispart. cpp is presented below. 

T^include <stdio.h> 
T^include <memory . h> 
T^include <stdlib.h> 
:^include <time.h> 

int tot , *part ; 

long unsigned int term = l; 

time_t init_time , end_time ; 
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void termgen () 

{ 

int freq , i ; 

for (i=0;i<tot ; i++){ 

freq=part [ i ] ; 
if(freq>l) goto end; 

} 
printf ("%ld: " ,term + + ); 

for(i=0;i<tot ; i++){ 

freq=part [ i ] ; 

if (freq) pr int f ( "%i (%i ) ",freq,i+l); 

} 
printf ("\n" ); 



end : : 
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int main( int argc , char *argv 

{ 

int i ; 

double delt a_t ; 

FILE >Kptr ; 

char filename [10] = " times" ; 

time(&init _t ime ) ; 
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if(argc != 2) print f (" usage : ./dispart <#partitions >\n" ) ; 
else { 

tot=atoi ( argv [ 1 ] ) ; 

part= (int * ) malloc ( tot * sizeof ( int ) ) ; 

if (part ^= NULL) printf (" unable to allocate array\n\n"); 

else { 

for ( i=0;i<tot ; i++) part[i] = 0; 
idx ( tot , 1 ) ; 
free ( part ) ; 

} 

} 
printf("\n"); 

tinie(&end_tinie ) ; 

delta_t = difftinie (end_time , init_tinie ) ; 

ptr=fopen ( filename , " a" ) ; 

fpr int f ( ptr , "Time taken to compute discrete partitions 

summing to %i is %f sees . \n" , tot , delta_t ) ; 

f close ( ptr ) ; 

return (0) ; 

The next code presented here determines the conjugate partition as de- 
scribed in Sec. 5. In evaluating conjugate partitions, the programme called 
transp casts the original partition in the form of a Ferrers diagram, but, 
instead of being composed of dots, the Ferrers diagram is composed of ones. 
Thus, the conjugate partition is determined by summing the ones in each 
column of the Ferrers diagram. 

/* This program evaluates the partitions and their 
conjugates for any integer greater than or equal to 
2. Conjugates are determined by summing the columns 
in the Ferrers diagram for each partition. */ 

T^include <stdio.h> 
T^include <memory.h> 
T^include <stdlib.h> 

int tot , *part ; 

long unsigned int term = l; 
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void termgen () 

{ 

int freq ,i ,j ,k,next_el , index, prev_el ,*part2 , 

*sum_col ,row_cnt=0; 
int * f errers ,** rptr ; 
/* The Ferrers array and its array of pointers */ 

printf (" Partition %ld is: ",term++); 
ferrers=(int *) malloc ( tot *tot *sizeof ( int ) ) ; 
rptr=(int **) malloc ( tot *sizeof( int *)); 

/* Get the pointers to the rows of ferrers */ 
for ( i =0; i<tot ; i++) rptr[i]= ferrers+(i *tot ) ; 
/* Here tot refers to number of columns */ 

for(j=0;j<tot;j++){ 

for (i=0;i<tot ; i++){ 
rptr [j ] [ i] = 0; 

} 

} 

/* Creation of the Ferrers diagram. Rather than 

being composed of dots the Ferrers diagram is 
composed of unit values since these will be 
used to determine the conjugate p artition . */ 

for(j=0;j<tot;j++){ 
freq=part [ j ] ; 
if ( freq ) { 

for ( i=row_cnt ; i<row_cnt+freq ; i++){ 
for(k = 0;k<=j ;k++){ 
rptr [ i ] [ k] = 1 ; 

} 

} 

} 

row_cnt=row_cnt+freq ; 

} 

/* Summation of the columns in the Ferrers 

diagram yielding a new array called sum_col . *, 
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sum_col = (int *) malloc ( tot *sizeof ( int ) ) ; 

/* Initialising the array elements to zero */ 

for ( i =0; i<tot ; i++) sum_col [ i ] =0; 

/* Summation of the columns now occurs */ 

for (j=0;j<tot;j++){ 

for(i=0;i<tot; i++) sum_col [j] = sum_col [j] + rptr [ i ] [ j ] ; 

} 

/* The array sum_col is reduced to the conjugate 
of the original partition through another array 
called part2 */ 

part2 = (int *) malloc ( tot * sizeof ( int ) ) ; 

/* Initialisation of the array elements to zero */ 

for ( i=0;i<tot ; i++) part2[i]=0; 

/* Now the conjugate p artition is arranged as 
in the order of part */ 

prev_el=sum_col [0] ; 

/* the highest part is the value of sum^col [0] */ 

index=prev_el —1; 

/* the index in the p artition array is one less */ 

part2 [ index ] = 1 ; 

for(i=l;i<=tot ; i++){ 

next_el=sum_col [ i ] ; 

if ( next_el==0) goto out; 

if ( next _el^=prev_el ) part2 [ index ] = part2 [ index] + 1; 

else { 

index=next _el —1; 
part2 [ index ] = 1 ; 
prev_el=next _el ; 

} 

} 
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out: for ( i=0;i<tot ; i++){ 
freq=part [ i ] ; 
if(freq) print f ("%i(%i ) ",freq,i+l); 

} 
/* The conjugate is actually the reverse order 

of part2 */ 

printf(" and its conjugate is: "); 

for( i=0;i<tot ; i++){ 

freq=part2 [tot— i —1]; 

if(freq) print f ("%i(%i ) " , freq , tot-i ) ; 

} 

free ( sum_col ) ; 
free ( part2 ) ; 
free ( f errers ) ; 
free ( rptr ) ; 
printf("\n"); 

} 

void idx(int p,int q) 

{ 

part [p-l] + + ; 

termgen ( ) ; 

part [p — 1] ; 

p -= q; 

while (p >= q) { 

part [q-l] + + ; 

idx(p ,q); 

part [q++ -1] ; 

} 

} 

int niain(int argc,char *argv[]) 

{ 

int i ; 

if(argc != 2) printf (" usage :./ transp <#partitions >\n" 

else { 

tot=atoi ( argv [ 1 ] ) ; 

part = (int *) malloc ( tot *sizeof ( int ) ) ; 
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if (part ^= NULL) printf (" unable to allocate 

array \n\n" ) ; 
else { 

for ( i=0;i<tot ; i++) part[i]=0; 

idx ( tot , 1 ) ; 

free ( part ) ; 

} 

} 

printf ("\n" ); 

return (0) ; 

} 

In Sec. 6 the partition method for a power series expansion was applied to 
an exponential form of the generating function of the partition function p{k), 
which yielded Eq. fll73p . Although this result represents a sum over partitions 
involving the special numbers q{i) called the discrete partition numbers, it 
incorporates much redundancy because these numbers are often zero except 
when i can be written as a pentagonal number or as (3j^ ± j)/2, where j is a 
non- negative integer. Then they are equal to (— !)■'. Appearing below is the 
program pfn which gives the partition function based on the properties of the 
discrete partition numbers (represented by Q[i\ in the code) in symbolic form. 
The actual values of the partition function can be evaluated by importing 
the output into Mathematica. 

T^include <stdio.h> 
7^1 n cl u d e <memory . h> 
T^include <stdlib.h> 
T^include <time.h> 
T^include <math . h> 

int dim ,* part , limit ,freq , first_term=0; 
long unsigned int term = l; 
time_t init_time , end_time ; 

void termgen(int p) 

{ 

int f , i , num_parts=0,j ,jval , dis_part_cnt=0,l , 

num_dis_parts=0; 

/* num_parts is the total # of parts in the partition 
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while num_dis_parts is the number of distinct parts. 
The latter is required for the multinomial factor. */ 
/* jv al =0; * / 

if (p^=dim) { 

printf("p[%i]:= " ,p); 
for ( i =1; i<=dim ; i++){ 
jval=0; 
freq=part [ i ] ; 
if (freq >0){ 

for (j=l;j<=limit ; j++){ 

i f ( i = (3*j— l)*j/2) jval=j; 
if(i = (3*j+l)*j/2) jval=j; 

} 
if((jval==0) && (freq>0)) goto end; 

first _term++; 

} 

} 

} 
else { 

for ( i =1; i<=dim ; i++){ 

j val =0; 

freq=part [ i ] ; 

if (freq >0){ 

for (j=l;j<=limit ; j++){ 

if(i = (3*j — l)*j /2) jval=j; 

if(i = (3*j+l)*j/2) jval=j; 

} 
if((jval==0) && (freq>0)) goto end; 

} 

} 

if (first_term ! = 0) printf("+"); 
if ( first -term ==0) f irst _t erm++; 

if (term %3 ^ 0) printf (" \n" ) ; 
term++; 

} 
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for(i=l; i<=dim; i++){ 
f=part [ i ] ; 
if(f ) { 

/* przntf("Q[%zJ",z); */ 

for (j=l;j<=limit ; j++){ 

if (i = (3*j-l)*j/2) jval=j ; 
if (i = (3*j+l)*j/2) jval=j; 

} 

printf("((-l)"(%i))"Jval); 



/oi 



f 



if(f>l) printf 
else printf (" "); 
num_parts += 1 = f ; 
num_dis_parts += (f>l); 

} 

} 

printf("(-l)" ); 

print f (( num_parts >1)? ""(%i) " : " " , nuni_parts ) ; 
if( num_parts > 1 ){ 

print f ("%i ! " , num_parts ) ; 
if(num_dis_parts){ 

printf ((num_dis_parts>l)? "/(" : "/"); 
for(i=l; i<=dim; i++){ 
f=part [ i ] ; 
if(f>l){ 

if ( dis_part _cnt++) printf(" "); 
/* dis_part_cnt counts the # of distinct parts with greater 
than unity frequency . It is required for inserting blank 
spaces in the denominator of the multinomial factor */ 

printf("%i!" , f ); 
} 

} 
if ( nuni_dis_parts >1) printf(")"); 

} 
printf (" " ); 

} 
end : ; 

} 
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void idx(int p , int q) 

{ 

part [p] + + ; 

termgen (p ) ; 

part [p]--; 

p -= q; 

while (p >= q) { 

part [q] + + ; 

idx(p ,q); 

part [q++] ; 

} 

} 

int niain( int argc , char *argv [] ) 

{ 

int i ; 

double delt a_t ; 

time(& init _t ime ) ; 

if(argc != 2) print f (" usage : ./pfn <#partitions >\n" ); 

else{ 

dim=atoi ( argv [ 1 ] ) ; 

limit= floor (1+sqrt ((1+ 24 * dim))/6); 

part = (int *) malloc ( (dim + l)*sizeof ( int ) ) ; 
if ( part^=NULL) printf (" unable to allocate array\n\n' 
else { 
/* for(sum = l; sum <=dim ; sum++){ */ 

for (1=1; i<=dim; i++) part [ i ] = 0; 
idx (dim , 1 ) ; 
/* printf ("\n" ) ; 

} */ 

free ( part ) ; 

} 

} 

printf("\n"); 

time(&end_time ) ; 

delta_t= difftimef end_time , init_time ) ; 
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printf ("Time taken to compute the coefficient is 

%f seconds\n" , delt a_t ) ; 
return (0) ; 

} 

Also in Sec. 6, the discrete partition polynomials q{k,uj) were found to 
be the coefficients of the generating function for by the product Q{z,u) = 
rifcli 1/(1 + 1^-2^) (see Eq. 01961) ). Later, it was found that these polynomials 
of degree i in co could be expressed in terms of the partition operator acting 
with the elements i equal to special polynomials 7j(w). The latter were 
referred to as divisor polynomials since their coefficients are divisors or factors 
of i. The relationship between both types of polynomials is given by Eq. 
), while another result below it relates the partition function polynomials 



p{k, u) to another sum over partitions acting on the divisor polynomials. 
Appearing below is the program called dispfnpoly, which prints out both 
the discrete partition and partition function polynomials for a specified order 
in symbolic form so that they can be handled by Mathematica ^12j. 

T^include <stdio.h> 
T^include <memory . h> 
T^include <stdlib.h> 
T^include <time.h> 

int dim ,* part , polytype ; 
long unsigned int term = l; 
time_t init_time , end_time ; 

void termgen(int p) 

{ 

int f ,i ,num_parts=0,dis_part_cnt=0,l ,num_dis_parts=0; 

/* num_parts is the total # of parts in the partition 

while num_dis-parts is the number of distinct parts. 

The latter is required for the multinomial factor */ 

if(j^=dim) printf ((polytype==l)?"Q[%i,-w_]:= ": 
"P[%i,w_]:= ",p); 

for(i=l; i<=dim; i++){ 
f=part [ i ] ; 
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if(f){ 



} 



num_parts += 1 = f ; 
num_dis_parts += (f>l); 

} 



,f); 



if (term ! = 1) printf (" + " ); 
if (term % 3 ==0) printf (" \n" ) ; 
term++; 

for(i=l; i<=dim; i++){ 
f=part [ i ] ; 
if(f){ 

printf ("DP[%i ,w]" ,i); 
if(f>l) printf(""(%i) 
else printf(" "); 

} 

} 

printf ((polytype==l)?"(-l)" :"" ); 

if (num_parts>l) print f (( polytype==l)?" "(%i )":"" , 

num_parts ) ; 
if(num_parts > 1 ){ 

if(num_dis_parts){ 

printf ((num_dis_parts>l)? "/(" : "/"); 
for(i=l; i<=dim; i++){ 
f=part [ i ] ; 
if(f>l){ 

if ( dis_part _cnt ++) printf (" "): 
/* dis_part_cnt counts the # of distinct parts with greater 
than unity frequency . It is required for inserting blank 
spaces in the denominator of the multinomial factor */ 

printf("%i!" ,f ); 
} 

} 
if(num_dis_parts>l) printf(")" ); 

} 
} 

} 



void idx(int p,int q 



158 



{ 

part [p] + + ; 
termgen (p ) ; 
part [p]--; 
p -= q; 
while (p >= q) { 

part [q] + + ; 

idx(p ,q); 

part [q++] ; 

} 
} 

int main(int argc , char *argv[]) 

{ 

int i ; 

double delt a_t ; 

time(&init _t ime ) ; 

if(argc != 2) print f (" usage : . / dispfnpoly 

<#partitions >\n" ) ; 
else{ 

dim=atoi ( argv [ 1 ] ) ; 

part=(int *) malloc ( (dim + l)*sizeof ( int ) ) ; 

if ( part^=NULL) printf (" unable to allocate 

array \n\n" ) ; 
else { 

for(polytype=l; polytype <=2;polytype++){ 
for (1=1; i<=dim; i++) part [ i ] = 0; 
idx (dim , 1 ) ; 

if ( polytype==l) printf (" \n\n" ); 
term = l; 

} 
free ( part ) ; 

} 

} 

printf("\n"); 

time(&end_time ) ; 

delta_t= difftime( end_time , init_time ) ; 
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print f ("Time taken to compute the coefficients is 

%f seconds\n" , delt a_t ) ; 
return (0) ; 

} 
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